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L. Hyperfine Structure in Paramagnetic Salis and Nuclear Alignment. 


By B. Buranery, F.B.S., 
Clarendon Laboratory, Oxford *. 


[Received February 14, 1951.] 


ABSTRACT. 


Abragam and Pryce (1951) have shown that the hyperfine structure in 
paramagnetic salts can be interpreted using a simple Hamiltonian involving 
components of the electron and nuclear spins. .On the basis of this 
Hamiltonian, general formule are given for the allowed transitions in 
strong magnetic fields, assuming that the paramagnetic ions have axial 
symmetry and that interactions with the electric quadrupole moment 
of the nucleus are small compared with those involving the magnetic 
dipole moment. 

As far as possible, general formule are also given for the energy levels 
in zero magnetic field, and the allowed transitions are indicated. The 
use of hyperfine structure in paramagnetic salts for the production of 
nuclear alignment is discussed, with some quantitative data. 


§1. INTRODUCTION. 


In a general survey of the theory of paramagnetic resonance and 
hyperfine structure, Abragam and Pryce (1951) have shown that it is 
possible to represent the behaviour of the energy levels by the following 
Hamiltonian, if the anisotropy has axial symmetry about the z-axis: 

H=B(g HS,+9,(H,8,+H,8,) }+D{$7—-38{8+])} 

+A8,1,+ B(8,1,+8,1,)+ Q{Z—gl0+1)}—yByH.I, - (1) 
B=Bohr magneton, 

By=Nuclear magneton. 
Here the effective electron spin S is defined by setting the multiplicity of 
the electronic energy levels equal to 2S+1, and I is the nuclear spin. 
g is a tensor, with principal values g, and g,. The term in D represents 
the presence of initial splittings of the electronic levels due to second 
order effects of the crystalline electric field. The terms in A, B correspond 
to the interaction between the nuclear magnetic moment and the magnetic 
field of the unfilled electron shell. The term in Qj is due to the interaction 
of the nuclear electric quadrupole moment with the gradient of the 
electric field at the nucleus, while the last term takes account of the 
direct effect of the external magnetic field on the nuclear magnetic moment. 


* Communicated by the Author. 
+ Unfortunately this nomenclature differs from that customary where Q is 


used for the actual quadrupole moment of the nucleus. The relation between 
them is that Q (this paper) =3eQ @?V/dz?/41(2I—1) (usual nomenclature). 
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This Hamiltonian is very convenient for the interpretation of 
paramagnetic resonance spectra, since it involves only the electronic 
and nuclear spins, together with a set of parameters representing the 
strength of the various interactions. These can be regarded as a set of 
unknowns to be determined experimentally by a comparison of the 
observed spectrum with that calculated from the Hamiltonian given 
above. This spectrum is normally relatively simple, but the anisotropy 
introduces a complicated variation with the direction of the external 
magnetic field. The computed variation has already been compared 
with experiments in a number of specific cases and good agreement found. 

The purpose of this paper is threefold. First, to give general formule 
for the allowed transitions in strong magnetic field, using perturbation 
theory carried to the second order. This is only possible with two 
limitations ; (a) the axial symmetry presumed in equation (1); (5) the 
magnitude of the quadrupole parameter Q must be small compared with 
A, B. Second, to give, as far as possible, formule for the energy levels 
in zero magnetic field, together with the allowed transitions. These 
transitions are important because they can give more accurate values of 
the parameters (in particular, the quadrupole parameter Q) than the 
strong field case. The energy levels of zero fields are important also 
because of their effects in experiments on adiabatic demagnetization. 
Third, to consider the use of hyperfine structure in producing alignment 
of the nuclear spin, as suggested by Gorter (1948) and Bleaney (1951a). 


§2. THe Srrone Frecp SPEcTRUM. 


In strong fields it is convenient to choose the axes so that elements in 
S,,, 8, do not occur in the major term g$ H..S in the Hamiltonian. When 
g is anisotropic this is equivalent to taking the axis about which the spin 
precesses as the new z-axis. Then, if the external field is applied at 
an angle @ to the axis of crystalline symmetry, the usual rules for rotation 
of the axes may be applied to the Hamiltonian if allowance is made for 
the non-commuting properties of the spin components. We will consider 
first the fine structure which arises in the absence of nuclear spin from the 
presence of a crystalline splitting of the electronic levels. 


2.1. The Fine Structure. 


In strong fields the allowed transitions are of the type 4M=+1. For 
the transition M+M—1, we have 


hv=gBH+D(M—4) {3 7 008? — oa 


Dgig, cos 6 sin 0\2 
D 2 ay 26 2 1 
+( a Waivers M(M~—1)=3}, 0 20. 22) 


where P=9, cos? 6+-g* sin? 0. JER ANAT AT Sy aS) 
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The physical basis of (3) can be seen from fig. 1. The magnetic field H 
has components H cos @ and H sin @ respectively along and normal to the 
symmetry axis. The components of the magnetic moment are therefore 
proportional to g,Hcos@ and gH sin@, and g is proportional to the 
length of the vector making an angle ¢ with the z-axis which is formed 
by compounding these two. This vector is parallel to the resultant 
magnetic moment of the electrons, but is not parallel to H. This is 
responsible for the factors (g,/g) and (g,/g) which appear as coefficients 
multiplying cos@ and sin@ respectively everywhere in equation (2) ; 
for cos 6=(g)/g) cos@: sin d6=(g,/g) sin 4. 


Rig, 1. 


AG 050 Kg 


41 C080 


“Hsin 8 gi sind Bg, sin 8 


. 
. 


. 
. 
. 
. 


“i (Kg) 


Oneiition of magnetic field (H), electron spin (g), and nuclear spin (Kg) when 
anistropy is present. The figure is drawn for g,>g), and B>A. The 
dotted vector (Kg) refers to the case where A is negative. 


The fine structure results in a splitting into 2S lines, which are equally 
spaced in the first approximation. The spacing varies with angle, falling 
to zero where (g,/9) cos = 1/4/(3). This equal spacing is disturbed by 
the second order terms which vanish in strong fields, or along the axis 
of symmetry (9=0). The perturbation denominator in these second order 
terms is taken as (g8H,), the average spacing of successive levels. The 
significance of H, is that in the absence of the fine structure all the lines 
would coincide and occur at this field. Writing hy=gfHp, the formula (2) 
can be re-arranged in the form 


D\? 1 
H=Hy— 3 Ot) A0-(G) gh - - - - 2a 


2H2 
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This formula shows that when the spectrum is displayed at constant 
frequency as a function of field, as is normally the case in paramagnetic 
resonance, its appearance is fundamentally the same (but reversed) as 
if it were displayed at constant field as a function of frequency. The 
value of the parameter D can be found from the separations (in field) 
between the lines. 

The intensity of these strong transitions is given by the formula 

Power absorbed by crystal=Z{S(S+1)—M(M—1)}, . . (4) 
where 
m2(gB'H')?_N 
ZL= Sree oe Vo, Av), 


N=no. of paramagnetic ions in crystal. 


F(v, vo, Av) is a “ shape factor ”’ depending on the line width. At the 
centre of the line where v=v,, FS 1/4v, where Av is the half-width at 
half-intensity. 

The formula for Z is only approximate because it has been assumed that 
the radio-frequency field H’ cos (27vt) inducing the transitions is normal 
to the axis about which the spin precesses, and g’ is the value of the 
spectroscopic splitting factor along this normal. In practice H’ is usually 
normal to H, and therefore makes an angle (@—d¢) with the assumed 
direction. The absolute intensities will therefore differ slightly from (4), 
but the relative intensities given here and later will be correct. A minor 
consequence of the anisotropy is that the intensities will not be quite 
zero if the r.f. field H’ is parallel to H ; and this can occur also through 
second order effects of the splitting in weak fields, even when g is isotropic. 

In finite fields transitions corresponding to changes in M of greater than 
unity can be observed. These are smaller in intensity than the main 
transitions by a factor of the order (D/H)?, and vanish only if the external 
field is parallel to the symmetry axis. The actual intensity formule are 
rather complicated, and show that the intensity does not vary greatly 
with the angle between the external and r.f. fields. We shall content 
ourselves with giving the formula for the positions of the lines 
corresponding to the transition MM —2 in the strong field quantum 
numbers. This is 


2 
hv=2gBH+2D(M—1) {3 MN 6632 9@— i} 


g 
— (Dg\g, 008 6 sin =) oe Biss ‘ay 
( —e (ser) {4S(S-+1)—24M(M—2)—33} 
Dg? sin? 6\2 / 1 ) if 
+( a ) (sar {28(S-+1)—6M(M—2)—9} . , . (5) 


(note that H in the perturbation denominators refers again to the mean 
field at which the 4M= +1 transitions occur). 


—— 
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2.2. The Hyperfine Structure. 


In the presence of hyperfine structure the strong allowed transitions 
are those in which 4m—0, where m is the nuclear magnetic quantum 
number. The position of the transition (M, m)+>(M—1, m) is found by 
adding to the right-hand side of equation (2) the quantity 


B2 /A?+K?2 
Km+ 4gBH, Ce) {L(I-+1)—m?} 
Beas 
1 A2— B2\2 2 
+ sor ( K y (4%) sin? 6 cos? 6 m? 
0 
Q? cos? 6 sin? 6 /AB 2 
Q? sint6@ /Bg,\4 
SRS) miOIl)—2mi—Vj, -. . - . . (6) 
where K2g2= A2g? cos? 6+ Bg? qa ee eer 


The physical significance of the relation (7), which is due to Pryce (1949), 
is similar to that of (3). The magnetic moment of the electrons has 
projections proportional to g, cos 6 and g, sin @ parallel and normal to the 
symmetry axis respectively. The components of the magnetic field acting 
on the nucleus are therefore proportional to Ag, cos@ and Bg, sin 8, 
respectively, and the nucleus is aligned parallel to their resultant (see 
fig. 1). It may be noted that in deducing these formule the axes of 
quantization of electron and nucleus are taken parallel to their respective 
magnetic moments. This is responsible for the introduction of factors 
such as (Ag,/Kg) multiplying cos 0, etc. 

In the first approximation, the result of the nuclear interaction is to 
split each electronic transition into 2I+1 equally spaced components, 
with separation K between successive lines. Each component has the 
intensity given by (4), divided by 2I+1. Where there is considerable 
anisotropy the splitting which, measured in field, equals (K/g6), may go 
through a maximum value at an angle to the axis, as in cobalt (Bleaney and 
Ingram 1951la); this is due to the different angular variation of K and g. 
It should be noted that there is no first order effect due to the quadrupole 
interaction so long as Q<B. The quadrupole interaction produces a 
first order shift of the energy levels (see fig. 2), but this is the same for all 
levels with the same value of m, and does not affect transitions where 
Am=0. This is true for all values of Q if the external magnetic field is 
parallel to the axis of symmetry, for the quadrupole term then contains 
only diagonal elements. If H is at an angle with this axis the quadrupole 
interaction will break down the ordinary selection rule when Q~K, 
because of the competition between the magnetic and electrostatic 
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interactions which are each trying to align the nucleus along its own 
preferred axis. Perturbation theory can then only be applied if Q<K, 
that is, Q<A,B. This is reflected in the fact that K appears in the 
denominator in the second order terms in Q?, which do not therefore 
vanish in strong fields. The remaining second order terms do. 


Fig. 2. 


M=+$ 


(a) (. 


Hyperfine structure and allowed transitions in strong magnetic fields for 
S=}, I=3/2. (a) Magnetic interaction only, (6) electric quadrupole 
interaction added. 


It is convenient to consider these two types of second order terms 
independently. From those which vanish in strong fields, one can separate 
terms in zero, first and second powers of m. The first of these produces 
& constant shift of all lines to higher frequencies (lower magnetic fields) 
which is unimportant except in accurate determinations of the g-value. 
The second produces a change in the separation of the hyperfine lines, 
which is different for different electronic transitions. This is important, 
for it can be used (except when S=4) to determine the relative signs of 
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the parameters D, A and B, which cannot be found from any of the first 
order terms. This point is considered more fully by Bleaney and Ingram 
(1951b), with reference to the manganese spectrum. The third terms, 
in m’, produce a linear change in the separation of the hyperfine lines 
within a given electronic transition. For 6=0 or 90°, the separation 
increases always towards higher fields, but in intermediate directions it 
may reverse. For example, in cobalt ammonium sulphate (Bleaney and 
Ingram 1951a), where there is great anisotropy, the reversal takes place 
only 13° from the perpendicular direction. Note that none of the terms 
is really independent of 0, owing to the variation in g and K with angle, 
unless both g and the hyperfine structure are isotropic. 

The second order terms in Q? are of importance because they form one 
method of detecting the presence of a small quadrupole interaction. 
They contain terms in m, which produce a constant change in the spacing 
of a group of hyperfine lines, and terms in m? which cause the spacing to 
be greater at the ends of a group than in the middle, when 6=90°. This 
effect is reversed as 6 decreases, owing to the opposite sign of the two terms 
in Q?, and both terms become zero at 0=0°. The effects are greatest for 
the transition M=}~—4. If A>B, and g,>g,, the unequal spacing is 
greatest at some intermediate value of 0 between 0° and 90°, but if the 
inequalities are the other way, it may be greatest at 0=90°. The 
quadrupole terms may be differentiated from the other second order 
shifts in two ways: (a) they do not vanish in strong fields; (b) the 
separation of successive hyperfine lines is greatest in the middle or at the 
ends, instead of showing a progressive increase or decrease. 

A second method of discovering a small quadrupole interaction depends 
on the weak lines corresponding to changes of +1 and +2 in the nuclear 
quantum number m, which have intensities of the order of (Q/K)? compared 
with the main hyperfine lines. The positions of these lines involve the 
first power of Q instead of (Q?/K), and we shall neglect small shifts 
of the order Q?/K, and (A, B)2/(g8H,). In the intensities we shall 
neglect contributions of the order (A, B)?/(g8H 9)? which vanish in 
strong fields, unlike the quadrupole effects. We have then lines 
(M, k+3)>(M—1, k+3), 


Energy =Kk+ pees BEN! ae ae eS) 
Intensity =4k?{(I +3) —k*\X, REaee ee peer reli ok. (9) 
where k takes on the values (I—4), (I—#), —(I—3), and for the lines 


(M, m+1)>(M—1, m=1) 
Energy =Km+{K(2M—1)+2Q’m—2y’), Nao er) 
Intensity= {(I+1)?—m?} {I?—m?)Y, pag ci ceeen (11) 


where m takes on the values (I—1), (I—2), .. .. —(I—]). 
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In these equations 
Z 


X= aren {aRMarT) (Ke 


Z Q (ges) sine, rire wetin 2 eb) 


Sr SKM 
2 

Q'=0 (Fa cos? 0-1), Dg ae Goan seen eae me mC 

y'=yByH(Ag, cos? 6+Bg, sin?6)/Kg. . . - . « « ~ (186) 


Fig. 3. 


eit. tu 
is | T | | \ | 


oe | | 
a 0) ee tee ek 


Quadrupole effects in spectrum for S=3, I=3/2, with no anisotropy and 
A=B=5Q. 


The significance of these formule can best be appreciated by reference 
to a simple case. Fig. 3 shows the spectrum at various angles for a 
hypothetical ion where S=}, I=3/2, with no anisotropy (A=B and 
9, =9,), and Q=0-2B. The term yfyH is neglected. As 6 increases from 
zero, four lines corresponding to 4m= -+1 increase in intensity and appear 
at symmetrical points between the two outer pairs of main lines. They 
collapse to two lines at cos @=14/(3), then separate out again and 
disappear. The lines 4m=-+2 are strongest at @=—90°, and are 
symmetrically disposed about the two inner main lines, coinciding with 
them at cos@=1/,/(3). It will be seen that the spectrum is always 
symmetrical, and no information can be obtained about the relative 
signs of Q, A and B. 


Q sp (Fant) cos?@sin?6, . . . (J2a) 
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At fields of the order of 10,000 gauss, the direct effect of the external 
field on the nucleus is no longer negligible, and important information 
about the signs of the constants A, B, and Q can be obtained. The 
term y8yH acts with the quadrupole interaction in increasing the 
separation of satellite doublets on one side of the spectrum, and against 
it to decrease the separation on the other side. If A, B, Q and y8yH 
are all positive, the separation is smaller at the high frequency (low field) 
end of the spectrum when cos @<1/,/(3), and greater for values of cos @ 
greater than this. (The separation is no longer zero exactly at 
cos 6=1/4/(3) but goes through zero at a smaller angle for the satellites 
on the high frequency side, and at a greater angle for those on the other 
side.) Ifthe sign of Q or yByH is reversed, the asymmetry is also reversed. 

When anisotropy is present the nature of the spectrum is not altered in 
principle, but the angular dependence and relative intensities are changed. 
The positions of the satellite lines are the same as in fig. 3, if we plot them 
against the angle % between the nucleus and the symmetry axis, for 
Q’ is simply Q(3 cos? 4—1). 

The term y’ in y8yH now contains a multiplying factor which is just 
the cosine of the angle (4—6) which the nuclear axis makes with the 
external field. If A and B are of opposite sign this factor passes through 
zero, since at some angle the nuclear axis will be perpendicular to the 
external field, as can be seen from fig. 1. 

The considerable intensity of the satellite lines 4m—=--1 at intermediate 
angles makes these generally the most sensitive method of detecting a 
small quadrupole interaction, if there is adequate intensity in the observed 
spectrum. If not, accurate measurement of the line separation, as 
described previously, is the best, but no information concerning the signs 
can be obtained. 

The question of what nuclear information can be obtained from the 
relative signs needs careful analysis. The experiments give the relative 
signs of (Ay), (By) and Q; but the signs of A, B themselves depend on 
that of y, since they involve the product of the nuclear magnetic moment 
and the field of the electrons. Essentially, therefore, the sign of the latter 
is compared with Q, which is itself the product of ( 0?V/dz*) and the nuclear 
quadrupole moment. The sign of the nuclear quadrupole moment can 
therefore be established if crystalline field theory gives the signs of the 
electronic magnetic field and. electric field gradient. The sign of the 
nuclear magnetic moment cannot be found ; this difference arises because 
use is made of an external magnetic field. 


§ 3. Enercy LEevELS AND SUSCEPTIBILITY AT ZERO MAGNETIC FIELD. 


3A. Levels for S=4. 

In zero magnetic field, evaluation of the energy levels is somewhat | 
more complicated than in the strong field case because in general we 
cannot apply perturbation theory. Thus in the case S=} the two 
electronic levels M=+4 and M=~—4 are degenerate in the absence of 
hyperfine structure, and the perturbation denominators arising from the 
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term B(S,I,+8,1,) in the Hamiltonian would be zero. The energy 
matrix splits up into a number of doublets in this case, however, so that 
we can give an explicit formula (14) for the levels by solving the quadratic 
secular equations 


W=— 7 +Q)P+4-HU+1)} 
+$/[{k(A—2Q)+ G}?+ B*{(I+ 3)?—#?j], - oy 9 4d) 
where k=M-+m takes the values (I+4), (I—4), ..., —(I+4); with the 


restriction that only the positive sign is to be taken before the square root 
when k—-+(I+4), and only the negative sign when k=—(I+3). Here 
G=g,BH, and the effect of a magnetic field along the symmetry axis has 
been included, since to do so causes no great complication. The 
‘significance of k is that, except for |k|=I+ 4, the levels are linear 
combinations of the states (M, m) and (M—1, m+ 1). In zero magnetic 
field they are all doublets, the states ‘ and —k being degenerate, except 
for the case k=0, which gives two singlets of separation B(I+ +4). The 
levels given by k= +(I-+4) form a doublet in zero field also, with nuclear 
magnetic quantum numbers m=—I and —I respectively. 

The allowed transitions in zero magnetic field are rather complicated 
when B40, owing to the admixture of wavefunctions in the different 
levels. A simple case arises when a radio-frequency field is applied 
parallel to the symmetry axis, which is the axis about which electron 
and nuclear spin precess. The levels which are found by taking the 
alternative signs outside the square root in (14) each consist of admixtures 
of the same states and transitions are allowed between them except in 
the special cases when either both A and Q are zero, or B=0. These 
transitions, which are all doublets, occur at 

hv=/[k?(A—2Q)?+ B2{(I+4)?—k7}], © 2... (18) 
This formula shows that an accurate value of the quadrupole parameter Q 
may be obtained from the spectrum in zero field in cases where no effect 
can be observed in the strong field spectrum. 

When an r.f. field is applied normal to the symmetry axis the allowed 
transitions are rather complicated. We shall content ourselves with 
pointing out that they can arise between any of the levels found by 
assigning values to k differing by one unit in (14). Here again the 
frequencies of the transitions may be used to find Q: in particular, a 
number of transitions fall into doublets where the separation of the 


components is 2Q(2k—1). This splitting gives Q directly, independent. 


of A, B. 
3.2. Levels for S>}. 


No simple formule can be given for this case except when the electronic 
splitting parameter D is large compared with the hyperfine structure 
parameters. In the first approximation we have, then, for levels where 
M++}, 

W=D{M?—-48(S+1)}+AMm+Q{m2—AT(I+1)}. 2. (16) 


ene 
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These are all doublets in this approximation (except when M or m=0), 
the states (M, m) and (—M, —m) are degenerate. 

For M= +4, the formule given in the last section still apply, in the 
first approximation, if the parameter B is replaced by B(S+4). All the 
levels will have shifts of the order B2/D, due to the off-diagonal elements. 

The allowed transitions in zero field may be outlined as follows. If the 
r.f. magnetic field is parallel to the symmetry axis, only transitions 
internal to the levels M= +3 are allowed, and in the first approximation, 
_ these are given by our modified equation (15). When the rf. field is 
perpendicular to the symmetry axis, the allowed transitions are of the 
same type as in strong fields, viz). AM=+1, 4dm=0, except where the 
levels M= +4 are involved. Here again the remarks of §3.1 apply for 
transitions internal to these levels. Transitions are also allowed between 
the level (3/2, m) and either of the levels given by k=(m-+-4), and between 
(—3/2, m) and either given by k=(m—4), in (14); these alternative 
levels are those corresponding to taking the alternative signs before the 
square root. 

These transition rules only apply when DSA, B. When they are of 
the same order the spectrum will be much more complicated. 


3.3. Susceptibility at Zero Magnetic Field. 


At high temperatures (kT >A, B, Q, D) the presence of a hyperfine 
structure produces almost no effect on the susceptibility of a paramagnetic 
salt. Bleaney (1950) has pointed out that if the susceptibility is expressed 
as a power series in 1/T, the hyperfine structure leaves the term in 1/T 
unaltered, and introduces no term in 1/T?, even for a single crystal. 
The presence of an electronic splitting gives an anisotropic term in 1/T? 
(equation (4) of Bleaney 1950) which averages to zero for a powder. 

At temperatures where kT is of the same order as the splitting 
parameters A, B, Q, D, the susceptibility varies in a complicated manner 
and no general formula can be given. One or two special cases are worthy 
of mention, however, being limiting cases where simple formule can be 
given. These illustrate the general behaviour of the susceptibility. 

The cases of S=4, B=Q=0 has been evaluated previously (Bleaney 
1951 b) since it approximates closely to several salts of cobalt. Reference 
to equation (14) shows that when B=0 the energy levels show a first 
order Zeeman effect when a small field is applied along the axis. The 
magnitude of the Zeeman splitting is the same for all the levels, and the 
susceptibility therefore follows Curie’s law, irrespective of the changing 
population of the hyperfine doublets. If a small field is applied 
perpendicular to the axis, only a second order Zeeman effect is produced, 
and the susceptibility is given by the formula 


Lane eA 
Ng2B? EN yA 17) 
C= area US Te 
i x cosh iS 
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where m takes the positive values I, (I—1),...4or0. (In the latter case 
the terms for which m=0 must be weighted by a factor } in the 
summation.) At very low temperatures y, passes through a maximum and 
then falls to a limiting constant value as T>0. 

For substances where S>4, the susceptibility will not differ, in the first 
approximation, from the case of no nuclear interaction, except for the 
levels M=-44. If these lie at the bottom (D positive) the remarks of the 
previous paragraph apply at temperatures where only these levels are 
occupied, and formula (17) holds if D>A. If D is negative, these levels 
lie uppermost, and they will not be populated at temperatures where the 
deviations due to nuclear interaction would show up. In this case we 
have simply when kT <D(2S—1) (i. e. only lowest levels occupied), 

NS*gii" 
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Note that x, is not really negative because this formula applies only when 
D is negative. 

Paramagnetic ions which do not have Kramers degeneracy in their 
electronic levels (7. e. where S is not half-integral), are not often used in 
adiabatic demagnetization work, and the levels in zero field are not of 
great interest. The remarks of this section apply equally to them, 
however, except when D is positive. The lowest level is then M=0, 
which has no temperature dependent susceptibility contribution in the 
first approximation. The hyperfine splitting of this level is also zero, 
except for second order and quadrupole effects. 


§4, NucLEAR ALIGNMENT. 


Three methods of producing nuclear alignment at low temperatures by 
means of hyperfine structure have been proposed. The first of these 
(Gorter 1948, Rose 1949) utilizes the hyperfine structure due to magnetic 
interaction in a paramagnetic salt. The method proposed is to demagnetize 
a paramagnetic salt to a small field, around which the electronic spin will 
precess. The nuclei are then aligned also parallel to this field by the much 
more powerful field (10°-108 gauss) of the electrons. The splitting of the 
nuclear levels produced by this latter field is of the order of 6-01 to 0-1° K., 
and at such temperatures an appreciable degree of nuclear polarization 
will be set up. 

The second method, due to Pound (1949), utilizes the electric quadrupole 
interaction to split the nuclear levels. This causes, not a polarization, 
but only an alignment of the nuclei, since the levels +-m remain degenerate. 
This is, however, sufficient for experiments involving radioactive emission, 
which depends only on even multiples of the angle with the axis of 
alignment (Spiers 1948, 1949). The latter is, of course, determined by 
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the crystallographic symmetry, and ideally a substance with one ion 
in unit cell is required. Since this method does not require the use of 
paramagnetic salts, it will not be further considered here. 

The third method, suggested by Bleaney (1951 a), is similar to that of 
Gorter in that the magnetic interaction between the nucleus and the 
electrons of a paramagnetic ion is used, but differs in that no residual 
magnetic field is applied. As in Pound’s method, an alignment, not a 
polarization, of the nuclei is set up. The nuclei precess about an axis 
determined by the crystalline field, and again, a substance with one ion 
in unit cell is required. This method has the advantage over Gorter’s 
that lower temperatures will be reached, since there is no external 
magnetic field. 

At first sight, on the other hand, Gorter’s method seems to be simpler 
in that single crystals are not required, and the method can be used with 
any paramagnetic salt. Closer examination of the effects of the crystalline 
electric field shows that this is not the case for small residual magnetic 
fields and almost the same conditions are required in the two methods. 
The reason for this is that, in general, a small magnetic field applied 
normal to the crystalline symmetry axis produces only a second order 
Zeeman effect, and the spins remain aligned along the crystal axis until 
fields large enough to overcome the crystalline splitting (if present) or 
hyperfine splitting are applied. Thus a powder, or a single crystal with 
many ions in unit cell, such as iron alum, will be unsuitable, and a crystal 
with only a single ion per unit cell is required, as in the other methods. 
The magnetic field must then be applied parallel to the axis of this ion. 
In the following paragraphs some attempt is made to calculate the degree 
of polarization and of alignment in zero field that can be achieved in 
limiting cases. 


4.1. Nuclear Alignment and Polarization with B=Q=0. 


For both methods the ideal case arises when the energy matrix of the 
lowest electronic level (assumed a°doublet) has vanishing elements off the 
diagonal. This can occur in two ways : (a) if the lowest levels are M= +43, 
and there is considerable anisotropy so that A>B ; (6) if the lowest levels 
are M4-}, in a salt with considerable splitting of the electronic levels 
in zero magnetic field (D>A, B). Examples of salts approximating 
to these two extremes are cobalt fluosilicate and manganese fluosilicate, 
both of which have only a single ion in unit cell. Then, if the off-diagonal 
elements are neglected, the lowest energy levels consist simply of a set of 
2T-+1 doublets, equally spaced by Ax|M| (cf. Bleaney 1951] a,b). Each 
doublet contains levels in which the electronic spins are opposed and the 
nuclear magnetic quantum numbers are m and —m respectively. The 
top and bottom levels have m=-+I, and the intermediate levels are in 
order of m changing by one unit (see fig. 4a), At a temperature T, the 
combined population of the two states with a given value of m is 
proportional to cosh (AMm/kT). 
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Fig. 4a has been drawn for the case of S=4, I=7/2, which is 
appropriate to the stable isotope, cobalt 59. The lowest temperature 
which could be reached by adiabatic demagnetization of such a cobalt 
salt is equivalent to kT=0-7A, and is represented by a vertical line. 
The populations of the nuclear levels +7/2 : +5/2 : +3/2 : +1/2 would 
then be 1 : 0:53 : 0-32: 0-21, showing a considerable preponderance in 
the states of high m. The radioactive isotope, cobalt 60, is thought to 
have a spin of approximately 4, and if its magnetic moment is about the 
same as that of the stable isotope, there would be a similar preponderance. 
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Energy levels in zero magnetic field for S=}, I=7/2, for the extremes 
(a) B=Q=0, (b) A=Q=0. 

Cobalt ammonium sulphate approximates to (a), and, for comparison, the nuclear 
levels of a diamagnetic cobalt compound in a field of 100,000 gauss are 
shown to scale in (c). The vertical lines in (a) and (b) represent kT for the 
lowest temperatures which could be reached by demagnetization from 1° K. 


If a magnetic field is applied parallel to the axis. a first order Zeeman 
effect splits each doublet by an amount 29,6] M|H. In each case the two 
levels originating from the same doublet have opposite signs of m; for 
the upper doublets the lower component is that with negative m (assuming 
A positive), but for the lower doublets this is reversed. If it be assumed 
that the degree of nuclear polarization produced is measured by the ratio 
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of the total nuclear magnetic moment o along the axis to the saturation 
nuclear magnetic moment o,,, (cf. Rose 1949), a simple expression for 
this ratio can be found. It is 


; M|A 
+'m sinh | 
(o/o eaten lou, sas , (19) 
Es i »' cosh Tea ee ae 
m kT 
where m takes the values I, [—1, ..., $ or 0. (In the summation, the 


terms for which m=0 should be weighted by a factor }.) 

As an illustration, we will suppose that a field is applied which 
would split each doublet by an amount equal to A. In cobalt this would 
require a field of about 100 gauss, and the temperature would rise so that 
approximately kT=A instead of 0-7A. Then (19) gives (a/o,.,,.)=0-3 
under these conditions. This calculation is, of course, based on the 
assumption that the cooling of the salt is due entirely to the 
paramagnetism of the cobalt ions. If the temperature could be lowered, 
for the same residual field, by the use of another paramagnetic salt in 
contact with the cobalt salt, the degree of nuclear polarization could be 
improved. 


4.2. Nuclear Alignment with A=Q=0. 


It is of some interest to examine the other extreme case which can occur 
for S=4, namely, when both A and Q are negligible compared with B. 
The energy levels in this case are given by the equation 


W=+By([(+3)—#'), 


where k takes all values (I+4), (I—4), ..., (—I+4). Each of these 
levels is a singlet composed of equal admixtures of the nuclear states 
m=k +1, but (except when k=0), the levels +k coincide to give doublets. 
The energy levels for [=7/2 are shown in fig. 4 b, the scale being such that 
B in fig. 46 is assumed equal to A in fig. 4a. The vertical line again 
represents the value of kT at the lowest temperature which could be 
reached by adiabatic demagnetization of such a salt. At this temperature 
the populations of the nuclear states +7/2 : +5/2 : 13/2: +1/2 would be 
6-34 : 0-62: 0-86:1-0. The predominance lies with the states +3, 
corresponding to the fact that in this case we have a “ plane magnet ”’, 
i.e. @ magnet precessing in the plane normal to the symmetry axis. 
The nuclear alignment is a good deal less favourable than in the preceding 
case of B=0 because of the concentration of energy levels near the 
bottom. 

As already pointed out, a small magnetic field applied in any direction 
produces only a second order Zeeman splitting, so that this is not a 
favourable case for Gorter’s method. 

~~ From these extreme cases of B=Q=0 and A=Q=0 it is easy to see 
what happens in intermediate cases, which are obviously less favourable. 
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If A is large but B is not entirely negligible and | M |=}, the energy levels 
in fig. 4a are displaced by amounts of the order B*/A, and each level m 
has admixed with it the state (m—1) or (m+1) to an extent (B/A)?. 
This is not serious until B becomes comparable with A. The case B=A 
corresponds to no anisotropy, and the energy levels fall into the usual 
two groups for S=4, corresponding to a new quantum number F with 
values I+4 and I—4. The splitting of this doublet (for I=7/2) is 4A=4B 
and each level contains an equal admixture of all the nuclear states m, 
so that no nuclear alignment in zero field is achieved at any temperature. 
This corresponds to the fact that with no anisotropy there is no preferred 
axis of alignment. As the value of B increases further relative to A the 
situation gradually changes to the extreme case A=0 which has already 
been considered. 

The remarks of this section can be applied when S>4, if the levels 
M= +4 are lowest and D> A, B, if the parameter B is everywhere replaced 
by B(S+4). Thus a substance of this kind will approximate to the case 
of A negligible, even if A=B, when S is large, since the essential criterion 


is A<B(S+4). 


§ 5. Discussion. 


From the point of view of entropy, the use of the same paramagnetic 
salt for cooling and alignment can be summarized as follows. When the 


salt is magnetized in a very strong field at 1° K., the electronic entropy 


is completely removed. On demagnetizing to zero field, the entropy 
remains constant, but is divided between the electron and nuclear spins. 
In the simple cases considered in §4, each level is a doublet, and in the 
two levels of every doublet the electron spins are oppositely oriented. 
Thus at all temperatures the “ electronic entropy ”’ will be R log 2 in zero 
field, with a corresponding decrease in the nuclear entropy, giving a degree 
of alignment. When a magnetic field is applied, the electronic entropy 
decreases, and the nuclear entropy increases, until at fields large enough 
to separate the electronic levels so that their hyperfine patterns do not 
overlap, the electronic entropy is again zero, and the nuclear entropy 
R log, (2I-+-1), with no polarization or alignment. 

For this reason the discussion of the merits of the Gorter-Rose method 
has primarily been limited to the case of small magnetic fields. The 
formule given by Rose correspond, however, to the case of an applied 
field so large that the hyperfine structures of the different electronic levels 
no longer overlap (cf. fig. 2). Rose assumes that only the lower electronic 
level is occupied, and further, that the temperature is so low that the 
nuclear levels corresponding to the splitting of this electronic level are 
not all equally populated. This can only be achieved by the use of a second 
paramagnetic salt to produce the cooling, but it would be necessary to 
adopt an experimental arrangement whereby the external magnetic 
field used to produce the nuclear polarization is not applied to the second 
paramagnetic salt. 


~~" 
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The requirement of a “dilute ” paramagnetic salt for these experiments 
has already been pointed out (Bleaney 1951a). The reason for this is 
that the random magnetic fields of neighbouring ions broaden the energy 
levels, and if these overlap considerably, the alignment. will be largely 
destroyed. The criterion for the degree of dilution required is that the 
hyperfine structure shall be reasonably well resolved in the paramagnetic 
resonance spectrum. Thus salts of the iron group will require “dilution” 
with isomorphous diamagnetic compounds, but in the rare earth group 
this will probably not be necessary if the ethyl sulphates are used. These 
have only one ion in unit cell, and hyperfine structure lines can be 
resolved in the neodymium salt without dilution (Bleaney and Ingram 
1949). The use of such dilute salts is also advantageous in cooling, since 
lower temperatures can be obtained; a disadvantage is that the rate of 
warming is larger owing to the smaller heat capacity. 


§ 6. CONCLUSION. 


In the first section of this paper general formulz are given for the 
allowed transitions in the paramagnetic resonance spectrum. These are 
used to determine the magnitude of the parameters in the Hamiltonian (1). 
When a hyperfine structure is present, a certain amount of nuclear data 
may be obtained from the spectrum, the simplest datum being the value 
of the nuclear spin I from the multiplicity. From the experiments, 
values of the hyperfine structure constants A and B, and Q if itlis 
comparable in magnitude with A and B, can be found with fair accuracy. 
The relative signs of the parameters A, B, Q and D are much more 
difficult to elucidate from the spectrum, since they can only be found 
from second order effects. Where a crystalline splitting of the electronic 
levels is present (D0), second order effects in the overall size of the 
hyperfine structure in intermediate fields can be used for this purpose. 
An example of this is manganese (Bleaney and Ingram 1951b), where, 
in addition, a subsidiary experiment on the anisotropy of the 
susceptibility gave the sign of D, and so determined all the signs. 

For a salt where S=4, relative signs can only be found when weak 
transitions involving changes in the nuclear magnetic quantum number 
are observable. In strong fields such transitions arise from the electric 
quadrupole interaction, but, of the salts so far investigated, only those of 
copper have shown such an effect (Ingram 1949). In other cases the 
relative signs of A and B can in principle be found from the second order 
transitions 4Jm=-+1, which occur in weak fields, using the small 
displacement of these transitions due to the direct effect of the external 
field on the nuclear moment. This method is not generally feasible, 
since this effect is too small to detect except in strong fields, and there 
the transitions vanish in intensity, if the quadrupole interaction is 
negligible. 

In cases where no weak transitions 4m= +1 are visible, or no anomalous 
spacing of the hyperfine structure is detectable, the formule given in §2.2 
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can be used to obtain upper limits for the magnitude of Q. A more 
sensitive method is observation of the allowed transitions in zero magnetic 
field, as emphasized in §3. This necessitates covering a wide wavelength 
range, generally in the decimetre wavelength region, where the sensitivity 
of the method, using a small single crystal to obtain the greatest accuracy, 
is rather low. Such experiments are under way in this laboratory. 

In this paper the author has not attempted to consider the relation of 
the parameters D, A, B, and Q to the electronic state of the ion and the 
crystalline field. The intricate theory required for this purpose has been 
the subject of a general review by Abragam and Pryce (1951a). It seems 
that the theory is now sufficiently accurate to enable quantitative data 
concerning the nuclear moments to be deduced from the paramagnetic 
spectra in many cases (see, for example, Elliott and Stevens 1951). 
Relative values of the moments of different isotopes of the same element 
can be obtained, of course, directly from the spectra. 
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SUMMARY. 


This paper contains an approximate method of calculating the communal 
entropy of an assembly of monatomic particles. The method of dealing 
with dense systems may be regarded as an extension of that used by 
Lennard-Jones and Devonshire in their theory of liquids and dense gases. 
The available volume is divided into cells and the communal free energy 
expressed in terms of a set of parameters related to the probability of two, 
three or more particles occupying a given cell simultaneously. Application 
to an assembly of rigid spheres leads to the conclusion that the communal 
entropy does not become appreciable until the available volume is over 
five times that of a close-packed assembly. This suggests that, for more 
accurate intermolecular potentials, the communal entropy is practically 
zero in the solid and liquid states, but that the extra terms in the free 
energy may have an appreciable effect on the calculated critical constants 
and vapour pressures. 


§1. INTRODUCTION. 


AccorDING to classical statistical mechanics, the thermodynamic functions 
of a system of N particles confined to a volume v can be derived from the 
free energy A given by 


ImmkT\2N2 pee 
exp (—A/kT)= (AS [--- [exp (—V/AT) do, ... dos, (1.01) 


where V is the potential energy of any configuration and. fr dv, indicates 


integration of the coordinates of particle 7 over the whole volume v. 

The only exact method of calculating the integral in (1.01) is that of 
Mayer (1937). However, this only applies to a limited range of densities 
and is exceedingly difficult to handle mathematically. Ofthe approximate 
methods used many (Eyring and Hirschfelder 1937, Lennard-Jones and 
Devonshire 1937) have been based on a model in which each particle is 
confined to a cell by its neighbours. Lennard-Jones and Devonshire 
calculate the field in such a cell on the assumption that the neighbours 
are in their equilibrium positions. An improved form for the cell field has 
recently been suggested by Kirkwood (1950). 

Although these cell models may give a reasonable picture of molecular 
environment at high densities, they will become increasingly inaccurate as. 
the cell size increases. At low densities restriction to cells will prevent 
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collisions occurring and the method gives no second virial coefficient. An 
attempt to remedy this situation and to extend the cell method to low 
densities has recently been made by Buehler, Wentorf, Hirschfelder and 
Curtiss (1951). In addition to the model used by Lennard-Jones and 
Devonshire, which they call the “empty centre model’’, they consider 
another in which the wandering particle moves in the field of particles 
fixed at the centres of neighbouring cells and one fixed at the centre of the 
cell in question. This is called the “occupied centre model’’*. The 
arithmetic mean of the two free volumes is then used to construct the 
partition function. This method reduces to that of Lennard-Jones and 
Devonshire at high densities and gives the correct second virial coefficient 
at low densities. No theoretical justification is given for this averaging 
process, however; in the intermediate region the arithmetic mean formula 
is no more than a convenient method of interpolation. 

When the cell method was originally developed, it was pointed out by 
Hirschfelder, Stevenson and Eyring (1937) and also by Lennard-Jones and 
Devonshire (1937) that the approximation of placing one particle in each 
cell, from which all others were excluded, constitutes an artificial restriction 
on the region of phase space over which integration is to be performed. It 
can be shown that, in the limit of low density, this corresponds to the 
omission of a factor e§ in the partition function and leads to an error of k 
per particle in the calculated entropy. This extra term is usually referred 
to as the communal entropy ; it has been discussed in some detail by 
Lennard-Jones (1940). The chief purpose of this paper is to examine this 
term and show how it is connected with the work of Buehler et al. 

Although we know that the communal entropy increases from zero to k 
per particle as the volume increases to infinity, no satisfactory theory of 
the point at which it becomes effective has been advanced. Originally, it 
was suggested by Hirschfelder, Stevenson and Eyring (1937) that it 
becomes available at the melting point. This hypothesis has been 
examined by O. K. Rice (1938) who concludes that it is without sufficient 
justification and that the communal entropy is unlikely to play a major 
part in the fusion process. We shall see that the theory of this paper 
confirms Rice’s conclusions. 

In this and subsequent work (1944) Rice has discussed the communal 
entropy of an aggregate of hard spheres in some detail. In the later paper 
he claims to show that for a gas of hard elastic spheres the communal 
entropy is fully excited for each direction of space, and amounts in all to 
3. R per mole. As these results are at variance with those obtained in 
this paper, it is well to find the cause of the discrepancy before proceeding 
further. The difference lies in the meaning attached to the phrase 
“restriction to a cell”. Rice interprets this as meaning that the whole of 
the molecule must lie within the cell. Although this has a precise meaning 
in the case of hard elastic spheres, it leads to difficulties when we deal with 


ee ee ee eee ee 
* In the published version of their work, Buehler et al. have used the terms 
“‘soft-centre model ” and “‘ hard-centre model ”’. 
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actual molecules which have no fixed boundaries. Here we shall say that a 
given particle occupies a certain cell if its centre (or some fixed point of the 
molecule if it is non-central) lies in that cell. 

In this paper we shall discuss the way we should set about calculating 
the communal entropy and free energy in a general manner. These are 
expressed in terms of a set of quantities w,, wz, .... which are related to the 
probabilitity of two, three or more particles occupying a given cell 
simultaneously. It is shown that w, is the most important of these and 
that a good approximation to the extra terms in the thermodynamic 
functions can be obtained if we neglect the others. This is followed by 
a simple model for the calculation of w,. With some further 
approximations, it is found that w, can be expressed in terms of the free 
volumes of the empty and occupied centre models used by Buehler et al. 
In this way we have obtained a formula for the communal free energy 
which is, in effect, a theoretical foundation for the use of these models. 
The most important result of the investigation is the modification of the 
arithmetic mean formula. If v,, and v,, are the free volumes of the 
empty and occupied centre models respectively, it is suggested that the 
true free volume should be calculated from 


Up O FeV) (2 psO pe) 1 ee Me) Caiky sy (1.02) 


At the end of the paper, this formula is applied to estimate the communal 
free energy of a system of rigid spheres, using the values of v,, and v,, 
calculated by Buehler et al. 


§ 2. A GENERAL METHOD FOR CALCULATING THE COMMUNAL ENTROPY. 


We shall base our work on the recent analysis of the cell method given 
by Kirkwood (1950), using a similar notation. If we write 


La it ‘ . | exp (Vil) due... dtx, . «= = (2.01) 


the free energy is given by 


A 3 QamkT 1 1 
<r = — 5108 —as— a a ie (2.02 
NET 5 log RB nos Z+ N log N (2.02) 
-Kirkwood now supposes the volume v to be divided up into a lattice of N 
equal cells and writes Z("1'-”s) for the restricted phase integral in which 


the cells are occupied by m, . . . my particles respectively. Z/N! can then 
be divided up as follows : 
eg eer en) gc, 2. (208) 
N! M1...My IIm,! 
=m;=N 8 


The usual cell method is concerned with an approximate evaluation of 
Zo---D. Tf we write 
Z 


Se es se ess (204) 
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the communal free energy may be defined as —NA&T log o and the 


§ oT . 


From (2.03) and (2.04) we see that the value of o will depend on the ratio 
of the phase integral with multiple occupation Z(™--""x) to Z@--). 
We shall suppose that this ratio can be expressed as a product of the form 


Z,(m +++ Mx) 


N 
7ilaet) ae aes N20 cal ee eat OY 


where there is an w-factor for each type of occupation (vacancy, single 
occupation, double occupation, etc.). Clearly we must have w,=1. Also 
we may take w)=1 without loss of generality, for if w,~1 we may replace 
it by 1 and at the same time replace w,(i 40) by w), | w, without altering 
IT Om, 

The approximation (2.05) is equivalent to assuming that multiply 
occupied cells do not interfere with one another. It will be accurate at 
high densities when there are very few multiply-occupied cells and at low 
densities when the cells are large and effectively independent. It therefore 
seems reasonable to use it throughout the full range of volume. w. may be 
described as the ratio in which the restricted phase integral is altered if we 
remove a particle from a singly occupied cell and move it to another also 
singly occupied, leaving a doubly occupied cell anda hole. ws, a4...... 
may be defined in similar ways. 

We next have to solve the problem of finding oc when we are given the 
values of the w’s. The full expression, 


N N 
poke ESy i (2), » \Doxuy aoa GC 


My++-My S=1 m,! 
2m;=N 


obtained from (2.03), (2.04) and (2.05) has to be converted into a less 
unwieldy form. 
Suppose that a given partition m,...my contains n;(m,... mx) 


cells with 7 and not more than 7 particles. Then, since we are dealing with 
N particles and N cells, we must have 


N 
+ in=N, 
i=0 
N (2.07) 
n, =N. 
i=0 
The number of partitions m, ... my with NM empty cells, n, singly occupied 
cells, etc., is 
N! (N—n))! N! 
tg! (N—no)! ny(N—m,—n,)!'°"  mg!nil...myl’ * (2.08) 
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Equation (2.06) can therefore be written 


N! aN 
of eae RUT ee eee 09) 
Noy... Ny O+ Ny: 2+. NN: s=0 8. 
where 2” indicates summation over those values of n,n... . 2 which 


satisfy (2.07). As we shall only be concerned with the logarithm of o we 
may follow the usual procedure and replace the series in (2.09) by its 
largest term. We therefore seek the unconditional maximum of 


N! N /w,\"8 A\ XN N 
OOP eNactgisartal LAUGH) at + (Pen) Zmctlowe) 3 ino 


* s=0 
yp ee ary) 


log (A/N) and log pw are Lagrange multipliers taken in this form for 
convenience in the subsequent analysis. 
Using Stirling’s formula we find that the maximum term occurs when 


n;=NadAp' a, /0! Koa Re reN) i alt ay meek V (oa) 


where the parameters A and p are determined from the conditions 


N Gin: ) 
aK (2.12) 
yee stan 
eee EN) 
If we define 
N 
A eee 1) eee | noe. 5: (2.13) 
i=0 
(2.12) may be written 
A t==f (ppl (pe) rh ptt seer (2.14) 


Finally we have 


N! NON? 1/N l 
== eae ase — — / ; : a 
a ie 22 Ny! a (=) pes Me f’(e) (2.15) 


The complete solution may therefore be written in implicit form 
o= f(wpaf'(u) . - . «+ + (2.16) 


§ 3. THE CoNTRIBUTION OF DoUBLE OCCUPATION, 


If all the w’s except w, vanish, equation (2.16) can be solved explicitly. 
In this case 


f(x)=1+2+40,27, ee tee ee ol) 
p=v/(2/w), 6=1-+-4/(20.). - . ~ « (3.02) 


This leads to a free energy 


A QarmkT 
N&T he 


whence 


= — Bog I og Z'-P—log {1+ Vi2eg)}. (8.03) 
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This is the solution we get if we replace the limitation of single 
occupation by a condition which prevents more than two particles 
occupying a cell simultaneously. We know that w, varies for 0 to 1 as the 
volume increases, so (3.03) shows at once that at low densities log {1+ 1/2} 
or 88 per cent of the communal entropy is due to double occupation of 
cells. It is also possible to get the equation of state from (3.03) by the 
usual formula 


DONS 6 SAIN fo erates 
BT duis, NEL)> 


To find w, at low densities we may neglect the effect of particles in other 
cells and write 


N2 4 ;A | 
te ee | exp [— E(r,9)/kT] dv, dv, 
N2 ;4 cA 
= ite om | dv | {(1—exp [—E(r,)/KT]} dvg, . . (3.05) 


where E(r,,) is the intermolecular energy of the pair of particles 1 and 2. 
As far as terms of order v-1 are concerned, we can replace integration over 
A with respect to v, by integration to infinity, so that we get 


oy=1— 2S | * (1 exp [—E(n)kTI} dr, 2 ay ee GE 


From (3.03) and (3.06) we now get the equation of state 


Ee 1 

foes ite se | “(1 —exp[—B(r)/kT]}dr-+-0 (=) _ (3.07) 
From this equation we see that the theory allowing for double occupation 
gives a second virial coefficient equal to \/(2)/[/(2)+1] or 58-6 per cent of 
the correct value. Thus, even in the limiting case, (3.03) corrects the major 
part of the deficiency of the original cell method. At higher densities, the 
value of ws, w4, etc. will be smaller relative to w, and (3.03) will become even 
more accurate. It seems very probable that, in the critical region, the 
terms ws, w,... will, in fact, be quite unimportant, so we shall not discuss 
them further. 


§4. MopELs FoR EVALUATING wp. 


Tbe next step in the theory is to set up simple models for the calculation 
of the multiple occupation factor w,. To find this we shall use a modifica- 
tion of the method of Lennard-Jones and Devonshire, in which each 
particle is treated as moving independently in a cell field ys (r). We shall 
make the further approximation of supposing that this field is unaltered 
when a neighbouring cell becomes vacant or doubly occupied. The errors 
in the two cases to some extent cancel. ws can then be expressed in 


terms of the phase intergal for two particles moving in the ordinary cell 
field y (r). If we write 


vi (r)=(r)—P(O)p- 2 Le. A. (4,01) 
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this expression for w, is, from the definition in § 2. 
i iF exp (—[¥r(t1) +¥1(t2)+E(r 2) /kT} dv, dog, 
Wo= a a ee i eee En aan BOL Dae 
le exp [—y,(r)/kT] dv \* 


Aves : 
where | indicates integration over the cell. Similar formule could be 


(4.02) 


given for ws, w4, etc. 
If the integral in the numerator of (4.02) is written as I,, we have 


T= [° Q(ry) exp [—(r,kT] dv, . (4,08) 


where Q(r)= iP exp {—[y(t2) +E (7,9) /kT} dvg, =e et4.04) 


Q(r,) may be described as the free volume of particle 2 in the cell, given that 
particle 1 is the point r,. 

In the case of a spherically symmetric cell, Q(r) can be expanded about 
the origin as follows : 


pri /Q2iQ 
Q(r,)= Zan (Sn), ee ( 4.05) 
The integral I, then becomes 
co a7'Q 1 ay 
i= a (S=) om | rexp[—x(r)/kT] dv. . . (4.06) 
0'Q 
The coefficients ( a *), can be written 
g2i A ( g2i 
(S),-/ a 5 exp [—Hiri. yer) exp [—w(7°9)/kT] dvy, (4.07) 
leading to 
“A 
Qo= | exp {alr HE()V/KT} de, 
oe me a), 4 [52 E 
7 (52), -- wl Eas. (z” — zr) | exp {—[hy(r) LE(r)/kT} do. 
; (4.08) 


A simple first approximation would be to neglect all but the first term in 
(4.05) and write Q=(Q)). Examination of equation (4.08) shows that 
this is just the free volume of the occupied-centre model used by Buehler, 
Wentorf, Hirschfelder and Curtiss (1951). In their notation v,, is the free 
volume in an empty cell and v,, is the free volume when there is a particle 
at the centre*. Thus to this approximation 


Q=%; 1, =e Y fo We=Vfo/ fe . ° ° ° (4.09) 
If we neglect w; for 1>2 the free energy is given by 
A 3, 2amkT @ 
NET 5108 a —log vs.— NUT —log {1+ +/(2v%/v;-)}, - - (4.10) 


* Buehler et al. use vy and vp, in the published version of their paper. 
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where (—®,) is the potential energy of the system when each particle is at 


the centre of its cell. 
This formula could be used to calculate the equation of state from the 


thermodynamic relation 
0A 
pa —, (a ') os ile, fey, tLe ate 
p ( ms (4.11) 


In their modification of the original work of Lennard-Jones and Devonshire, 
Buehler, Wentorf, Hirschfelder and Curtiss use the formula 


A 2amk'T 


®D 
—log Us— NUT 


where 
OF =2(Ue+Yfo) MEE ar 


This formula, however, is only an empirical way of interpolating the 
equation of state, and does not give the correct entropy changes. From 
the theory of this paper, it appears that a more satisfactory procedure 
would be to replace the arithmetic mean formula (4.13) by 


Us=VUpet V/ (2UF¢ Vso) . . . . . . . (4.14) 


Buehler, Wentorf, Hirschfelder and Curtiss (1951) have found vw, and 
UV for a gas of rigid spheres. Using their figures we can easily get w, 
and the communal part of the free energy as functions of the volume. The 
results are shown below. v, is the least volume into which the spheres 
of diameter 7, can be packed. 


_ TABLE. 


The communal free energy of a gas of rigid spheres. 


V/V Vpe/4875 Vfo/ 4875 We —(A/NkT) communal 
1-000 — — — — 
2-828 -0102 — - — 
4-887 -0550 -OO18 ‘032 “044 
7-761 1143 ‘0272 -238 -290 
10:0 +1473 ‘0600 ‘407 -454 
15-0 -2210 a BEY ‘605 -618 
20-0 -2946 -2073 -704 ‘691 


Examination of this table shows that the communal entropy does not 
become appreciable until the total volume is about five times the volume 
of the close-packed lattice, It seems probable that the results would be 
similar if other, more accurate, inter-molecular potentials were used. The 
general conclusion is, therefore, that the communal entropy is practically 
zero in the solid and liquid states, but may form an appreciable part of the 
entropy of evaporation. Cell theories of critical phenomena, such as that 
of Lennard-Jones and Devonshire, should be corrected to allow for it. In 
the application of their theory to vapour pressures and boiling points, 
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Lennard-Jones and Devonshire (1938) included the additional factor e% 
in the liquid partition function. At the end of their paper, however, they 
stated that if this factor was omitted, the calculated vapour pressures 
would be reduced by a factor of e and the calculated boiling points by about 
9 per cent. The present work shows that these corrections should be 
applied. 


The author is indebted to Professor Sir J. Lennard-Jones for the loan of 
an. advance copy of a paper by Buehler, Wentorf, Curtiss and Hirschfelder, 
sent to him by Professor Hirschfelder. 
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ABSTRACT. 


By an extension of a recent analysis of the problem of grain boundary 
diffusion in metals, the expression 


aes : . 2 (wDt)!? cot? « 

has been derived relating the diffusion coefficients D’ and D for grain 
boundary and for volume diffusion, 5 the effective “‘ width ” of the grain 
boundary and « the angle at time ¢ between a line of constant concentration 
and the grain boundary at the point where it meets the grain boundary. 
Some remarks are made on the use of the equation and when it is applied 
to a recent experimental observation on grain boundary diffusion of copper 
in nickel, the result is deduced that D’=8 cm.2/day at 1000°C. 

A discussion follows on some aspects of the relation of other work on grain 
boundaries to grain boundary diffusion. Evidence is presented, based on 
models of the grain boundary, which suggests that the activation energies 
for grain boundary and for volume diffusion may not differ appreciably. 


$1. INTRODUCTION. 


Work on preferential grain boundary diffusion in metals has in the past 
been hampered by the lack both of a quantitative treatment of the process 
and of reliable experimental methods for studying it. It is the purpose of 
this note to extend a recent analysis (Fisher 1950) of the problem in such a 
way that the simple qualitative method of study reported from this 
laboratory (Barnes 1950) can be made to yield quantitative results with 
very little extra work. A discussion is then given of the importance of such 
quantitative studies insofar as they may give information on the nature of 
the grain boundaries themselves. 


§2. QUANTITATIVE METHODS FOR GRAIN BOUNDARY DirFrusIon STUDIES. 
(a) Fisher's Analysis. 

Fisher (1950) has recently obtained, with certain assumptions, an 
expression for the concentration of diffusing material in the neighbourhood 
of a grain boundary in terms of the time ¢ during which diffusion has been 


* Communicated by the Author. 


On the Grain Boundary Diffusion in Metals 469 


taking place down the grain boundary, the coordinates 2, y, (see fig. 1) 
and the coefficients D’ and D for grain boundary and for volume diffusion 


respectively. 
Ben (Mv (2) eia\ 
C= Exp ( Ae . Erfe an) oy Rie eked ose (1) 
Dt x y (D' 
where eS oy > t= 5° Oi 5 (=) 4 


The grain boundary commences at the point z=0, y=0 and 8 is its effective 
““ width ”’. 


igs lie 


LINE OF CONSTANT 
CONCENTRATION AS 
REVEALED BY ETCHING. 


DIRECTION 
OF 
DIFFUSION 


GRAIN BOUNDARY (xc =0) 


Diffusion down a grain boundary. 


As a confirmation of the validity of equation (1), Fisher has successfully 
applied it to an interpretation of the results of an experiment on the grain 
boundary diffusion of radioactive silver in polycrystalline silver, by the 
standard technique of measuring the activity of layers of silver 
perpendicular to the diffusion direction, machined from the diffusion zone. 
He obtains with the assumption 5=5x10-§, 7.e. about two atomic 
spacings, 

DYD=9-4% 108. 
‘(b) Author's Extension. a 

More recently a publication from this laboratory (Barnes 1950) has 
described how the contours of concentration in the neighbourhood of a 
grain boundary, down which preferential diffusion (in this case of copper in 
the grain boundaries of nickel) is taking place, can be revealed by etching. 
It was suggested therein that the angle « between the grain boundary and a 


2 % 
* Erfe c=1— a4 e dt, 
0 
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concentration contour at the point where it meets the grain boundary could 
be used as a measure of the ratio of the rates of diffusion down the grain 
boundary and into the grain. This ratio of course is only of secondary 
importance but by development of equation (1) measurements of the angle 
a can be used to calculate the ratio of the actual diffusion coefficients for 
grain boundary and for volume diffusion. 
Simple differentiation of equation (1) gives the slope at any point (x, y) of 

a line of constant concentration c at time t. From this is readily deduced 
the equation 

tes  o(nD1)12 cot? « (2) 

D735 b. oth ote) ae 
connecting the ratio of D’ and D with the angle «. When this equation is 
applied to the previously published work on copper-nickel, taking «—30° 
(see fig. 1, Barnes 1950) D~10-> cm.?/day and assuming again 6=5-10~, 
we find 

D’/D=8-10* at 1000°C., 


a ratio the order of magnitude of which is the same as that obtained for 
self diffusion in silver. 
Inserting the known value of D into this equation gives 


D’'=8 cm.?/day. 


Fisher’s analysis provides the first real attempt to study grain boundary 
diffusion in metals in a quantitative way and it is suggested that the 
extension of his relation to give equation (2), together with the experimental 
method for providing values of «, constitute a simpler and quicker means of 
obtaining grain boundary diffusion coefficients in metals than a sectioning 
technique such as used by Fisher. Possibly many micro-photographs. 
showing diffusion in metals may already exist in which a preferential 
grain boundary penetration is apparent and to which the relation (2) 
could be usefully applied. 

A calculation of D’/D from equation (2) requires some assumption as to. 
the magnitude of 5, the width of the grain boundary. This is usually stated 
as being a “few atomic diameters”, this estimate being based on the 
assumption of the transitional lattice theory of the grain boundary and the 
assertion of very short range forces between atoms. A precise evaluation. 


of 6, and so of D’, is not possible at present, but tor practical purposes we - 


may use D’d in place of D’ to characterize grain boundary diffusion rates, 
for D’é clearly represents the amount of material diffusing in unit time 
under unit concentration gradient per unit length of the grain boundary 
in the “‘z’’ direction (fig. 1) and perpendicular to that length. D’S can be 
obtained by measurement of only «, D and t, and requires no assumptions 
about the structure of the grain boundary. The results for copper in 
nickel give 
D’s/D=0-04 cms, 


nt — ie 


ee 
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There are two points that need bearing in mind in any use of equation 
(2). The first is to ensure that the grain boundary is normal to the plane 
in which « is being measured, and that the grain boundary has not migrated 
during the diffusion anneal. The second point is to remember that the 
equation (1) is derived on the assumption that D and D’ are independent of 
concentration, so that the value of D’/D obtained will only be a mean value 
for the concentration range covered in the experiment. It should be 
noted that any variation of « with the position of the point on the grain 
boundary at which it is measured will indicate a variation not of D’/D but 
of D’/D?'? with concentration. (Such a variation was in fact found in the 
experiment on copper andnickel.) It can probably be assumed that each 
value of D’/D corresponds, although probably only roughly when D’/D is 
strongly dependent on concentration, with the concentration at the point 
where the corresponding value of « was measured. 


§3. Graty BounpaRy MopELs anpD DIFFUSION. 


It is interesting to consider how other work on the properties of grain 
boundaries might be related to grain boundary diffusion studies. 

In view of the difference between D’ and D it is important to 
know whether it is due to a decrease in the activation energy “Q” for grain 
boundary diffusion over volume diffusion, or to an increase in the constant 
term D, in the expression usually used to express the temperature 
dependence of diffusion rates. 


DE Dre CORT eee nt (3) 


We may perhaps draw some conclusions by considering some of the 
grain boundary models that have been proposed. 

To explain his result that the activation energy for grain boundary slip 
in several metals is the same as that for their volume diffusion and for the 
constant rate creep of single crystals, Ké (1949) has suggested the concept 
of “disordered groups ” of atoms as the units responsible for all three 
processes. These exist in certain concentration in the body of the crystal, 
but in much higher concentration at the grain boundary which therefore 
consists of regions of good fit separated by regions of disorder, the * dis- 
ordered groups”. The rearrangement of atoms within a group normally 
constitutes a diffusion process, but in the presence of a shear stress leads to 
slip. The equality of the activation energies for the three processes 
considered by Ké is thus satisfactorily explained. But the same 
rearrangement of the atoms of a disordered group at the grain boundary 
will lead to grain boundary diffusion, so that on Ké’s model, the activation 
energies for volume and grain boundary diffusion will also be the same. 

A study of the temperature dependence of D’ and comparison with that 
ot D would both settle the relative values of the corresponding activation 
energies Q’ and Q and provide useful confirmation or otherwise of such grain 
boundary models as this one of Ké. In such studies one must remember 
that any impurities which segregate in the grain boundary will most likely 
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have an effect on Q’ out of all proportion to their total concentration. As 
an example of the great importance of impurities in grain boundary 
studies, we may quote Ké (1949) who found that even very small quantities 
of impurity may completely block grain boundary slip. 

Dislocation models of the grain boundary have been proposed in the 
past (Burgers 1940, Bragg 1940) and recently studied by Read and Shockley 
(1950), some of whose conclusions have been verified by Aust and Chalmers 
(1950). On this model the grain boundary is made up of a suitable array 
of dislocations. With certain assumptions it can be shown that the 
measured value of D’8;D can be explained on the further assumption that 
Q’ = Q. 

Dislocations exist also, of course, within the crystal, and we will suppose 
that for our present purpose the grain boundary differs from the crystal 
only in that it contains a much higher concentration of dislocations. Seitz 
(1950) has suggested that it is dislocations which provide the very high 
density of vacancies in crystals required to explain certain phenomena, 
particularly the Kirkendall Effect. We will assume then that the main 
sources of vacancies are the dislocations, and further, that the density of 
vacancies in any region is proportional to the density of dislocations. The 
vacancies provided by the dislocations will then lead to diffusion by the 
normal vacancy mechanism, so that we can write 


D= const. p exp (—Q/RI),. 5) eae ee 


where D is the diffusion coefficient in a region where the density of dis- 
locations is p. Q is the activation energy for the diffusion process and does 
not include the energy of formation of a vacancy, as it does when it is 
supposed that the vacancies exist in thermal equilibrium with the lattice. 
In the present case the energy required to form a vacancy, from a disloca- 
tion, is negligible. If we suppose first that Q and Q’ are not equal, we can 
write for D’5/D, from equation (4), 

D'S _ pg, exp(—Q’/RT) i. 

D STEIN VERE 

where p,, and p, are the densities of dislocations per unit length of the 
grain boundary and per unit area of the crystal respectively. We should 
expect the excess density of dislocations at the grain boundaries to provide 
an excess density of vacancies for some distance from the grain boundry, so 
that for diffusion the effective width 5 of the grain boundary will be greater 
than the distance over which the normal lattice structure is disturbed. We 
avoid having to make an estimate of 5 by considering in equation (5) the 
ratio D’d/D instead of D’/D. 

With estimates of p,, and p, we can calculate the difference between Q’ 
and Q that would be necessary to account for the measured ratio of D’8/D. 
We may take for p, the value of 108, an estimate for fully annealed metals 
based upon the dislocation theory of their strength (see Cottrell 1949). The ~ 
value of p,, will vary considerably with the orientation difference between 
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the two grains adjoining the grain boundary. As a maximum value of 
Pg» We may take 1/2a (a—lattice constant), for we should not expect the 
number of dislocations in the grain boundary to exceed about one every 
two lattice planes. With a=3 x 108 cm. this gives P gp(Max.)=1-5 x 107/em. 
A minimum value for p,,, can be estimated from the expression given by 
Read and Shockley for the density of dislocations in a grain boundary 
defined by the angles 6 and ¢ (see fig. 2), 
Sin $6 
a 


{Cos ¢+Sin ¢}. 


Taking ¢=0 (p,, is not very sensitive to the value of ¢) and @~2° we 
find p,, (min)=1-3 x 108, 


Pgop=2 : 


IY TTT 
“Ene. 


Characterization of a grain boundary in terms of the angles @ and ¢. 
(Read and Shockley 1950.) 


Inserting these values of p,, and of p, into equation (3) and using the 
value of D’5/D=0-04 (Cu/Ni) we find, when 

Pop (max.)=1:5X107 Q—Q’=—2-7 K.cals/mole, 

pay (min.)=1-3X108  Q—Q’=+2-2 K.cals/mole. 

It is also readily shown that Q and Q’ would be equal if we supposed the 
grain boundary to have characteristic angles 6=6°, 4=0, which are not 
unlikely values. 

We see then that the measured ratio of D’S to D can be well accounted 
for on the assumption that Q and Q’ are equal, and therefore from equation 
(3) that the whole of the difference in the diffusion rates for grain boundary 
and for volume diffusion arises from a difference in the D, term in equation 
(3). But if the activation energies do differ, they will not do so by more 
than about 2-5 K. cals/mole. 
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} for T=1000° K. 


474 On the Grain Boundary Diffusion in Metals 


Similar calculations based on Ké’s model of the grain boundary are 
.mpossible in the absence of estimates of the grain boundary and volume 
densities of disordered groups. 

The experimental results at present available in the literature on the 
relative values of Q and Q’ are very few and conflicting. Measurements of 
the diffusion coefficient of molybdenum in tungsten by Van Liempt (1932) 
and of thorium in tungsten by Fonda et al. (Fonda, Young and Walker 
1933) show that D, increases with decreasing grain size, but that the 
activation energy is independent of the grain size of the tungsten, indicating 
the Q=Q’. On the other hand (Langmuir 1934) finds that Q—Q’ for 
the diffusion of thorium in tungsten is about 25 K. cals;mole and that 
D, for grain boundary diffusion is less than for volume diffusion. (See 
Barrer 1941.) 

On the other hand, we have seen that theoretical considerations lead us 
to expect that Q and Q’ might be equal and experimental work to confirm or 
disprove this would be very interesting. 

Another important line of work would be a study of the dependence of 
D’/D upon the relative orientations of the crystal faces forming the 
boundary. Any marked dependence would provide further evidence in 
favour of the “ transitional lattice ’ theory of the grain boundary. 

A study of grain boundary diffusion rates is then of importance, not only 
as a diffusion study but also in investigating other properties of grain 
boundaries. It is hoped that this new method of investigation, requiring 
only micrographical techniques, will facilitate such studies. 
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SUMMARY. 


The rectification of random narrow-band noise in the presence of a 
C.W. or modulated signal is analysed. The detector considered is of the 
type providing a rectified voltage which is a function of the instantaneous 
amplitude of the input wave. The smoothing present in the detector 
circuit is such that the high-frequency components of the applied wave 
are removed but the low-frequency variations of its envelope are faithfully 
transmitted. 

The mean (or D.C.) component of the rectified voltage and its r.m.s. 
fluctuation about that mean (or L.F. noise output) are calculated as a 
function of the input signal/noise ratio and particular attention is paid 
to linear and square-law detectors. 

The output signal/noise ratio for an amplitude-modulated signal input 
and for the audible beat reception of a C.W. signal are calculated. 

The spectrum of the L.F. noise output from a detector supplied with 
random noise is shown to be closely similar for linear and square-law 
detectors. 

The discrimination of a weak signal in noise is shown to be not critically 
dependent upon the law of the detector, and in particular the difference 
for a linear and a square-law detector is negligible. 

The effect of receiver bandwidth, meter time-constant and integration 
time in improving the discernment of a weak signal is considered. 


More Important SYMBOLS. 


A=instantaneous amplitude of the wave applied to the detector. 
B=bandwidth of input noise. 
E=r.m.s. input noise voltage. 

F(w)=input noise spectrum. 

G(w,,)=output noise spectrum. 

H=amplitude of heterodyne voltage. 
M=r.m.s. improvement due to meter time-constant. 
N=vumber of observations. 
Q=—detector discrimination. 

R(r)=an autocorrelation coefficient of the input noise. 
Samplitude of C.W. signal. 
T=time constant of meter. 


* Communication from the National Physical Laboratory. 
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T,—duration of period of integration. 
X, Y=“ in phase ” and “ quadrature ” amplitudes of noise voltage. 
f,,(@q)=normalized component spectra of the L.F. output noise. 
m=depth of signal modulation. 
n=exponent of detector law. 
u=instantaneous rectified voltage. 
v=—instantaneous input noise voltage. 
x—S8/4/(2)E=r.m.s. signal/noise ratio. 
y=A/r/(2)K. 
B=27B=angular bandwidth. 

é—instantaneous meter indication. 
7=time interval in autocorrelation. 

w,—angular frequency in L.F. output. 

w =—central angular frequency of input noise. 


| 
§1. INTRODUCTION. 


Tue calculation of the transmission of signal and noise voltages by a 
linear amplifier is a familiar and simple process for it is known that the 
output from an amplifier supplied with a mixture of signal and noise is 
the sum of the outputs if each were applied separately. However, this 
principle of superposition cannot be applied when the signal and noise 
pass through a non-linear circuit such as a rectifier or a limiter. 

The rectification of noise and of a mixture of noise and signal is of 
practical importance in relation to 


(i.) its effect on the signal/noise ratio in reception, 

(ii.) the measurement of noise, 
(iii.) the measurement of a signal in the presence of noise, and 
(iv.) the measurement of receiver sensitivity. 


The analysis of the rectification of noise presents difficulties which 
arise partly from the non-linear characteristics of the circuit and partly 
from the inherent characteristics of the noise which can only be specified 
statistically. (It should be stated at this stage that only fluctuation noise 
is considered in the present paper.) 

Franz (1940) considered the case of a detector whose current/voltage 
characteristic could be represented by a power series. The results of the 
analysis were not applicable to the case of the discontinuous characteristic 
of the linear detector nor did they refer to the practical condition in which 
the detector output follows the envelope of the applied wave rather than 
the instantaneous voltage. 

Franz (1941) later published an analysis of the rectification of noise 
by a detector with a discontinuous linear characteristic. In this paper 
the method of using a singular integral to represent the characteristic was 
successfully applied to derive the mean rectified voltage, and the spectrum 
of the low. frequency output. 


a as 


_- 
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Ragazzini (1942) considered the rectification of a mixture of a modulated 
signal and noise by an ideal linear detector in which the output voltage 
accurately follows the envelope of the applied wave. In order to make 
the analysis tractable it was found necessary to approximate in the 
evaluation of the components of the envelope and for this reason some of 
the equations are not exact. 

In the present paper it is proposed to analyse by elementary methods 
the rectification of a mixture of signal and noise by linear and square-law 
detectors and to compare the results so obtained with those appearing 
in the papers mentioned above. The practical application of the formule 
to the evaluation of the sensitivity of reception of modulated signals and 
of C.W. signals using a heterodyning oscillator is also discussed. 

It is convenient to draw attention at this point to the significance of 
the time constant of the detector circuit. Fig. 1 shows a typical 


Fig. 1. 


LINEAR AMPLIFIER DETECTOR (ENVELOPE FOLLOWER) 
BANDWIDTH 4 CENTRED ON Wo LAW: U = A" 


SIGNAL LOW-FREQUENCY 
NOISE @ A ‘, RECTIFIED OUTPUT 


INSTANTANEQUS AMPLITUDE A. 
NOISE SPECTRUM F (0) 


INSTANTANEOUS VOLTAGE UW 
NOISE SPECTRUM G (“a) 


Schematic circuit of a typical receiver. 


arrangement of a receiver ; the amplifier having a central frequency w, 
and bandwidth f passes the signal and noise to a detector having a load 
circuit RC across which the rectified voltage is developed. There are 
three possible cases of the time constant RC to be considered : 


‘ 1 
(a) RCK—. 
ial 
In this case there is practically no smoothing of the rectified voltage so 
that the latter consists of the positive half-waves of the applied voltage 
with the result that a considerable high-frequency component is present 
in the output and rectification is inefficient. This is the case considered 
by Franz, but it does not conform with the usual circuital practice. 
1 
tas eee 
In this case the high frequency components are filtered from the output 
and yet the time constant is small enough for the output to follow 
accurately the envelope of the input voltage. This is the case of practical 
importance in receiving technique. 


(ie ROS 
Wo 
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1 | 
= <nvU 
Seiad 

In this case the time constant is so large that the fluctuations of the 
envelope of the applied wave are smoothed out and the output voltage is 
steady. This case is relevant to the measurement of signal and noise by a 
diode voltmeter, and will be dealt with in a further paper. 

The law of the detector can be defined in terms of the relationship 
between the detector output voltage w and the instantaneous amplitude A 
of the input wave. In an ideal detector free from distortion the relation 
would be linear; in other cases, either through non-ideal behaviour in 
the detector or by design the dependence can be non-linear. In particular, 
for small applied voltages the rectified voltage from a diode detector 
tends to be proportional to the square of the amplitude while for large 
applied voltages it tends to equality with the amplitude. (See 
Appendix IV.). From the point of view of the present analysis there is 
no essential difference between a linear detector followed by a non-linear 
aperiodic amplifier and a non-linear detector. 

The detector characteristic considered in this paper is assumed to have 
the general form 

UA nS yee ca 


It has been assumed for simplicity (but with only trivial loss of generality) 
that the constant of proportionality is unity. 


§2. SraTIsTIcCAL PROPERTIES OF A MrxTURE OF A C.W. SIGNAL AND 
Narrow-Banpd NOISE. 


Since fluctuation noise arises from a succession of random and 
independent elementary events it obeys certain statistical distributions 
which can be conveniently specified in terms of the r.m.s. value E of the 
noise voltage, defined by 


ete 
F?=. Lim erty. 
T+>o Le 0 
where v(¢) is the instantaneous noise voltage. 

Since the amplifier supplying the noise is assumed to have a narrow 
band * B centred on wo, the noise at the input to the rectifier can be 
written 

v(t)= X(t) cos wot-+ Y(t) sin wot CO eS) 


in which X and Y are relatively slowly varying functions of time, and 
can be regarded as the ‘“‘in phase” and ‘“ quadrature” amplitudes of 
noise relative to a C.W. signal A cos wo. 


* For the purpose of setting a limit on the term “ narrow band” it may 
be arbitrarily defined as applying to the case in which the bandwidth is not 
greater than one-fifth of the central frequency This sort of restriction is 


necessary for the concept of “ instantaneous amplitude” to be applicable to 
the noise wave. 


Cars 
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Now v, X and Y have a normal or Gaussian distribution about zero ; 
thus for », 


P(v) . dv= Wee sees. ea) 


I 
Vn). Bo? be om): 


From this distribution the following averages may be formed : 


2— X2— Y2— E2 
XY=0, 
[el =[X]=[Y]=veme. 
The amplitude or envelope of the fluctuation noise voltage is given by 
the positive value of 
r=/(X?2+ Y?) 
and its distribution function is of the Rayleigh type 
r v2 
P(r) -dr= paexp (— sm) CO) SE ta (D) 


from which 
722K? and r= 4/(7/2)K. 
If now a sinusoidal signal S cos wot is added to the noise the r.m.s. 
signal/noise voltage ratio is given by 
Sas 
V/(2)E 


or 20 log,)x in decibels. 


The instantaneous amplitude A of the mixture is given by 


eV eee Yo [tne aeneb os 3. * (8) 
whose distribution is readily shown (Appendix I.) to be given by 
A A?-+-$? AS 
P(A); dA= 7eexP (- ae ibe (Fs): Ame = (1) 


I, being the modified Bessel Function of the first kind, zero order. Le 
the parameter 


Soe: 
Y= TRE 
is introduced we may write (7) in the form 
P(y) . dy=2y exp (—y?—2?) . Iy(2yx) . dy. are, (oS) 


For weak signals (x<1) the amplitude distribution of the mixture is 
still of the Rayleigh type 


9 2 
P(y)  dy= sgex (— 5) a, Be eee (9) 
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which shows that the addition of the weak signal effectively increases 
the r.m.s. noise voltage from E to E(1+-2?)?, which is equal to (H?-+-- 48?) # 
the r.m.s. voltage of the mixture. This property will be found useful 
later in considering signals below noise level. 

In the case of large signals (v>1) we have from the asymptotic 
expansion of I, ? 


Py). dywv/(y/mx) exp[—(y—2)*].dy  . . . « (10) 


which shows that the amplitude distribution of the mixture is centred 
on the signal amplitude and is approximately Gaussian with r.m.s. 
fluctuation equal to the r.m.s. noise voltage E. 


§3. THe Mran ReEcTIFIED VOLTAGE. 


If the voltage output V from the rectifier is related to the amplitude A 
of the applied voltage by the nth-power law, the d.c. or mean rectified 
voltage which would be indicated by a voltmeter of long time constant is 


mS | Awp(a) dA 
0 


I 


2(4/(2)E)” ike y"*1 exp (—y’—2*)I,(2yx) . du, 


0 


where, as above, the signal/noise ratio is denoted by x. This integral can 
be evaluated by a general relation due to Hankel which leads to 


w= A"=(y(2)B)" (5) B(—5.h —a), A, cate 


in which ,F, is the hypergeometric function expressible by the Kummer 
series 


mY, (-92)(—2?) (= 2) —n/2) (— 02)? 
1B, ( gels at) =14 iaghs Diet 12 Smo 
(This function is denoted by M in Jahnke—Emde’s Tables.) 

It is clear that if m is a positive even integer the series terminates 


and the average of A” can be simply expressed. Thus for a square law 
rectifier n=2 for which 


At—28(1-Px)y ee 


showing that the mean rectified voltage is accurately proportional to the 
mean square values of the applied voltage, as in the case of the 
thermocouple for measurement of coexistent signal and noise. The 
anode-bend and the diode rectifier approximate to this condition for 
small input voltages. 

The mean square output voltage is obtained in the general case by 
replacing by 2n in equation (11): 


U2= A — (22) | iF \(—n, 1, —2?). otters obhe) 
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A case of practical importance is the linear rectifier (2=1) for which 
Ree (2M (dl a). cs dieu, (14) 


Fig. 2 shows A| 4/(77/2)K as a function of the signal/noise ratio x. For 
a small signal/noise ratio (v<1) the Kummer series gives. 


= 4 6 
A=4/(n/2)E (14424 - ¥ a ae :) 


==4/((7/2)A(4-a*)| approx, . . . . «+. « (15) 


showing that the linear rectifier tends to give an indication proportional 
to the r.m.s. value of the mixture of signal and noise when the latter 
predominates. 


Fig 2 


en Ma Gare nee so) 


— = RMS. NOISE VOLTAGE 
A = AVERAGE OF ENVELOPE 


R.M.S. SIGNAL / NOISE VOLTAGE RATIO 


Average amplitude of the envelope of a mixture of signal and noise 
(as would be indicated by a linear rectifier). 


For a large signal/noise ratio (v>1) the asymptotic expansion of the 
hypergeometric function 


17 n2 
ee er os) see 


gives 


A=vie)Ee(1+Z5). Sir ea eT) 
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In the above it has been assumed that the time constant of the rectifier 
circuit is such that the R.F. or I.F. is removed from the output, but 
sufficiently small for the envelope to be accurately followed, the long 
time constant of the meter providing the averaging. In Franz’s analysis 
(1941) of the linear rectifier, however, it was assumed that no smoothing was 
present, so that the rectified output consisted of the positive half-waves 
of the applied voltage. Since the average positive voltage is related to 
the average amplitude by 

Tie Asin 8. d0==, | ed 


7) 0 


it is seen that Franz’s rectified output should be 1/7 of the value calculated 
above, 2. é. 


7, — — iF (eh 1, —-2)a i ee 


Franz (1941) obtained the average in the form 


Les INE ; ; 


It may be shown that equetions (20) and (19) are identical by using 
Kummer’s relation for the Bessel function and the recurrence formule 
tor EF.” For: 


e-#'($0%) =F, (5, 1, —2%), | 
e~ #1 9(32") +1, (32")]=1F (3, 2, —2*), 
FAs, 2, — 2?) + Fi, ], —x?)=,F,(—3, 1, a); 


(21) 


§4. Mran Square VOLTAGE OF THE OvuTPuUT NOISE. 


The noise output from a linear rectifier which accurately follows the 
envelope of the applied voltage can be expressed in terms of the 
fluctuation of the envelope about its average value. Thus the mean square 
noise output is given by 


wR (Aye=et| oat) 3 |, ee heh 


where F denotes the hypergeometric function in equation (14). 
If the signal/noise ratio is small («<1) the noise output is 


Oi (2-5) B+ a2) <0-43(1 4 29, ad. ea aod 


which shows that the mean square noise output is proportional to the 
total mean square value of the applied signal and noise. 


For large signal/noise ratios (x1) the output noise tends to the limit 


given by 
Ug tehat.. i SSRN eee em 24) 


which shows that in the presence of a large C.W. signal the r.m.s. values 
of the input and output noise are equal. 
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To illustrate equation (22) the ratio w,/E has been plotted as a 
function of x in fig. 3. It is of interest to compare this evaluation with 
the approximate analysis by Ragazzini (1942) who finds the simpler 


but less accurate value 
i 
7} —— = i? 
Uaae He || Mie howe ten of ot OO) 


which gives }K? instead of 0-43E? for the output in the absence of a 
signal, but tends to the same upper limit of E? as above. For comparison 
Ragazzini’s equation is plotted as a dashed curve in fig. 3, and it is seen 
to give a very good approximation to the exact value of x>1. 


Fig 3 
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EQUATION 22 OF TEXT 
——=—=— RAGAZZINI'S EQUATION 


R.M.S. OUTPUT NOISE/R.M.S. INPUT NOISE 


o Os Lo =i rs 20 2u5 30 
INPUT SIGNAL /NOISE RATIO 


Response of a linear rectifier to a mixture of noise and C.W. signal. 


For the linear rectification of noise in the absence of a signal and 
with no smoothing Franz (1941) found for the mean square output noise 
voltage ‘ 


Now in this series the ratio of the (r+1)th term to the rth term is 
(r—4)2/(r-+1)? showing that the series is hypergeometric and is in fact 
equal to 4[F(—4, —4, 1, 1)—1] which by Gauss’s theorem= 4[(4/r)—1] 
(Franz gave 25/2816 instead of the correct value 25/4096 for the fourth 


term). Thus Frianz’s expression becomes 


1 7 
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which is 1/7? of the value given by equation (23) for the output from an 
envelope follower with no signal. This reduction factor occurs since in 
Franz’s case of no smoothing the noise output arises from fluctuations in 
the average of the positive voltage (0,,.) which was seen in the last section 
to be 1/7 times the corresponding average of the envelope A. 

The square law detector can be considered in a manner similar to 
that used for a linear detector. In this case the mean square noise output 
is given by 


w=w—m=[At—(A?)?] 
—E4 E (: 4202+ >) —4(1 +27)| 
=—454(1-1 29%) —4H2(E2+8%), 2. . . . . . . (26) 


Fig 4 


Uno = R.MS. NOISE OUTPUT VOLTAGE WITHOUT CARRIER 
Un = RMS. NOISE OUTPUT VOLTAGE WITH CARRIER 


SIGNAL / NOISE RATIO 


Increase of noise output from a linear and from a square-law rectifier on the 
introduction of a C.W. signal. 


which shows that the noise output rises uniformly as the signal is increased 
instead of tending to a limit as in the case of the linear rectifier. The ratio 
of the r.m.s. output noise voltage with signal present to that with signal 
absent is shown in fig. 4 as a function of the signal/noise 2. 


$5. FREQUENCY SPECTRUM OF THE OuTPUT NOISE. 

The formule given in the last section for u? are for the total mean 
square output noise voltage ; in point of fact the noise is spread out in a 
continuous spectrum whose integral is equal to u2. The form of this 
spectrum is of importance in determining the effect which post-detector 
selectivity will have on the noise voltage and on its ratio to a signal. 


AG Ee hy 


nag 
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The frequency distribution of the output noise can be conveniently 
expressed in terms of the spectrum G(w,) where w, is the audio 
frequency. This spectrum is such that the mean square noise voltage in 
a frequency interval dw, is given by G(w,,) dw, with 


0 
iE Gay Wee te ROE 7) 

Clearly the spectrum of the output noise is a function of the spectrum 
applied to the input of the rectifier ; since the noise impulses at their 
origin are extremely short, the input spectrum is determined by the 
frequency characteristics of the amplifier between the noise source and 
the rectifier. The rate at which the envelope of the noise varies from 
moment to moment is inversely dependent upon the bandwidth of the 
amplifier since the narrower the bandwidth the greater the lengthening of the 
impulses and hence the closer the correlation at successive instants of 
time. The correlation function of a random disturbance and its power 
spectrum are connected by the Fourier cosine transform. 

The frequency spectrum of the input noise to the detector denoted 
by F(w) is defined in a manner similar to that used above for the output 
noise, namely, that the mean square voltage in a frequency interval 
dw is F(w).dw. Then 


| F) .do=E», ree... (28) 
0 


where E? is the mean square noise voltage. 
If w, is the mid-band frequency, it is possible to specify an effective 
bandwidth £ for the noise at the input of the rectifier : 


[ Flo). de pe | 
ee gy. «3 (29) 
B Fw) F(w9) ( 
Consider first the response of the square-law rectifier to a mixture of 
noise and signal. If the rectifier characteristic is given by 


u= A2 


it is seen that when two sinusoidal voltages A, sin w,t and Ay, sin wet are 
applied, the rectified voltage can be written 


u—[A?+ A2+2A,A, cos (wy—wy)t], 


which has a pure modulation product of amplitude 2A,A,. Hence, if 
any number of sinusoidal voltages are applied to the square-law rectifier, 
the low frequency products will be linearly superposed, so that the 
components of w have the general form 

u=2 2 A,A, cos (w,—a,)t. Me ae ee (SU) 


Uae an 
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The mean square output voltage thus has the form 


Wins? 2 AAS ee ee 
y & 
In these summations the products are only to occur once so that A,A, 
is taken as being included in A,A,. 

If now fluctuation noise of spectrum F(w) and a sinusoidal signal of 
amplitude S and frequency w, are applied to the rectifier the mean 
square noise output in an interval w, to w,+dw, has two components— 
that arising from the cross-products of the noise terms and given by 


2 | Ona) Dace te dercdar 
0 


and that arising from the product of the signal and the noise : 
28°[2F (wo+w_) +2F(w 9—w,) Jder,. 


Thus the output noise spectrum is given by 
G(w,)=4 [2] F(w)F(w+a,) « dea + SHE (cop beg) + Foo ea) | : 
0 


For simplicity it will be assumed that the input noise spectrum is 
symmetrical about wy, and thus 


Glo)=8| |" Fe)Fa+e,) .da+S*F(oy bo) |. -. (82) 


This formula will be applied to a rectangular input noise spectrum 
corresponding to an ideal band-pass characteristic. 

If the pass band is f it is seen from equation (29) that within this band 
the input spectrum is given by 


Bi(w)s=HA Be 7 ee eis eee (33) 
and is zero outside it. 
Hence equation (32) gives for the output spectrum 


wy , 222 1 
G(w,,)=8E* [= 3 + 7 for Wa Se > | 
Glo.) =8E for’ <w, <P, | oa nt ee) 
G(w,)=0 for B<a,. J 


The integration of G(w,) gives 
eaene “B 
Un= j} G(@,) . deog= | G(w,) . dw, =4E4(1+ 22?) 
0 


in accordance with equation (26) for the mean square noise output from 
a square-law rectifier. 
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Fig. 5 shows the output spectrum for a number of values of signal/noise 
ratio. It is seen that in the absence of a signal the spectrum is triangular 
and of width f at the base which is twice as wide as the input spectrum 
about its central frequency w,. As the signal increases the output spectrum 
progressively tends to coincide with the bisected input spectrum. These 
results can be easily understood from a physical point of view by 


Fig 5 
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X = SIGNAL/NOISE RATIO 
AT RECTIFIER INPUT 
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Wa 
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Output noise spectrum from a square-law rectifier when the input 
spectrum is rectangular. 


consideration of the relative contribution to the output by the various 
cross-products of the noise components with each other and with the 
signal. 

The results derived above apply to the square-law detector and are of 
practical application, since for small values of applied voltage detectors 
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tend to behave as square-law rather than linear devices. It is of interest, 
however, to know how the output spectrum from the linear detector 
differs from that from the square-law detector ; the difference between 
the output noise spectra of the two detectors is least when a large 
sinusoidal signal is present, since the input spectrum is faithfully 
reproduced. On the other hand, if the detector is supplied with noise 
having a rectangular spectrum and no signal is present, then, as is shown 
in §8, the output spectrum of the linear detector is very little different 
from the triangular spectrum given by the square-law detector. Therefore, 
a good approximation to the frequency distribution is obtained if 
equation (22) is used for the total output noise instead of equation (26) 
and the output spectrum shape is assumed to be the same as for the 
square-law detector. 


§6. Sr@naL/NoIsE Ratio FOR RECEPTION OF MODULATED WAVES. 
If the signal of carrier frequency w,g and amplitude Sp is modulated 
at frequency w,, to a depth m, the wave can be written 


So(1--m cos w,,t) COS wot 
m 
=Sy E wot+ = 008 (wy—W» t+ F 008 (obo) | ; 


whose r.m.s. value is $83(1-+-4m?). 
It will be assumed that the input spectrum is rectangular and of 
width 8 and the detector is square-law, and furthermore, that 


IDS Ses 


since otherwise the sidebands would not be passed to the detector. 
The fundamental component of the rectified output at frequency w,, 
has the form 


2 
2mSp COS w,pt 
the mean square value of which is 
uz—=2m?8$—=8m2a7iEt 


taking x)=S)/1/(2)E, the r.m.s. carrier/noise ratio. The ratio of the 
r.m.s. values of the modulated signal to the noise is then a=a)(1-+-4m?)}. 
Now the mean square noise output due to the components of the noise 
mixing with each other and with the carrier is given by the value already 
found : 
4H4(1+ 22°). 


The mean square noise output due to the components of the noise mixing 
with the two sidebands is easily shown to be 


4 E4m?2?. 
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The total mean square noise output is therefore 
a GBA (Da mtae toy (35) 
Thus the output signal/noise ratio is given by 
Uy maz 
U, VRP aR yn)] 
which is the ratio of the r.m.s. audio output voltage of frequency w,, 


to the total r.m.s. output noise voltage. Fig. 6 shows the variation of 
u,/u,, with the carrier/noise ratio x) at the input for m=1, 0-5 and 0-2. 


(36) 


Fig. 6. 
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m = DEPTH OF MODULATION 
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-20 -10 fo) 10 20 
INPUT SIGNAL CARRIER/NOISE RATIO (2g) 
IN dbs 


Relation between the signal/noise ratio at the output and input of a 
square-law detector for an amplitude modulated signal. 


‘Even the best case (m=1, 2» large) the loss in signal/noise ratio is 
1-8 db., but it can be much greater than this in conditions of small input 
signal/noise ratio x, since then w,/u, 1s proportional to oe This 
suppression of the signal by the noise when the latter preponderates is 
clearly seen in fig. 6 as a change of slope as x, decreases through Ce) (0 db). 

In this analysis any contribution of the second harmonic 2w,, to the 
useful signal output has been ignored, since its inclusion would only 
increase u, by a factor of (1-+-m?/16)? which is always negligible. 


o 
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§7. SrGNAL/NoISE Rat1Io FOR HETERODYNE RECEPTION OF C.W. Siena. 


A commonly used method of reception of a C.W. signal consists in 
producing an audible beat tone by heterodyning the signal at the detector 
by means of a local oscillator. Let the local oscillation be of amplitude H 
and the signal of amplitude S and let their frequencies lie with the 
rectangular pass-band (w)—f/2) to (w»+f/2) of the amplifier preceding 
the detector. 

The amplitude of the beat note produced by a square-law detector is 
2 HS and so the mean square signal output voltage is 


i= 2S? eee eer 


The mean square noise output voltage can easily be derived from the 
results of the earlier sections : 


. u2—41(B?-H2--S*)) ees 
Thus putting 
eee 
V(2)E° 


the output signal/noise ratio is found to be 


CP hx 
u, = V/ (a+) . eo. ro (39) 

Fig. 7 shows the signal/noise ratios at the input and at the output of 
the square-law detector as a function of the heterodying voltage. If this 
voltage is large compared. with both the noise and signal voltages, the 
signal/noise ratio is seen to be unimpaired by the detector ; in fact the 
detector behaves as a linear frequency changer and leaves the spectra 
and the relative magnitudes of the signal and noise components 
unaffected apart from a common translation in frequency. 

This equality of signal/noise ratio before and after detection for large 
values of heterodyning voltage also applies to the linear detector; in 
fact for any type of detector the transfer ratio will tend to unity if the 
heterodyning voltage is made sufficiently large. 

These remarks apply to the frequency-changing stage of a super- 
heterodyne as well as to the audible beat method of C.W. reception. 
However, it should be remembered that the present analysis is not 
concerned with the effect of any noise introduced by the heterodyning 
detector itself. 


§8. THE Response or A Detector to Norse or To NoIse WITH A 
Weak SIGNAL. 


The preceding sections have dealt with the detection of a mixture of 
signal and noise with particular emphasis on square-law detectors: and 
a general ratio of signal to noise. However, when the signal is weak 
compared with the noise as is the practical case in many problems of 
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marginal reception we have seen (equation (9)) that the amplitude of the 
mixture is statistically similar to that of random noise having a mean 
Square value equal to the sum of the mean square values of the signal and 
the noise. Thus results may be derived for the output from an nth power 
envelope detector which are reasonably accurate and indicate the 
behaviour of such a detector in the observation of signals below noise level. 

In order to derive the output spectrum of a detector to which random 
noise has been applied it is quickest to calculate the autocorrelation 
function of the output voltage and then form its cosine transform. Now 
in the nth power detector the output voltage is given by 


u=A” 


Higa: 


_ R.M.S. HETERODYNE VOLTAGE H 


“RMS. NOISE VOLTAGE 


Us 
Un 


OUTPUT SIGNAL/NOISE RATIO 


INPUT SIGNAL / NOISE RATIO 


Relation between the signal/noise ratio at the input and at the output 
of a square-law detector for heterodyne reception of a C.W. signal. 


and hence the autocorrelation function 
Ug A,"A,” 


where the suffices 1 and 2 refer to instants of time separated by an 
interval 7. Then it is shown in Appendix 2 that 


nN nN 
=H En) 1? F (— 5, = ark R’), - + (40) 


2L2 
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where F is the hypergeometric function and R is a function of 7 having 
the form 
sin $Br 
R()=—, . . --.- = « = i 
=, (41) 
in the case of a rectangular input noise spectrum of width f. 
The audio output spectrum then has the form 


= Sa fice 
Gum | (tigi) cos 17) dr 
7/0 
Now since %%, is a series in R* it may be written 
UU,g=A)+a,R?+a,.R4+ ne ele . . . . . (4 3) 


which will terminate if n is an even integer and be an infinite series 
otherwise. Each term in this series can be regarded as producing its own 
component spectrum 


2 rae i 1 2m 
G(w,)=- | 2: Om (= -) COS w,T . dt 


TJQ 1 2Br 
= 2nd am(a)s ot oe PO eee ee eee 
where 
2° /sin 4Br\? 
flo = ( TBr ) 008 to, Ts OT. feo cee eT 


is the pth normalized component spectrum each of which is a function 
of w,/B alone and has unit content, 7. e. 


| Solera) Bes aD 


The properties of these spectra have been investigated by Franz (1940) 
and some of these are given in Appendix ITI. 
From equation (43) the d.c. component of the output voltage is 


U=Up=V (Go); ue » Sow @ 2 48) 
while the mean square output voltage is 
w—=a)ta,+a,+ eer a erence . . . . . . (47) 


and the mean square L.F. output voltage is therefore 


w—uw—=a,tant.... =| Glo.) dae kr. eeeraes 
0 
For an nth law detector 


Uo (4 (2) H*(tn) lo ee ee roy 


actly cease Mare 


——— 
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and the mean square noise output is hence 
u®—ua—[n 1—(4n) 12(2E", . . |... (5D) 


The last three relations are consistent with equations (11) and (13) 
derived earlier and can in fact be obtained directly from the Rayleigh 
distribution. The spectral composition of u will now be considered in 
more detail. The low frequency limit of the L.F. spectrum is given by 


G(0)= Say fan(O)- 


This must not be confused with the d.c. component of the detector output 
which corresponds to a spectral “ line ” at zero frequency of infinitesimal 
width but finite content. 

The square-law detector will be discussed briefly first, since it has 
already been dealt with in an elementary fashion in §5. For this case 


Ogee SET) eo cele fet Se) ~ (52) 
whence the mean square L.F. noise voltage output is 
peat) od Ve ed OA 8... (53) 
and the spectrum therefore has only one component : 
8H4 
Cle.) =ABYf,(o.)= (1-2), Pee res) 


which agrees with the earlier equation (34) when the signal is put equal 
to zero (x=0). 

The effect of varying the bandwidth 6 of the predetector circuits on 
the detector output spectrum is seen to be to increase both the height 
and the width of spectrum proportionally to 8 since E?/8=F (the input 
noise spectral density) is assumed to be constant. 

It is noted from equations (52) and (53) that the autocorrelation 
coefficient of the output of a square-law detector (i.e. of the square of 
amplitude of a noise voltage) is given simply by 

NET 
px(7=22—"=R. 2 1... (5) 
w—wuP 
‘In the linear envelope detector (n=1) the output voltage is equal to 
the envelope of the noise current and its autocorrelation function is 


7 


sn fog ob R? Rt BR ) hey ued afi) 
=75 (4 y+ tat (56) 


The mean square L.F. noise output voltage is 


w= 2 BYR(—}, —4, 1, 1)—1)= (2-5) B -. (87) 
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which is consistent with the results of § 4. : 
Now. the spectrum of the L.F. noise output is obtained from 
equations (45) and (56) 
- 7, &(m—3/2) !7 
G(o,)= 2A nd 2m= 2 ee mit Fama). Te 622 (58) 
Since the content of the first spectral component a,f, is 16 times 
that of the following component a,f,, (for a,=a,/16 and all the Fn have 
unit content) it is seen that this component provides a close approximation 
to the noise output spectrum 


7 7 2 Wa 
Gla) S 3 Eee) = TE; (1— =) . . ° e (59) 
Fig. 8. 
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Spectral distribution of noise output from linear and square-law detectors 
when input noise has rectangular spectrum of width f. 


up to 6, and zero beyond. This is in fact of just the triangular form 
which is realized exactly with a square-law detector. The higher order 
component spectra represent the small departure from the triangular 
shape but their combined content is under 8 per cent of the total L.F. 
noise given by (57). 

One point of difference from the square-law detector is that in the 
linear detector, when the bandwidth 8B of the input noise is increased 
(leaving its spectral density E?/8 constant), the height of the spectrum 
does not change although its width increases proportionately with £. 
Fig. 8 illustrates the very close similarity in form of the output noise 
spectra from the linear and square-law detector. 


——— 
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From equation (56) it is seen that the autocorrelation coefficient of 
the envelope of a noise voltage is given by 
Ute — UR. at amy cae 
The numerical factor 7/4(4—7)=0-91 and it is seen that the autocorrelation 
coefficient of the amplitude will always lie between 0-91 R? (for R small, 
7 large) and R? (R nearly unity, 7 small), and is therefore only slightly 
different from that of the square of the amplitude. 


§ 9. SrenaL/NoIsE DiscRIMINATION OF A DETECTOR. 


In a receiver used for the observation or measurement of weak signals 
the law of the detector plays a part in determining the discrimination 
against noise since upon it depends both the incremental output for a 
small applied signal and the fluctuations of output due to the applied noise. 

In the case of meter indication a weak signal would be observed by 
first noting the reading in the absence, and then in the presence of the 
signal, the difference of reading providing a measure of the signal 
intensity. The initial reading due to the noise can be “ backed off” 
(i.e. compensated electrically or mechanically) and the sensitivity of 
measurement so increased since a more sensitive meter may now be used. 
The limit to this process is set by the fluctuations of the meter readings 
about their average positions arising from the noise. These fluctuations 
may be reduced by the use of a meter of long time constant but 
correspondingly the duration of a measurement is prolonged ; the effect 
of the meter time constant is discussed in the next section. 

The signal/noise discrimination of a detector can be objectively 
expressed in terms of the ratio of the increase of the average (or d.c.) 
output on applying a certain signal to the r.m.s. fluctuation of the 
output (¢.e. the L.F. output) about the average. The effect of long-term 
observation and averaging on the discrimination is considered in the 
next section. 

Usually in the condition of limiting sensitivity the signal/noise ratio 
at the detector is small and this enables a simpler parameter to be obtained 
to express the discrimination of a detector in terms of its law. If the 
noise of r.m.s. voltage E is accompanied by a signal S cos wot and is applied 
to an nth power-law detector the d.c. output voltage is given by 


where Te baer i) ee iy hs -3*) Mere Bie, hae LA ) 
* 4/E 


is the signal/noise ratio as before. wie 
The total mean-square output voltage from the detector 1s given by 


ur2=2"E2"y | F\(—n, 1, —2?). 
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For the signal/noise ratio « small the mean-square fluctuation of the 
output will be given approximately by 


u2=(2E2)" E = (5) | 


as already given by equation (51). 
The increase of d.c. output voltage due to the signal is found from 
equation (11) to be 


U(X) —Up(0)=[ V/ (2) E]"(32) 5a ol fuente erent? 


for « small. It is noted that the detector behaves as a square-law device 
as regards the increase of output due to a weak signal owing to the 
suppression effect of the stronger noise. 

A convenient parameter for expressing the detector discrimination for 
very small signals is thas : 


QS es a 


Uc" n! why 3 
(nj2) 12 


It is found that Q has a maximum of unity at n=2 but this maximum 
is very flat, so that, in practice, superiority of a square-law detector in 
this respect is of no significance. For example: 


n=}, Q=0-90 
nai) | G220-96 
n=2, Q=1; 


n=4, Q=0:89. 
A greater practical advantage of the square-law detector is that the 
d.c. output for a signal and noise input is given by 


Up =2E*(1 +2?) =2(E2+ 48%), 
showing that whatever the signal/noise ratio x, a change of indication 
can be directly interpreted as a change of mean-square input voltage. 


§10. INFLUENCE OF THE Meter Time Constant AND INTEGRATION 
TIME ON SENSITIVITY. 


If the time constant of the indicating instrument connected to the 
detector output be increased the disturbing effect of the fluctuation of 
indication is reduced, so giving enhanced sensitivity at the expense of 
speed of measurement, or in statistical language, we would say that the 
random variations of the measured quantity are smoothed out by 
sampling over a longer time. 

If the time constant of the meter is T the mean-square fluctuation of 
the output voltage is reduced in the ratio 

© G(w,) . dw, 
ldo labo Ale) 


Sr = verv nearl 
M? ue Tut y 
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if the receiver bandwidth is large compared with the reciprocal of the 
meter time constant. It is assumed that the signal is small compared with 
the noise so that the output spectrum is substantially the same as that 
found earlier for a noise input only. 

The improvement factor due to the meter time constant is thus 


2Tu2 
m= /[=3 |. Mee eee AVE TG 3) 


which for a square-law detector is found from equations (53) and (54) to be 
M,=-v (BT/7)=~+/(2BT). Me ee a OL) 

While for a linear detector from equations (57) and (58) 
Mr== 10347 (2B Tandon eit? erin 2 (68) 


which is negligibly different from the value for the square-law detector. 
Thus combining these results with those of the previous section the 
ratio of the increase of indication to the r.m.s. fluctuation is 


Up (x) —Uo(0) 


MP OM=eVeBT) . . . . (66) 


to a close approximation for hoth square-law and linear detectors. 

Since the signal/noise power ratio x? is inversely proportional to the 
receiver bandwidth B, the discrimination y is seen to be proportional 
to 4/(T/B). Hence, halving the receiver bandwidth or doubling the 
meter time constant is equivalent to a 1-5 db increase of signal voltage. 

Another method of improving the sensitivity of measurement of a weak 
signal in noise is by taking the average of meter readings over a time T, 
long compared with the meter time constant T and so extending the period 
of integration. Clearly it is pointless to take readings too frequently 
in a given period since the closeness of correlation between successive 
values will prevent any further improvement in accuracy being obtained. 

Consider the case in which N readings of the meter indication @ are 
taken. Their average is 


N 
pee ery 


1 
NRT 


and the standard deviation o of this averaye due to the random 
fluctuations in @ about its true average 4) is given by 


=| 2 (0,00) | 


Noa 
= ga[ Not+2ot 2 w—pr, | ae) cage (67) 
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where o, is the standard deviation of a single reading and 1, is the 
correlation coefficient between readings spaced » apart. In the case of 
readings evenly spaced by a time s and a meter time constant T 


r,= exp (—ps/T) 

and 

ox VIN =rj)—2r7, 01-77) 

Go) N(I—r,) 
where 7,= exp (—s/T) is the correlation coefficient between successive 
readings. Naturally when this tends to unity corresponding to readings 
taken in very close succession oy=c, and there is no improvement due 
to averaging. 

At the other limit ox tends to o,/1/(N) as the correlation tends to zero 

and this factor ,/(N) represents the maximum improvement obtainable. 
Fig. 9 shows the smoothing factor as a function of N for various degrees 


Pe Pap ees Geog eS) 


Fig. 9. 
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NUMBER OF READINGS AVERAGED N 
Reduction of fluctuation by averaging N regularly spaced readings. 


of correlation. However, it is more useful to know what improvement 
can be obtained in a given observing time T, than with a given number of 
observations. 

Imagine an integrating circuit which at any moment continuously 
provides the average value of the indications over the previous period of 
length 'Fo, that is, it forms the quantity 


l “to 
gata | 0. dt 
to—To 


at the time ¢) of observation. The fluctuations in z can be catculated 
directly or as a limit to the summation given by equation (67) when 
discrete observations are made at infinitesimal spacings. It is found 
that the mean square deviation in z is given by 

2 -T 

on, = T2 (T—t)r(t) 3 dt, 5 . ° . . . (69) 
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which for an exponential correlation r(t)= exp (—t/T) of indications gives 


onlays? iat [ly ay it 
= fa (+ gow (TwD- er) |, ope 


which is shown in fig. 10. 

Now from equations (63), (64) and (65), o,=w,/M, and thus the overall 
improvement due to the meter time constant T and continuous averaging 
over Ty which is assumed large compared with T is given by 


Un To 
sa aN 0) = 
a J Gr) 4/ (BT). ee Se tL) 
Fig. 10. 
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Reduction of fluctuation by continuous averaging over time T». 
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AVELP WeNeD dT X51: 
PROBABILITY DISTRIBUTION OF C.W. StanaL Pius Narrow-Banb 
NoIseE. 
The instantaneous amplitude A of the mixture of signal and noise is 
given by 
A?=(X--S)?+ Y?. 


500 R. E. Burgess on the Rectification and 


Now (X-+8) and Y have independent normal distributions each of 
standard deviation E and centred on S and zero respectively. Putting 
é=(X-+S)?, n=Y? and «=A?=£-+7, we see that 


P(g)= sEaap e | - vw | + exp | — Eelees \ 


l E+8? SV/(é) 
~ E/Qz€) exp(—§ 2H? ) os E 


whence 
1 


* 1 
Peo)= arms |, Rae? (aE) cosh 
= or (= Se) i cosh eo] . a0 


1 a+S? (2) 
= sex (— ser)! »| a 


Thus the probability distribution for the amplitude A=«? is 


A A?--S? SA\ 


which is equation (7) of the text. 


Aree ND oe 


AUTOCORRELATION FUNCTION FOR THE NTH POWER OF THE 
Noise AMPLITUDE. 


If V14=X, C08 wot; + Y, sin wot, 
Vg= Xo COS Wola +Ye SiN wol, 
are the noise voltages at times Sencar by (¢.—t,)=7, their autocorrelation is 


E2r(r)=0,0,=}(X,Xo+¥ 1 Y_) COs wor 


when the spectrum is symmetrical, 7. e. X,Y,=X,Y,=0. In the case 
of a rectangular spectrum F(w)=E?/B and 


1 "Wo +38 sin 4Br 
r(7)= 3. es gi wt .dw= TBs COS WoT. 
Thus 
XX, ok we a <FR(z), 
where BGs aBr 
tBr 


is the “ envelope ” of the correlation function r(r). 


—S 
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The joint distribution of X,Y,X,Y, is thus the product of two bivariate 
normal distributions 


P(X, Ya: xe Y,) 


ul I oxp | — eit X-2RX,X.+V2+ V3—-ORY,Y, 
~ 472 R41 —R) XP 2H2(1—R2) ‘I 


Now we require the autocorrelation of 


An = (X24 Y2)n2 
which is given by 


MAR ArAg= | i { hin (X24 Y?)""(X24 -V2"2P(X., V,, X,, Y,)dX, dV, dX, dVp. 


In order to evaluate this, a convenient substitution is 
X7=qe"’E?(1 —R2) cos? ¢,, Y?=qe?E?(1—R2) sin? ¢,, 
X3=qe-?E?(1—R?) cos? ¢,, Y3=qe-"E?(1— R?) sin? do, 


which for the transformation to the variables p, g, 4, and ¢, has the 
Jacobian 


J |=gE(1—R2), 
whence 
—— #2] —R2)r+1 2 
At At — eo iN aq dp | dd, ie dps . qr 
x exp [—q cosh p—q R cos ($,—¢4))] 
—F2"(] — R201 | g"*1K o(q)Iq(qR) dq 
0 


—(22)" Ran (5 5)iee (Sp oe —R') 


Lapog ved LEV ical mae, _R2 
=n (5) 8 ae oe 


which is equation (40) of the text. 


A ELPAN Dt Xx = LL 


THe COMPONENTS OF THE OuTPUT SPECTRUM DUE TO A 
RECTANGULAR INPUT SPECTRUM. 


The component spectra of the detector output voltage are defined by 
2 (° /sin 48r\? 
f(@Od= =) ( a ) COS wat . dt. 


= 


502 R. E. Burgess on the Rectification and 
They have the normalizing property 
Ir Sp(@a) dw =1. 


Each f, is represented by a polynomial in w,/B with coefficients changing 
at each integral value of this variable. Furthermore, each f,,, cuts off 
at mp Berond which it is identically zero. The first three even spectra 
are as follows : 


ussmaale-3) +09) r4(-3)] oem 
aleaC-ay] te 


The higher order normal spectra tend rapidly to the Gaussian form since 


i Ih ~ exp (—pz*/6) for p large, 


1s $(Q)eo(-9) 


APPENDIX IV. 


whence 


RECTIFICATION EFFICIENCY OF A DiopE DETECTOR. 


For small applied voltages (i. e. such that the anode is not more positive 
than +0-2 V. with respect to the cathode) the current vs voltage 
characteristic of typical diodes is of the exponential form opr RE as 
to the retarding-field region with negligible space charge : 


i=c exp (KV), 

where K =(e/kT.) is usually about 10 volt, 
«electronic charge, 
k=Boltzmann’s constant, 


T.=cathode temperature. 
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If now a sinusoidal voltage A sin wt is applied to the detector and 7 is 
the mean current over one cycle 


V=A sin wt—iR, 


es 2Qa0 sa 
ine [ exp [K(A sin wt—éR)] d(wt)—cexp (—KiR)I,(KA). 
“0 


In the absence of an applied voltage there is an initial diode current i, 
given by 


ets exp(—Kzi,R). 
The rectified voltage is thus given implicitly by 
w=(i—i))R=ipR[exp (—Ku)I(KA)—1]. 
If KA <1, I,(KA)1+4K?A* exp (—Ku)1—Ku, 


a) BKB - 
whence Ur I+KiR) 4 


demonstrating the square-law behaviour of the detector. 
At large amplitude KA>1, 


: ~ exp (KA) 
Eo, KAls 
‘ U 1 w log, 27KA 
and thus 1 F= gles. (14+ 7h) + 


showing that the rectification efficiency u/A approaches unity as KA 
becomes large. When A becomes very large the exponential law no 
longer holds for the diode and the current tends to be more nearly 
proportional to the voltage. It is easy to show that in this condition 
the rectified voltage approaches more and more closely to proportionality 
with the amplitude of the applied wave. 


[ 504 ] 


LIV. A Radio Echo Apparatus for the Delineation of Meteor Radiants. 


By A. AsprnatL, J, A. CLEGG, Ph.D., and G. 8. HawkIns, 
Jodrell Bank Experimental Station, University of Manchester ™*. 


[Received February 13, 1951.] 
[Plates XVIII. & XIX.] 


ABSTRACT. 


This paper describes a radio echo apparatus for the continuous recording 
of meteor activity, and the determination of the coordinates of active 
meteor radiants. It has been in operation since September 1949, and the 
occurrence of well known night time showers during this period have 
provided an opportunity to estimate its performance. The results 
obtained during the Geminid shower of 1949 and 1950 are described as 
an example, and its accuracy and resolving power are discussed. 


§1. INTRODUCTION. 


THE radio echo technique has been used extensively during the past 
five years to detect meteor ionization in the atmosphere, and with the 
development of methods of radiant determination (Hey and Stewart 
1947, Clegg 1948, McKinley and Millman 1949) and velocity measurement 
Hey and Stewart 1947, Ellyett and Davies 1948, Manning, Villard and 
Peterson 1949), has proved valuable in the observation of meteors under 
conditions which render visual watching impossible. The earlier radio 
investigations were concerned primarily with the well-known night time 
meteor showers, but the experiments of Hey and Stewart (1947) and 
Prentice, Lovell and Banwell (1947) revealed a high level of daytime 
activity during the summer months of 1945 and 1946, and more compre- 
hensive surveys in 1947 (Clegg, Hughes and Lovell 1947) and 1948 
(Aspinall, Clegg and Lovell 1949) confirmed the existence of a number of 
rich meteor streams incident on the sunlit side of the earth between May 
and July. This work has now been extended, and an apparatus has been 
developed to record meteor activity continuously throughout the whole 
year. It is designed to determine the radiants of the major daytime and 
night time showers occurring in the Northern hemisphere, to detect the 
minor streams, and to measure the hourly rate of sporadic meteors. 

The method of radiant determination is similar to that described by 
Clegg (1948) ; but in order to obtain a more complete record the single 
rotatable aerial used in earlier experiments (Lovell, Banwell and Clegg 1947) 


* Communicated by Dr. A. C. B. Lovell. 
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has been replaced by two fixed arrays which operate independently, 
producing low elevation beams, one directed towards the North of West 
and the other towards the South of West. The radiant coordinates are 
calculated from the times at which echoes appear at maximum range 
on the two systems. 

The equipment has been fully operational since September 1949, and 
the results of the continuous survey and of the work on the summer 
daytime streams which has been carried out since that -date will be 
discussed in detail elsewhere (Aspinall and Hawkins 1951). The purpose 
of the present paper is to describe the apparatus and to discuss the 
accuracy and limitations of the method of radiant determination employed. 


Fig. 1. 


TRANSMITTER 


A brief account of the observations made during the Geminid meteor | 
shower is included, however, in order to illustrate its performance. 


§2. APPARATUS. 


A schematic diagram of the apparatus is shown in fig. 1. It operates 
on a frequency of 72 Mc./s., and comprises two independent beamed 
aerials A, and A,, which are directed at low elevation, on azimuthal 
‘bearings of 242° and 292° respectively *. These arrays are common to 
the transmitter and receiver, and the pulses from the transmitter are 
radiated simultaneously by both. The received signals are fed through 


a 
* The azimuthal directions being measured in degrees East of North. 
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the transmitter—-receiver switches, S, S,, into the preamplifying stages, 
P,, P., and after further amplification and detection are applied to the 
grids of the intensity modulated display tubes, D, and D,. The sequence 
of operations is initiated from a master control unit by the unit T,, 
which provides pulses to trigger the transmitter, the multivibrator M, 
and the common time base TB of the two display tubes. The square waves 
from the multivibrator are used to suppress the preamplifying stages 
P, and P, alternately, and are also applied to the grids of the cathode 
ray tubes D, and D,, so that the signals from each aerial are displayed 
separately side by side, and are photographed on a film which moves 
continuously in a direction perpendicular to that of the time bases. 

The unit T, is capable of providing pulses with a number of different 
recurrence frequencies, and the apparatus can be triggered in a variety 
of ways. For normal observations the transmitter is usually operated at 
a pulse recurrence frequency of 150c./s. and the time base unit at a 
frequency of 75c./s. Under these circumstances the transmitter pulses 
appear twice on each display, at the beginning of each time base and at a 
range of 1000 Kms. This serves as a useful check on the range calibration, 
although subsidiary range markers can be supplied from the triggering 
system at intervals of 200 or 250 Kms. along the time base. One of the 
most serious difficulties encountered during the early experiments was 
that of differentiating between short duration echoes and the random 
noise impulses which present a similar appearance on the photographic — 
record. This form of interference was considerably reduced by inserting 
a discriminator unit, similar to that used by Davies and Ellyett (1949), 
in the output of the receiver, but it has also been found advisable to trigger 
the transmitter with a pair of pulses separated by approximately 300 psec., 
so that a true echo appears asa double. A portion of a record, taken with a 
film speed of 12 cms. per hour, is shown in fig. 2 (Pl. XVIII). The continuous 
lines A, B, A’, B’ correspond to the double transmitter pulses, and a number 
of echoes are visible as at C, D and E. 

The photograph in fig. 3 (Pl. XTX.) shows the general arrangement of the 
two aerial systems, which are situated at longitude 2° 18’ W., latitude 
53° 14’ N., and are placed symmetrically on either side of the building 
housing the transmitter and receiver. Each array consists of 6 Yagi 
aerials mounted at horizontal distances of 1-2 apart and at a height of 
1-57X above the ground. They produce identical beams of elevation 
8-5° and of half amplitude width -+-5° and have a power gain of 165 over 
a half wave dipole. The horizontal and vertical radiation patterns are 
shown in figs. 4a and 4b. 

The minimum detectable signal at the receiver is 7< 10-1 w., and in 
normal operation the peak power of the transmitter is 5 Kw., with a 
pulse length of 8us. Simultaneous visual and radio observations have 
shown that under these conditions the echo rate produced by an active 
shower, when the radiant is at 90° elongation to either of the aerial beams, 
corresponds closely to the rate of occurrence of visual meteors. 
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§3, Tae DETERMINATION oF METEOR RADIANTS. 


The method of radiant determination has been described previously 
(Clegg 1948) and only a brief description of its application to the present 
apparatus will be given. 


Fig. 4a. 
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It follows from the condition of specular reflection (Hey and Stewart 
1947, Lovell, Banwell and Clegg 1947), that when the earth passes through 
an active meteor stream, the radio echoes observed on frequencies of the 
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order of 72 Mc./s. are returned from a narrow zone of the atmosphere, of 
limited depth, which lies in a plane passing through the observing station 
and perpendicular to the direction of the radiant. The rate of occurrence 
and ranges of the echoes depend on the position and orientation of this. 
plane relative ‘to the aerial coverage, and vary with time as the radiant 
moves across the sky. It can be shown that with an aerial directed at 
low elevation on an azimuthal bearing of 270°, echoes attain a 
maximum range when the radiant is close to the central meridian. For a 
radiant close to upper culmination, the passage of the echo zone through 
the beam is accompanied by the initial appearance of echoes at short 
ranges, followed by a gradual increase in range to a maximum value 


Ryex, 2nd a final sudden decay in the echo rate. 


Declination of radiant. 


0 50 100 150 200 250 6-6 
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If the aerial characteristics are such that echoes can be detected down to 
the horizon, R,,,,, is equal to the distance at which the meteor zone cuts the 
horizon circle of the station, and for a typical shower of medium velocity, 
for which the mean height of the zone can be taken as approximately 
95 Kms., R,,,,=1100 Kms. Under these circumstances the time of 
maximum range, T, is the time at which the azimuth of the radiant differs 
from that of the beam by 90° and corresponds to the transit of the radiant. 
For aerials directed along azimuths 270°+-6, and 270°+6,, the echoes 
reach a maximum range at times T+¢, and T+#,. The time difference 
t,—t, is a function of the declination, 5, of the radiant, 0,, 5, and the 
latitude of the station. This function is plotted in fig. 5 for the case of 
6,= — 28°, 0.= + 22° and lat.=53°N., and from this curve the declination of 
the radiant may be found. From the curve labelled 1100 Kms. in 
fig. 6, t; and t, may be determined separately and the time of transit then 
found from the observed values of T-+-t, and T+-é,. 
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It may be noted that ¢, is negative for radiants which transit south of the 
zenith and echoes appear first in aerial 1 and then in aerial 2. For 
radiants which transit north of the zenith (dec. >53°) this order is reversed, 
and ¢,—t, becomes negative. If the declination is greater than 76° the 
echo zone never enters the coverage of aerial 2, but remains for many hours 
in the beam of aerial 1. A radiant position can still be determined however 
by fitting theoretical range curves to the range/time plots. 

In practice, owing to the conditions governing the reflection of radio 
waves from the surface of the earth, the aerial sensitivity falls off rapidly in 
directions close to the horizontal, and echoes are never observed out to 
ranges exceeding 1000 Kms. Under these circumstances the times of 
occurrence of echoes of maximum range are advanced by an amount 
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depending on the declination of the radiant and the observed value of R,,,,.- 
This has been taken into account in the curves labelled 900 and 800 Kms. 
in fig. 6, from which the time of transit is found, but since the displacement 
is to a close approximation the same for both aerials, the curve of fig. 5 
and determination of declination is not affected. Finally, a small correction 
has to be applied if the mean height of the meteors in a stream is not 95 Kms. 
This correction, which affects only the time of transit, is given in Table I. 


$4, Resutts OpTainED Durine Tur GEemMiInip MerEor SHOWER. 


The range-time plots obtained during the Geminid shower of 1949 and 
1950, which are shown in fig. 7, provide a typical example of the results 
obtained with the apparatus. Times are shown as absicisse, and the 
ranges of individual echoes are indicated by the lengths of the vertical 
lines. ‘The results obtained on each aerial, between 0000 and 0400 U.T. on 
Dec. 13, 1949 and Dec. 13, 1950 are shown separately. 


510 Ona Radio Echo Apparatus for the Delineation of Meteor Radiants 


These plots indicate clearly the general variation in echo rate and range 
which occur as the echo zone sweeps through each beam in turn. The 
occurrence of echoes at maximum range, and the subsequent sharp fall in 
rate, appear at approximately 0110 U.T. on the first aerial, and at approxi- 
mately 0240 U.T. onthe second aerial. The radiant positions obtained from 
these plots are as follows :— \ 


Time of local transit. | Radiant coordinates. 
Dec. 13, 1949 0214 U.T. R.A. 111-5°, Dec. 32°5° U.T. 
Dec. 13, 1950 0218 U.T. R.A. 112-2°, Dec. 32:5° U.T. 


When the radiant position has been determined, theoretical range-time 
curves can be constructed. These are shown as dotted lines in fig. 7. In 
the case of an ideal point radiant, in the absence of sporadic meteor 
activity, these envelopes should contain all the echoes, and they provide a 
useful indication of the diffuseness of the radiant. 
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§5. AccuRACY OF RADIANT DETERMINATION. 


The results obtained with this apparatus during the major daytime and 
night time showers have shown that when the echo rate, during the passage 
of the echo zone through the beam, exceeds 30 per hour, times ¢, and ¢, can 
be measured to within -+-5 minutes. This corresponds to an error of 
+1° in right ascension and +1-5° in declination. In the exceptional case 
of a major stream with a radiant lying within 20° of the celestial pole, the 
position must be determined by the curve fitting method, and the accuracy 
is necessarily reduced. 

For weaker streams the accuracy of radiant determination falls off 
markedly with decreasing echo rate, particularly if the shower rate does 
not exceed that of the sporadic background by a sufficiently high factor. 
Indeed, in the case of minor streams occurring during times of high sporadic 
activity, it is often difficult to distinguish between true and spurious 
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radiants, and for this reason it has been found necessary to classify showers 
for the purposes of radiant determination in the following way : 


Ciass A, which comprises the most active and prominent showers, and 
which is defined by the following criteria : 

(1) During the passage of the echo zones through the two beams there 
must be at least 3 echoes of range greater than 700 Kms. appearing 
on each aerial. 

(2) Two-thirds of the total number of echoes observed must fall within 
the theoretical range envelopes. 

A statistical analysis shows that the probability of a radiant of this class 

being spurious is less than 1/100, so that all showers which fulfil the fore- 
going conditions may be taken to be real. 


Cuass B. This is defined as follows :— 


(1) During the passage of the echo zone through the two beams there 
must be at least 3 echoes of range greater than 600 Kms. appearing 
* on each aerial. 
(2) One-half of the total number of echoes observed must fall within the 
theoretical range envelope. 


A range-time plot fulfilling these criteria is not by itself considered to 
afford definite proof of the existence of a shower, but only to provide 
confirmatory evidence in cases where a Class A radiant appears on the 
previous or the following day. 


Ciass C. This Class is used to monitor the meteor activity from a region 
of the sky containing a suspected radiant for which the plots 
do not fulfil the conditions required by A and B. 
In estimating the radiant position of a shower which is active for several 
days, a mean of individual radiant positions is taken, weighting the plots of 
Class A, B and C in order of importance in the ratio 2:1: 0. 


§6. Resotvina Power. 


It occasionally happens, particularly during the summer daytime 
streams, that two radiants are simultaneously active ; and if the apparatus 
is to resolve between two such centres, their range-time plots on at least 
one of the two aerials must be separated by approximately 30 minutes. 
This time separation depends on the relative positions of the two radiants, 
and is most conveniently expressed in terms of their angularseparation, s, and 
their relative position angle, ¢, as indicated in fig. 8. This figure represents 
a small portion of the celestial sphere, viewed internally, with two radiants 
situated at R, and R,. If ABC is the meridian bisecting the great circle arc 
R, R, at B, then the position angle ¢ is defined as the angle ABRg. 

Fig. 9 shows the angular separation s required for resolution on each of 
the two aerials, plotted as a function of 4 for different values of the 
declination of the radiant R,. It is evident from these curves than any 


ot 
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pair of radiants whose angular separation is greater than 20° can be resolved 
by at least one of the two aerials, while for more favourable cases the 
degree of resolution is considerably higher than this. 


Fig. 8. 
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LV. Transit-time Phenomena in Electron Streams.—-III. 
The Electron-I[on Plasma and Beam Fluctuations. 


By D. K. C. MacDonatp, 
Clarendon Laboratory, Oxford *. 


[Received February 13, 1951.] 


ABSTRACT, 


Continuing the investigation started in earlier papers, we analyse here 
the equilibrium (isothermal) fluctuations (“noise”) in a travelling 
electron beam taking into account full space-charge interaction. The 
results are then considered in relation to a number of practical electron 
valve problems. Because of its relevance to the general problem we also 
present first a short discussion of the characteristics of electron-ion and 
electron plasmas based on a “ lattice ’? model. 


§1. INTRODUCTION. 


In two previous papers in this Journal (MacDonald 1949, 1950) we have 
studied systematically the problem of fluctuations (‘‘noise’’) in a travelling 
electron beam and their growth. In the first paper we neglected 
interaction entirely while in the second an average, or ‘‘ smoothed ’’, 
potential-variation due to space-charge (e. g. Child’s law) was taken into 
account. In the present paper we derive the equilibrium shot noise in a 
travelling beam where the complete fluctuating interaction is included 
insofar as this is expressed by Poisson’s law at any instant. A general 
expression emerges which enables a number of very suggestive comparisons 
to be made with a range of valve problems. 

When undertaking this work we were naturally led to consider models 
for an electron-plasma and in particular found a “ lattice’ model 
illuminating in several respects. We have therefore thought it relevant 
at this time to present first a brief discussion of such a model although it 
should be emphasized that the fluctuations analysis which follows is 
essentially independent of this model, being in fact based on quite general 
premisses. 

§2. Puasma MODELS. 


In the past it has been customary (e.g. Tonks and Langmuir 1929, 
Bohm and Gross 1949, Twiss 1950, Gabor 1950) to consider an electron 
plasma or electron-ion plasma as a quasi-continuum of essentially gas-like 
(cf. also v. Laue 1918) properties. If we assume the average inter-charge 
distance to be d then the potential energy per charge will be in the order of 
e/d, while at temperature T the random kinetic energy~kT. Since 
generally in electron-tube problems kT > e2/d(=(pe®)"/8) the gas (or rather 


* Communicated by the Author. 


516 D. K. C. MacDonald on the 


perhaps fluid) approximation is essentially justified. However, it appears 
useful, as has often been done in the corresponding metallic ‘‘electron-gas” 
problem, to consider the properties—particularly as regards the frequency 
spectrum—of a ‘rigid ” charge lattice of inter-atomic distance d. The 
well known difficulties associated with three-dimensional lattices restrict 
us to a linear lattice, following closely Brillouin (1946). 


2.1. Electron-Ion Lattice. 

Consider then a linear chain of ions (mass M, charge +e) and electrons 
(mass m, charge —e). In order to maintain this chain in a stable state it 
will be necessary to assume some repulsive force in addition to the Coulomb 
attraction *. Born assumed a potential law of the form A/r” when dealing 
with the alkali-halide solid and we shall follow him here. 

Let d be the distance between electrons, and let U be the potential 
energy between an adjacent electron-ion pair. Then it may readily be 
shown (cf. Brillouin, loc. cit., Cap 3) that two separate modes of vibration of 
the lattice can occur. In the lower-frequency branch, generally known as 
the ‘‘ acoustical ” branch, neighbouring electrons and ions will essentially 
be oscillating in phase with one another while, after a ‘‘ forbidden-zone ~ of 
frequencies, we arrive at an upper branch—the ‘ optical ’’ branch—where 
electrons and ions will effectively be in anti-phase. 

Tf, as a first approximation, we consider only “ nearest-neighbour ” 
interaction then the frequency/wavelength relations for the two branches 
are given by 


. gare 1/2 
» (@U I 1 pe Tee ee, | 
—\or), |m* M (=r) as 0 
if _ .md 12 
cls eel ee Lapel see ae | 
“o-\ ae), (mt Mm (+m) ~~ inven | Fe? 


where 7 is the distance between an electron-ion pair, w, and w, are the 
(angular) acoustical and optical frequencies and the differentials are to be 
evaluated at r=d/2. 

Since m/M <1, these reduce to 


ow [2 (22),] ant 

is a (ae), | sim To 2 ape (3) 
2 /aU 

and Wo= Js (s), | F . « . “ 4 ey eee 4 (4) 


Now with the foregoing assumptions we may set 
(0?U/Or?)>=8(n—1)e?/d?=8(n—1 ep, 


_——S$.S 

* It should be emphasized that this is merely a convenient ad hoc assumption 
to enable us to treat the problem effectively as a lattice. In fact, the thermal 
energy of the plasma could be regarded for our purpose as playing in part the 
réle of this repulsive potential. If T is too small, then electron-ion recombina- 
tion will simply occur. 
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where p is the density of electrons (or ions) in the corresponding three- 
dimensional “ lattice’. For A>d, (3) may then be written 


16(n—1 ad ' 
a=, /[-ae Sse Rak ear od ge (3a) 


corresponding to travelling waves with a phase—(and group)—) velocity 


d 16yep 
Satr! wil M ) 6 re aks Par eT s, iol Meni d's (5) 
where y=n—1. 


For A~d, dispersion becomes evident and a limiting frequency, @,, 
occurs for A=2d where the group-velocity becomes zero (standing waves) : 


ao 16yep 
a= /(?) . . . ° ° . . . (3b) 


This may be compared with the “ plasma-ion (angular) frequency ” (ef. 
Tonks and Langmuir, loc. cit ) given by 


ap=, | (AR). Bay CAD nite 81 Sc) 


From (4) above it is clear that the “ optical” frequency branch has 
degenerated into a single standing-wave vibration with 


ov | (=X), Met tp sida) 


which then corresponds directly to the “ plasma-electron frequency ”. 
To our approximation the wavelength may then assume any value greater 
than 2d which is the common limit for both classes of oscillation. 


2.2. Electron Lattice. 
With a linear lattice consisting purely of electrons, only the acoustical 
branch can exist and we have now, to the “ nearest-neighbour ” 


approximation, 
bis < 4 (PU _ 7d (6) 
w= E a | sim X Z nh es = (OAs 


Taking now U, the potential energy between two adjacent electrons, as 
simply ¢2/d in this case we see that for A>d we have travelling waves with 
phase —, and group —, velocity 


BR 


2 : 
imv2= = ° . e . ° ° . . . a (7b) 


d 


Again a maximum frequency (standing waves) occurs for 


a ea 


and Nene (i eee ee eee (0D) 


which may again be compared with the familiar plasma-electron frequency. 


(7a) may also be written : 
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2.3. Effect of Thermal Energy. 

So long as kT is less than, or the order of e?/d—the electrostatic energy— 
we shall not expect in our model appreciable deviation from these results 
within their degree of approximation. ihe however, kT exceeds e?/d 
significantly we should expect the velocity in (7) to be given now rather by 


kT 
Impnw 1. € e. ype ~ MS ). . (8) 
kT 
Ma Ne ey: on (=) “/ (s) |: a (2), ae 
4 
If in fact we set a= |(= =), 


the plasma electron frequency, we have 


im 2) () | arf) = = 0 


which is essentially just the Debye characteristic length (cf. e. g. Bohm and 
Gross, loc. cit.) 
3. Noise in Travelling Electron Beams. 

We turn now to a beam travelling as a whole with velocity vy such that 
dmv2=®, say. We wish to determine the equilibrium observable 
current fluctuations in such a beam due to thermal agitation. The 
“atomicity ” of electricity (cf. the “ shot effect ’’) will introduce itself 
naturally as a consequence of the frequency characteristics of the discrete 
lattice discussed above. 

Consider then an electron beam of mean density, p9, velocity v), which 
has been travelling “ for a sufficiently long time ’’, to be regarded as being 
in thermal equilibrium at a temperature*, T. Let us also take a length, L, 
of the beam of “ sufficient ” extent for the end conditions to be neglected 
and let us assume an arbitrary disturbance of electron density, p, and of 
velocity, v, which has been analysed into Fourier components over this 
interval. Let the component at (angular) frequency, w, say, of electron- 


density be 
Pn =P. sin (« (— =~) +6), RPMs RAL 


java sin (w(t =) +4). eee Ae Acai) 
% 


The corresponding component of current is given by 


1, =A(p.%+po®.,) | 


and of velocity be 


(where this is essentially a vector addition) ul) 


=I, tI, say 
where A is cross sectional area of the beam. 


ee SS 
* The question of the actual temperature appropriate to a beam is returned 
to briefly later. 
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From (10a) using Poisson’s equation we derive the corresponding potential 
variation, and hence that the perturbation energy corresponding to (11) is * 
__ 722 AL * pov AL m 


E oe (12) 


Ww 


We can then write down immediately the probability distribution function 
for the beam current components and phases 


S (Ty, 1, ; 8, 6) dl, d1, dd dd 
aLI,? mLI,? 


iby Lm ~ A@ekT ~ deg, AkT : 
Ba (scan) (<u) Tee ae dI,d1,d0dd. (13)+ 


And, in particular, 


=F 2kTA = 4ep,AkT 
Peed hye ceee ie 
[ it alee IF 7.4 mL ; (14) 
Since, further, in a bandwidth B there will be BL./vy) components, we may 
write 
ae BL /2-E2 2 i 
sE| = > (432) € kTA me B, 5) 
Up 2 277V Mv 


which may finally be arranged as 


a PET Oe 
2a —— 
1G = 53 1+( -)  2e1yB, Ae oo w= (1 6) 


noting that I, (the d.c. current)=p)v)A, and that w,, the electron plasma 
frequency, is given by w>—4zep9/m. 
3.1. Discussion of Derived Formula t. 

In an actual thermionic valve one may expect that the above expression 
for 81? should lead to a lower limit for the average value of current noise 
which will be observed. Let us consider first a parallel-plane diode at low 
frequencies (thus w<w,) ; in the strongly space-charge limited condition 
let us assume Child’s Law to hold as a sufficient approximation 


iPod Oaks 
1 2 tt dl, 
therefore v, = 5 guile ( Ia= ww) : 


* The omission in the first integration of a term linear in 2 corresponds of 
course to the “ removal” of a d.c. potential variation which is irrelevant. 

+ We are using here the well-known result that for a system in thermal 
equilibrium characterized by canonical coordinates p;, 4%, say, the distribution 
function for these coordinates is given by 

f (Dis %) Woe dqy ~ €~ PP WET dp, dqy. 

+ This comparison of our derived formula with particular valve problems is 
to be considered as essentially heuristic in nature rather than in any way 
rigorous; we mention below some aspects which require further careful 
consideration. 
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Now eV, represents the final anode energy and it might seem therefore 
rather more appropriate to insert in (16) some effective average value for 
®. A direct time-average, for example, for ® yields in fact D.g=1@ 
whence, substituting, we find 


=; 5 
él7 =4kTq, (5) B. e e e e e e e (1 7) 


This may be compared with the classical detailed analysis of North et al. 
(1940) and Schottky and Spenke (1937) of the noise in a diode at low 
frequencies which under extreme space-charge limitation yields 


5H} —aKTy, [ 3 (1-3) |B. aioe eee 


We might therefore suggest that the diode wnder strong space-charge 
limitation behaves very much as a current stream in thermal equilibrium. 
That is, we may consider that the available accelerating energy is used 
efficiently in a quasi-reversible manner through the medium of the space- 
charge barrier to maintain the noise at practically its lower thermal level. 
Again, the electron-current in a diode under purely retarding-fleld 
conditions has its z-component of velocity given everywhere by 4mv2=kT, 
and it seems not unreasonable in this case to define ©, the effective 
transport energy, by }mv2—kT/2, =kT/2, whence from (16) we have 


CL FN: a er To 7 G2) 


in agreement with other approaches to the problem and with experiment 
(MacDonald and Firth 1947). ; 

On the other hand, in the case of a saturated diode (i. e. having a high 
enough anode voltage to prevent the formation of any potential minimum) 
where experiment and theory show that 6I?7=2eI)B it is clear from 
comparison with (16) that the system is then very far from thermal 
equilibrium, as is otherwise quite evident. 

The discussion above of the space-charge limited diode suggests that a 
cylindrical, or, a fortiori, a spherical-structure would exhibit lower noise 
since in that case the final energy is effectively achieved in transit more 
quickly than in the parallel-plane valve. A detailed analysis by Wainstein 
(1947) of the cylindrical diode based on the concept of a single-valued 
velocity-equivalent of the random thermal energy as first introduced by 
Rack (1938) proposed this and his analysis indicates that for r,/r, 20, the 
reduction factor over the parallel-plane diode would be as low as 
~(0-4)?=0-16. On the other hand, a direct application of our equations 
(16) and (17) suggests that one might expect a lower limit to be given by 


= 1 
5 akg, (5) B. Sa eee eee oe 


Turning to the general travelling beam problem without restriction on 
frequency, it is clear from the first part of this paper that (16) will only be 
valid up to an (externally ‘‘observed’’) frequency f=V9/Ayin» Where Amin iS 
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defined by (96). v9 is given by 2mvy—=® or—more strictly it appears— 
for a beam originally emitted from a cathode at temp. T, by 
4mve= ®—kT/2. 


~ 2 2 
Thus tbh aan Agila 
one AMG 
a2 
therefore 1+ —= 29, 
W 


and thus at this frequency limit we have 
s2—= 2e1,B, ° ° e e ° . e e (21 ) 


characteristic of a purely random electron flow (pure “ shot effect ’’). 
Thus, generally, 


es 2 
SH 55 (1+ 5) 2108 HiMeuiltiarhs's §(LGa) 
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The introduction of the limiting frequency and the approximation of 
spectral variation in (16a, b) are closely analogous to the limiting frequency 
involved in the Debye treatment of the specific heat of a crystalline solid. 

It is very interesting to notice here also that Pierce, Smullin (1949) 
and Robinson (1950) have recently derived an expression for the noise 
current at high frequencies emerging from an electron gun, again on the 
basis of Rack’s approximation in combination with a detailed analysis of 
the gun as a diode at high frequencies (cf. Llewellyn and Peterson : 1944). 


They find 
—_ ON kL, 
ol? = (4—77) a) (53) 2e1,B, Ot. | ee gece (22) 


which may again be compared with our (16). 

The analysis and experiment (Cutler and Quate 1949) also indicate that 
if the beam then drifts freely the noise magnitude will show an oscillatory 
variation along the drift-space about this value. Ultimately, however, this 
variation must decay with transit time, 7, due to the thermal velocity 
spread ; the decay time will be given approximately by t~1/27f (kT/2®) 
(cf. MacDonald 1949). . 


3.2: Outstanding Considerations. 

A more precise application of our results to valve structures would 
require a rather careful consideration of the influence of the boundary 
conditions at the cathode in terms of the effective degrees of freedom of the 
system. The situation then is of course very similar to that encountered 
in a discussion of the thermal energy and specific heat of a crystal, 
2N 
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but in the latter cases we are not indeed interested in the actual behaviour 
in the neighbourhood of the physical boundaries but only in so far as it 
may affect our specification of the normal modes of vibration. Conse- 
quently, an idealized artifice such as Born’s “Cyclic crystal” is adequate 
in that situation. It appears in fact to the writer that if one specifies 
the current emerging from a thermionic cathode by I—Apv then only v 
may be regarded as an independent degree of freedom, “ near ” the cathode. 

Since also electron beams in practice presumably correspond more nearly 
to adiabatic than isothermal systems care would also be necessary in 
ascribing an equilibrium temperature in a particular situation and more 
work remains to be done on both these aspects of the problem. 
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SUMMARY. 


A method of estimating the cross-sections for meson production in 
proton—proton collisions near the threshold is outlined, together with 
detailed calculations for spin zero meson fields under certain simplifying 
assumptions concerning the nuclear forces. It is found that, for the final 
continuous neutron—proton system, a low energy for the neutron—proton 
relative motion is favoured; this, together with a large contribution from 
transitions to a final bound deuteron, leads to a meson spectrum well 
peaked at the highest allowed energies in agreement with recent experi- 
ments. While the total cross-section is dependent more on the shape and 
size of the inter-nucleon potential than on the meson type, the spectrum at 
a given angle and, more particularly, the angular distribution are critically 
dependent on the parity of the produced meson. Comparison with 
experiment favours scalar mesons. 


§1. INTRODUCTION. 


THE problem of the production of mesons in nucleon-nucleon collisions has 
previously received treatment by Heitler (1943-5) and his collaborators, 
and more recently by Morette (1949) and by Foldy and Marshak (1949). 
The emphasis of earlier work is on the effect of strong radiation damping of 
the meson field in reducing the otherwise divergent cross-sections calculated 
by perturbation theory. In the present paper we do not intend to deal 
with the case where the nucleon kinetic energies become greater than their 
rest mass, so that some approximate non-relativistic treatment should be 
applicable to the nucleon motion. The mesons, however, are given 
relativistic treatment. Recent measurements by Cartwright et al. (1950) 
and Peterson et al. (1950) at 345 MeV. are such as to come within our 
treatment. It is natural to make the further step of replacing the inter- 
- nucleon interaction by a phenomenological one. In doing so we ignore the 
fact that nuclear forces are, at least in part, due to m-meson exchange. 
This approach has already been adopted by Foldy and Marshak, who 
discuss briefly the errors likely to be introduced as compared with a full 
field theoretic treatment, such as that of Morette. (These questions have 
been more fully investigated by one of us (E.A.P.) and will be reported 
elsewhere.) The advantage of the phenomenological approach is that the 
internucleon wave functions can then be described with an accuracy 
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limited only by the inadequacy of existing data on scattering to provide 
nuclear potentials. These scattering experiments have been very much 
studied recently, and the analysis of the experimental results by various 
authors, see for example Jackson and Blatt (1950), has led to several at 
least qualitative conclusions regarding the internucleon potential. 

In the present paper we take as a working assumption that :— 


1. The interaction in states of odd parity is so small that it may be 
neglected. This assumption has been introduced by Serber and is based on 
the approximate symmetry about 90° of the angular distributions in 
n—p scattering. 

2. The interaction in states of even parity is charge independent and 
may be best represented by a singular form of potential with a long tail. 

No account is taken of tensor forces, or of possible spin orbit forces. 
Some such addition is clearly necessary to account for the observed p—p 
scattering and will be considered in later calculations. We further 
confine attention to the production of positive mesons in proton—proton 
collisions, limiting ourselves to scalar and pseudo-scalar meson fields. 


§2. GENERAL FORMALISM. 


The Hamiltonian of the interaction between the nucleon and meson 
fields will generally be represented by 


H’= ft(x)Qd(x)b(a)dz+adjoint, . . . . . (1) 


where 7 is the field operator of the nucleons and ¢ the field operator of 
the mesons. The operator Q differs according to charge and spin 
dependence of the meson theory. The matrix-elements of the Hamiltonian 
(1) may be evaluated following a method discussed by Becker and Leibfried 
(1946), the term appropriate to a transition from an initial state o to a final 
state feach with A nucleons being 


Hi fay... fdeyphley-.-ag){ 2 (Qebled+O1stea)h vyler «- -e9) 
(2) 


Here x and ys, are the properly antisymmetrized and normalized wave- 
functions of the initial and final states of the nucleon system, expressed as 
functions of the space, spin and charge variables. The term with ¢f 
corresponds to the annihilation, the term with ¢ to the creation of a a+ 
meson. 

Calculations have been carried out in detail for scalar and pseudo-scalar 
meson fields, using a scalar coupling in the first case and a pseudo-vector 
coupling in the second. The equivalence theorem shows that, in the 
approximation used here, pseudo-scalar and pseudo-vector coupling should 
give identical results in the pseudo-scalar meson theory. With the usual 
Fourier expansion of the meson field we then find for the matrix-element 
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describing the emission of a positive 7-meson with momentum & and 
energy w from a two nucleon system : 


2 
of=J9 Jz f fdr dr bi(r;,7,;,8; 1 2X [TB exp (tk . ro} by (04,7 48;), 
i=1 
e e e (3 S) 


for scalar mesons and 


t=1 
tee eet (SEP) 


Here and in the following S denotes scalar mesons, P pseudo-scalar ones ; 
B®, o®, p\?, are the usual Dirac matrices applicable to the ith nucleon, 
IT is the charge operator changing nucleon i from a neutron into a proton, 
7; and s; are the isotopic and spin variables of the ith nucleon, yu is the rest 
energy of the meson and we work in natural units with A=c=1. 

The frame of reference will be chosen so that in the initial state the total 
momentum is zero, 72. ¢. if p,, p, denote the momenta of the two nucleons 


(P1)o= —(P2)o= Po- 
In the final state it follows from the conservation of energy and momentum. 
that 

(Pit P2)s=P=—k. 


The motion of the nucleon mass centre corresponding to P is fairly slow 
(7. e. non relativistic), so that we may suppose an approximate separation 


&(01,0%2;3 8;,7,)=exp (1P.R)u,(r,s;,7,, «© - « (4): 


where R is the position of the mass-centre and r the relative position. 
vector of the two nucleons. The isotopic factor in uw, and u, may also be 
separated out. In the cae considered the initial state contains two- 
protons so that %,—4,(r,s;) (7), ; the final state may be either the charge 

singlet 1(7), or charge triplet Sr ) state. Using the properties 


1 


Te 1(7)9= Ade? 1(7)o= oy 2 3(7)1; 
Fr) H%r)y—= 5%) 


of the change of charge operator J7, one may write for the metrix elements 
(3) 
Hij=9-V(7/w) fdru}(r, s;){-F B® exp (dik . r) +8 exp (— pk . 4) }uj(r, 8;) 
doeetee 11020) 
Hig=tfleV (7/@) J druj{F (6® . k—wp{?) exp (30k . r) 
+(6®.k—wp) exp (—dik.r)}uy, .« . (5 P) 
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In these expressions the alternative — and + signs hold for the final 
charge singlet and charge triplet states respectively. In accordance with 
approximate non-relativistic treatment of the nuclear motion the wave- 
function wu) and uw, may be reduced to their large components. To the 
first order in the velocity this corresponds to replacing 8 and o (which do 
not mix large and small components) by 1 and the 2x2 o-matrices 
respectively. On the other hand 


<_ —> 
uj pup t/2Mvjio .(V—V)v,, M-=nucleon mass, 
where v, and v, denote the large components of uw) and u; and V is the 
gradient with respect to the relative coordinate r. Introducing these 
expressions into (5S) and (5 P) we arrive at the final approximate form for 
the matrix elements 


Hij=9V/(7/w) fdr v}{F exp (4k . r)+ exp (—3ik .r)}v,, . (68) 
Hi p=tflev/ (m/w) fdr vf {Fo . [kk — (éco/M) 7] exp (37k . r) 


+o) . [k-+ (tw/M)V] exp (—$ck.r)}o, . . . (6P) 


for scalar and pseudo-scalar meson theories. The — and + signs hold 
for the final charge singlet and triplet states respectively. 

The cross-sections for the various processes by which a pair of protons 
may produce a positive meson must now be evaluated by using equations 
(6). The chief problem that remains is the selection of suitable wave- 
functions to describe the motion of the nucleons. In this we shall be 
guided by the information acquired in the fitting of the deuteron and two 
body scattering data by internuclear potentials. We shall find that the 
cross-sections are in fact rather sensitive to the choice of potential. 


§3. TRANSITIONS TO THE DISCRETE STATE. 


We first consider the production of a meson when the transition of the 
nucleons leads to the ground state of the deuteron. Separate consideration 
must be given to the two possibilities that the protons are initially in the 
spin triplet or spin singlet states. In the case of an intial spin triplet the 
wave function vy will be of the form 

%="(0) molt) 


so that the matrix element in the pseudo-scalar case becomes 


Hoyp=tflev/ (7/w) dr 3(o)m Vy{—k . [6 exp (sik . r)—o® exp (—iik . r)] 


+iw/M[o™ . V exp (jik . r)+ 6°. V exp (—dik . r)]}5(0),, Uf 


The evaluation of this matrix element will in general require the 
consideration of integrals such as 


foo exp (dik .r)v,dr, f Vos exp (dik . r)v; dr, mtd 
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where ¥ and v, are solutions of Schrédinger-equations, say, 
VU + (po +MU,)%=0, 
V2v,+(p?-+MU,)0,=0, 


The potential functions U, and U, will in general be different on account of 
the spin and charge dependence of nuclear forces. We shall, however, 
always restrict ourselves to central forces of the same range and shape, so 
that 


| Mnclon mass 7 oe aed 2) 


U,p=Jow(r) and U,=J,wi(r), 


where the J’s are constants of the dimension of an energy and w(r) is a 
dimensionless function of r depending only on one parameter — the range 
of the forces. 

It can immediately be shown from equation (9) that 


(vi—pj) Svyvydr=M J(U,—Up)ojvpdr . . . . (10) 


and this identity allows one to restrict the space integration to an interval 
of the order of the range of the nuclear forces. There is, however, no 
transformation directly available for the integrals (8). It is here that the 
assumption of ‘no interaction” in odd parity states leads to some 
simplification, for according to this assumption, Uy=0 in the initial spin 
triplet states, and we may write 


U=[eXP (ipo . T)—eXP (—%pp - r)]//2. 
The matrix element (7 P) can now be written 
Hop= —tfley/ (7/e) Pom, | CP+ 5 |Fom,,) » [{1(Po+ 2k) —1(Po— 3k) 5 
+ (/M)potl(Pot+ 2k) +1(po—3k) 5], - + + (ITP) 
where as an abbreviation we have put | 


I(po)= Jexp (tpo . r)v, dr. OE IT 4) 


If the initial state is a singlet state its parity will be even, so that 
U, 40 and no simple approximation will be possible to the integrals (8). 
However, as an estimate of the order of magnitude of the even parity 
contribution we shall make use of the approximation 


ia [(2) 2522 Cen | o=—0 [%o)ny) Poles + - (18) 


eg 


which is equivalent to assuming that of all states with even parity only the 
2 proton S-state contributes, and that the momentum of the meson in the 
integrals (8) can be neglected. The expression (13) can then be expected to 
give the right order of magnitude and to represent an upper limit in the 
sense that consideration of the meson momentum will tend to decrease the 


value of the integrals (8). 
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For the calculation ofthe cross-section we need 2'| Ho, |? summed over 
the final spin states of the deuteron and averaged over the orientations of 
the original spin. This is now given by 


(63) 2 
E]HG P= | {itp +34)—Uep— 1)-+ © py(l(po+4k)+1(p,—B1)} 


AJ\2 
+2141(p,)* (F) i 


where 4J=J,—Jo, and we have written J for J,. 
The differential cross-section can be determined from 
do 
dQea ae 
where V is the initial relative velocity of the nucleons and dQ is the 
element of solid angle for the meson. 

For the evaluation of I some assumption has to be made about the inter- 
nucleon potential. Present scattering evidence on the whole favours a 
fairly singular, long tailed potential for which we may use as a simple 
analytical expression the potential suggested by Hulthén (1942) : 


U(r)=J exp (—«r)/[1—exp (—er)].  . . . (18) 


The normalized wave function of the ground state of the deuteron is then 


(eb b2— 
r= Ned e: exp [—Hb—ler][] —exp (—«r)] 


with b= MI 


The depth and range of the potential well must be adjusted so as to fit the 
deuteron data. We shall use the values 


I/ke=1.17X10-¥cm., &J=—46.6 MeV., “J—27.2 MeV., 


where “J and (J refer to the spin triplet and singlet states respectively. 
The integral I from equation (12) may now be easily evaluated yielding 


2arb3(b2-+-1)x5 
1 pe : oa 
(Di MW p\8-+-(S=) 
In the energy region which we can compare with experiment we find 
Po ~~ 400 MeV. and k= 90 MeV., and, therefore, in good approximation 
T(po+3k)+1(po— 3k) ~21 (pp), 
I(po+3k)—I(pyp—2k) ~ —(4k/pp) cos 0 I( pp), 
where @ is the angle between k and py. This gives finally an approximate 


expression for the cross-section describing the formation of a deuteron 
under meson emission 


do f® — B3(b2—1) ke F/AT\2 wPy\? 8c 
eee peeee wampap a (7) +2(ae) ares}. «aD 
4p 


P23 | Hi Pavers es 


> 
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where f? denotes the square of the coupling constant in ‘‘ordinary” units. 
A numerical discussion of this result will be given in §5. 
The corresponding result for scalar meson theory is 
do 12g% 63(b?—1 k35 
a [tna a 8 Ae 6 (178) 
(cee 
4p} 

In this case no nuclear spin changes are allowed owing to the assumed 
central character of the forces. The only effective initial states, therefore, 
are triplet states. For the transition to be allowed the mesons must be 
emitted in states of odd parity. The approximate formula (178) represents 
the emission of a p-wave meson, though some contributions from 
interference with higher angular momentum waves are included. The cross- 
section for scalar mesons increases more slowly just above the threshold 
than that for pseudo-scalar mesons which may be emitted in s-waves. 
However, with 345 MeV. protons (in the laboratory system) the two cross- 
sections have become comparable. If g?=/? the pseudo-scalar total cross- 
section is about twice the scalar at this energy. In the forward direction, 
however, the differential scalar cross-section is larger than the pseudo- 
scalar one. The angular variation cos?@ of the cross-section for scalar 
mesons is more in accord with present experimental evidence than that for 
pseudo-scalar mesons, where interference between s- and d-waves leads to 
an approximate cancellation in the forward direction. 


§4, TRANSITIONS TO CONTINUOUS STATES. 


The transitions to the continuous final states of the proton—neutronsystem 
can be treated in a manner very similar to that applied in the discrete case. 
In the energy region considered the relative motion of neutron and proton 
in the final state is very slow, so that it will be sufficient to deal only with 
final S-states of the neutron—proton system. For pseudo-scalar mesons the 
final S-states of the n—p system can be reached either from the odd initial 
triplet states under emission of an even parity meson, or from the even 
initial singlet states with emission of an odd parity meson. The !S final 
states can only be reached from initial triplet states, the transition 18S 
being totally forbidden. The contribution to the total cross-section from 
transitions to the final nucleon 1S-states will be found to be negligible. 
For scalar mesons the only allowed transitions are triplet—triplet and 
singlet-singlet transitions, the latter giving a negligible contribution. 

Turning to the detailed treatment of the case of pseudo-scalar mesons 
we deal first with transitions to the final triplet states. The matrix elements 
may be calculated as in the previous section. The integrals can be taken 
over, provided that v, is normalized so that 


vp: (Iipry sin praeo) "i a2 1 eh: (18) 


asymptotically, where p denotes the momentum of relative motion of the 
nucleons in the final state. In the evaluation of I a transformation of the 
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type (10) is always made so that only the behaviour of v, near the origin is 
important. In this region v, can be adequately represented by the form of 
the ground-state solution 


v,=(N/r) exp [—3(6—1)xr] [1—exp (—xr)]. 


The normalization factor N has to be derived from the asymptotic 
behaviour of the true continuous s-waves in a Hulthen potential. It is 
then found that 
ated br sinh 2770 iF hg (19) 
x |_a {cosh 27%—cos 27 /(b—a*) } 
where «=p/« and b has been defined in equation (16). In the region of 
interest we may simply put 


N=(bz/xp)?, LE ULC Se hoe 
and it is then found that 
4 oh 
I(q)= — Bq)? et Soa ae) 
(q) al 2 (q) 
where 
a! Po 901 
B(q)= 2 ee Ed ae 


: (b+1)?k? : 
(po—4") (vi+ ae 
The differential cross-section for a transition to a final 3S state with the 
production of a meson in the energy interval w to w+dw is then given by 


do f?b? kMk3 
ABET RTS | Pp! 
AJ\2 


+ $y Pol (pot 1E)+F(p.— 2) +2 (>) ipo |. (1B) 


Using the same approximation as in the treatment of the discrete state one 
obtains 


d 4f? bs kn k 2k?M /AJ\2 
nicer aes | ee ae de oben nee eo 
dwdQ V {is (6-+1)« V p28 |M pe prow \ J 
4 p5 
The result exhibits the same tendency towards cancellation in the forward 
direction as was found for the cross-section in the discrete case. 
Transitions to the final 'S states can be similarly treated. The resulting 
cross-section, however, as is evident from equation (22 P), is proportional 
to the cube of the well depth of the final state. This reduces the contri- 
butions from singlet transitions by a factor of at least 8 and as a result they 
may be neglected in calculations of the present approximate character. 
In the case of scalar mesons the same method leads to a cross-section . 
da 3g? bFkM. 


dwdQ  2V pe {F(po+$k)— F(py— 3k) }, otis hela 
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which can be approximated by 


do _yg@ 12M 
Ado = Viteme eye pw ce” ae is (22 S) 
ey 
‘Po 


§5. NUMERICAL RESULTS. 


We have first calculated the variation of the total cross-section for the 
production of scalar and pseudo-scalar mesons in a transition to the ground- 
state of the deuteron as a function of the energy of the incident proton in 


Fig. 1. 


-29 
off?or g* x0 z ay 


10-0} 


‘Sao 


300 400 Ep Mev. 


The cross-sections for meson production, scalar and pseudo-scalar, in proton— 
proton collisions against energies of incident proton, in cases where final 
nucleons form a bound system. 


the laboratory system. The results are shown in fig. 1. Corresponding to 
the possibility of production of the mesons in an s-state the cross-section 
for pseudo-scalar mesons is considerably higher near the threshold 
(E,=290 MeV.). A maximum is reached at E,=—370MeV. At 
E. =350 MeV., o,, 592 X 10-*°cm.?, op, ~ 12f? x 10-29em.2. It has already 
been noted that the results are rather dependent on the internucleon 
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well shape. For comparison these cross-sections have also been evaluated 
for a square well potential 


@)J—41 MeV.; a=1.85x10- cm. 


and yield results smaller by a factor 1/25. 

Fig. 2 represents the energy dependence of the integrated total cross- 
section for transitions to the continuous part of the deuteron spectrum. 
The result in each case naturally increases more slowly above the threshold 
than in the discrete case. At E,—350 MeV. approx.—the energy of the 
Berkeley experiments—<,,,, is about half of o4;,. for scalar mesons and 
about 2 for pseudo-scalar mesons. The fall in the continuous cross-section 


Fig. 2. 
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The integrated cross-sections for meson production in proton—proton collisions 


against energies of incident proton, in cases where final neutron—proton 
system is in a continuous state. 


indicated in fig. 1 should not be taken too seriously. At these energies 
higher angular momentum contributions as well as higher order corrections 
to the meson field must be considered. The ratio of the continuous to the 
discrete cross-section is consistent with experiment for both types of 
meson. 

The transformation to the laboratory system has the effect of throwing 
more mesons into the forward direction. We shall now denote by k’, w’, 
6’, the momentum, energy and the angle at which the meson is emitted in 
the centre of mass system and by k, w and 6, the corresponding quantities in 
the laboratory system. Taking, for example, the differential cross-section 
for the production of scalar discrete mesons to be 


do/d2'=A cos? 6’ 


es 


Production of Mesons in Proton—Proton Collisions 533 


in the centre of mass system, the corresponding cross-section in the 
laboratory system is 


_,.. # (k cos 8—vw)? 
ee ys (k—vw cos 6) ’ 


where v is the relative velocity of the two frames of reference and 
y=(1—v?)-*. The results have been evaluated for incident 350 MeV. 
protons for the production of both types of mesons. They are shown in 
fig.3. The cross-section for scalar mesons is strongly peaked in the forward 


Fig. 3. 
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Variation of the differential cross-section for meson production with the angle 
of emission of the meson, in laboratory system, is shown for incident 
350 MeV. protons and mesons of maximum energy. 


direction where it obtains its maximum value of 499? x ieee om.*/sterad. 
The cancellation in the forward direction for pseudo-scalar mesons is still 
apparent in the laboratory system. The maximum of the angular 
distribution of 39 f2 x 10-%° cm.2/sterad. is reached at an angle of between 


30° and 40°. 
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The determination of the continuous contributions in the laboratory 
system is more laborious. In fig. 4 do/dwd22 is plotted against the kinetic 
energy for scalar mesons—assuming a proton energy of 350 MeV. for 
observations at 0° and 30°. The differential cross-sections are both 
peaked towards the upper end favouring high meson energies. Integrating 
over all energies we find at 0°, 


ue. = om.? aoe ee -4g? x 10-%° cm.?/sterad. 
70 829" x 10-* cm.?/sterad. and at 30°, IO 3°49? X ai! 


Fig. 4. 
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The cross-section per sterad. per unit energy is shown against the kinetic energy 
of the emitted scalar meson in the forward direction and at 30° to the 
beam of incident 350 MeV. protons. The differential cross-section in 
angle falls off rapidly from the forward direction as the results show, e.g. 


do) 1 (do) 
dQ) 59° 10 \ dQ) e° 


The cross-section falls away very rapidly from the'.forward‘ direction. 
Pseudo-scalar mesons show a very differently shaped differential cross- 
section, approximately proportional to (w—) * so that there is not the same 
favouring of high meson energies. The approximate integrated contri- 
butions to the differential cross-section are for pseudo-scalar mesons 


d 
at 6=0, 0 =30f2 x 10-89 em.2/sterad, 
> da 
6=30°, —~ =20f? x 10-*0 cm.2/sterad. 


dQ 
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Recent experiments have all been carried out with 345 MeV. protons 
and so the correct dependence of the cross-sections on energy cannot yet be 
discussed. The angular distributions at this energy can, however, be 
compared with experiment as cross-sections have been reported at angles 
0°, 18° and 30° to the beam. Agreement appears possible only with 
a cos?@ type of distribution at these angles, and so favours scalar rather 
than pseudo-scalar type mesons. Also the marked peak in the continuous 
spectrum at high meson energies is predicted by scalar theory. 

The forward cross-sections appear to require rather a large coupling 
constant (g?~3, f%~5). However these cross-sections are rather 
susceptible to fine changes in shape and size of the nuclear well. It is 
perhaps also of interest to note that Brueckner (1950) found f? 0.2 to 
fit the observed meson production by photons (assuming the mesons to be 
pseudo-scalar). 


§6. CONCLUSIONS. 


It has not been the purpose of this paper to give the most complete 
description possible of meson production with the use of all the available 
data on the »—p and p-p interactions. Rather, by means of relatively 
simple analytical approximations for the potentials, and by using 
approximate methods, we have shown how important it is, in the low 
energy region (say up to 500 MeV.), to take accurate account of the nucleon 
wave functions. In this region indeed allowance for the detailed behaviour 
ot the nucleons is likely to be more important than the inclusion of any 
field theoretical refinements. Increasingly at higher energies the neglect of 
higher order reactive terms from the meson field will affect the results and 
simultaneously the concept of internucleon potential will lose its validity. 
At the same time multiple meson production will begin affecting the results. 
We believe, however, that in the low energy region the methods applied in 
this paper should be expected to give at least qualitatively correct results. 

If the use of Serber forces can be regarded as satisfactory then a 
considerable difference has been established between the behaviour of 
pseudo-scalar and scalar mesons. The angular distributions are quite 
different : for scalar mesons the distribution goes with cos?@ in the centre of 
mass system; for pseudo-scalar mesons the angular distribution though 
isotropic near the threshold has come nearer to a sin? law at the maximum 
of the total cross-sections. These results can only be fitted to present 
experimental evidence with the assumption of scalar mesons. It is 
interesting to compare this conclusion with that suggested by other meson 
processes. In making this comparison we disregard the information 
about the meson character which can be obtained from processes in which 
they are only involved as virtual particles—magnetic moments, nuclear 
torces, etc. This leaves the evidence from the production of mesons by 
photons, and from the various capture processes in hydrogen and 
deuterium. The first of these has recently been considered by Brueckner 
(1950), and his conclusions are decidedly against scalar mesons. The 
evidence from the capture processes is perhaps less certain. On the other 


536 On the Production of Mesons in Proton—Proton Collisions 


hand, we have made no exhaustive effort to discover how, by alteration of 
the internucleon potential—say by addition of spin orbit forces—the 
results for the pseudo-scalar meson theory in the present investigation 
might be affected. 


One of us (E.A.P.) is indebted to the Department of Scientific and 
Industrial Reserach for a Maintenance Allowance during the tenure of 
which this work was carried out. 
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ABSTRACT. 


A review is made of experimental evidence on the energy losses of 
protons and deuterons in materials of low atomic number, at energies 
below about 350 keV. The discrepancies between results of different 
experiments are discussed, and an attempt is made to arrive at a plausible 
curve of energy loss as a function of energy for deuterons in heavy ice at 
energies below 700 keV. The region 0-100 keV., which is of considerable 
interest in the study of the D-D reaction, is given special attention. 


§1. INTRODUCTION. 

In this paper we shall attempt to deduce, from existing experimental 
evidence, the rate of energy loss of deuterons in a D,O (heavy ice) target, 
for the special case of very low bombarding energies (10 to 100 keV.). 
The determination of this quantity is of obvious importance in obtaining 
D-+D reaction cross-sections from measured values of thick target yields. 

There are not, in fact, any published data on energy losses in ice. The 
nearest approach to what we need is contained in a paper by Crenshaw 
(1942), who studied the energy loss of deuterons in water vapour. 
Unfortunately, however, the lowest energy at which he worked was 
60 keV. ; for information on the rather critical region of lower energies 
we have to look elsewhere. We propose to review the low-energy data 
first, obtaining from them a tentative energy-loss curve which can then 
be compared with Crenshaw’s results over the region common to both. 

Setting aside Crenshaw’s work for the moment, the data most pertinent 
to our problem are contained in several papers by Gerthsen (1930 a & b) 
and his collaborators (Eckardt 1930, and Reusse 1932), concerning the 
energy loss of slow protons in various media (air, hydrogen and celluloid). 
Their results may be applied in the present case, provided that one makes 
the following assumptions :—. 

(a) That the rates of energy loss of a proton and a deuteron having 
the same velocity are identical in any medium. 

(6) That the atomic stopping powers of H and D are identical. 

(c) That the molecular stopping power of D,O can be found by suitably 
combining the atomic stopping powers of two atoms of D and one atom 
of O. 


* Communicated by the Authors. 
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(d) That the stopping power of D,O is independent of its physical state. 

Although there has been no experimental verification of assumption 
(a) for the very lowest energies, Crenshaw (1942) found it to hold good 
between 60 keV. and 300 keV. deuteron energy (=30 keV. and 
150 keV. proton energy) for energy losses in hydrogen. Moreover, it is 
a general principle that a proton and a deuteron having the same velocity 
will lose energy through electronic interactions at equal rates, whatever 
the exact nature of the energy loss process*. We shall therefore take it 
that this assumption is rigorously justified. In general we shall speak 
in terms of proton energies ; any relation arrived at will then be true for 
deuterons of twice the energy. 

With regard to (b), measurements have been made by Crenshaw (1942) 
and others (Joos 1942, Jussuf 1942, and Koops 1938) of the relative 
total amounts of ionization produced by protons and deuterons in 
hydrogen and deuterium gas. The sum total of the evidence is that a 
difference of ionizing effects, if it exists at all, is very slight, and becomes 
noticeable only at very low incident-particle velocities (~6-10’cm./sec.) 
corresponding to 2 keV. protons or 4 keV. deuterons. Thus, for 
proton or deuteron velocities greater than 6-10" cm./sec., one eed) consider 
the atomic stopping powers of H and D to be equal. 

Assumption (c) has purposely been written in a rather vague form. 
At high energies the suitable combination appears to be simple addition. 
The molecular stopping power of a chemical compound is equal (to about 
2 per cent) to the sum of the atomic stopping powers of its constituent 
atoms. This result, however, is based upon measurements with alpha 
particles, and for such energies the difference in stopping power due to 
differences of chemical binding is expected on theoretical grounds to be 
less than 1 per cent. (The problem is more fully discussed in a review 
article by Gray (1944).) The situation is markedly changed at the 
energies in which we are here interested, and according to Crenshaw’s 
measurements the stopping power of D,O is significantly less than one 
would obtain by adding the component atomic stopping powers. We 
shall defer a full discussion of this matter until later in the paper. 

The evidence to support assumption (d) is practically non-existent. 
There are some grounds for believing that the stopping powers of a 
substance in solid and gaseous forms are equal, but this result again is 
for alpha particles, for which no significant difference is expected 
theoretically. (References and discussion may be found in Gray’s 
paper}.) By assuming that the stopping power of heavy ice is equal to 
that of water vapour we therefore subject ourselves to some uncertainty. 

After this general introduction to the problem we will proceed to 
a survey and analysis of the experiments of Gerthsen, Eckardt and 
Reusse. 


ea a en ee en, SIA ee 
* See, for example, Livingston and Bethe (1937), Rev. Mod. Phys. .. 9, 271. 


+ See also Wilkinson, 1948, Proc. Camb, Phil. Soc., 44,114 ; Appleyard, 1949, 
1949, Nature, Lond., 163, 526. 
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§2. THe EneRcy Loss or Stow Protons 1x Matter. 

The experiments on the loss of energy of protons in matter can be 
classed into two types, integral and differential. In the first type of 
experiment the total range in a gas of a proton of given initial energy 
is measured. This is the experiment described by Gerthsen (1930 a & b). 
His results are shown in fig. 1 for the case of protons in air. The ranges 
R are expressed in cms. of air at 1 mm. pressure Hg, and the proton 
energies Kin keV. Gerthsen found (cf. fig. 1) that R as a function of E 
could be very closely represented by the formula 

eet uy wert cess. a ee (1) 
We have applied a least-squares analysis to his data and have found 
that the results are better fitted if the power of E is 0-773 ; this difference, 


Fig. 1. 
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however, is scarcely significant Now the quantity of interest in 
determining the cross-section for a nuclear reaction is —dH/dz, the 
energy loss per unit distance in the target at a given energy If we accept 
Gerthsen’s formula, we have 
Fie earn: De ees, Nees oe we) 
dx 3a 
d#/dx as a function of velocity is plotted in fig. 2. 
In the differential experiments the loss of velocity of protons in 
traversing a very small quantity of matter is directly observed. This 
is the type of experiment described by Eckardt (1930) and Reusse (1932). 
A proton beam was passed through celluloid films of various thicknesses 
Az. For a given film the energy loss 4E in passing through it was 


202 
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measured as a function of the initial energy E. The results are shown 
in fig. 3. By considering the results for the various films one can find 
Fig. 2. 
ee vs. V_ for protons in oir 
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the magnitude of 4E/Az as a function of Ax for a given E. The relation 
between them is found to be linear, of the form j 


AE / 7 
— Fe U8’ Ae. bd el 
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By extrapolating the line to zero film thickness one finds the value 
of dE/dx at the initial energy E. By performing this analysis for a series 
of values of E one can obtain dE/da as a function of velocity. We have 
treated all of the data of Eckardt and Reusse in this way, using least- 
squares solutions throughout, and find that they are consistent with a 
relation 


pis pe eel, mem, <. fy C3 BS) 
This is plotted in fig. 4. 
Fig. 4. 


oon vs. v for protons in celluloid. 
(Eckarct & Reusse) 


6 vy (msec) 


2:0 30 40x 108 


The discrepancy between the results of Gerthsen (equation (2)) and 
of Eckardt and Reusse (equation (4)) is very striking. From the one 
experiment one finds that dE/dx is proportional to the velocity, from the 
other to its square root. It is true that the values of 4E/4x, obtainable 
from the measurements of Eckardt and Reusse, lead more directly to 
a value for dE/dx than do Gerthsen’s measurements of total range, but 
both experiments are open to criticism and merit more detailed 
consideration. 

§3. THe EXPERIMENT OF GERTHSEN. 


A proton beam generated in a canal-ray tube was magnetically analysed 
and passed through a thin celluloid window situated at the centre of 
curvature of a hemispherical ionization chamber. The thickness of the 
window was 80 mu (1 mu=10-*mm.=10 A.U.). The ionization chamber 
could be filled to any desired pressure (up to a few cms. Hg) with air 
or H,. The positive ions produced by the protons in traversing the 
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chamber were collected on an electrode which was made negative with 
respect to the walls of the chamber. The charge collected in a given 
time was registered on an electrometer ; a second electrometer recorded 
the charge carried by the primary beam entering the chamber. The 
ionization current was found to increase steadily with increasing gas 
pressure up to a certain critical pressure p, and thereafter remained 
constant. At the pressure p, the protons are just failing to reach the walls 
of the chamber, so that for this and higher pressures they lose their whole 
energy in the gas. If the radius of the chamber is p the proton range r at 
unit gas pressure is pp,. 

Gerthsen’s experiment consisted in measuring p, for air and H, for 
several values of E between 27 and 57 keV. He found that the range 
in H, at all energies was 2-50 times the corresponding range in air, so that 
the same form of the energy-range relation must hold for protons in air 
and in H,. It is necessary to point out that the value of E was measured 
before the protons passed through the celluloid window. Gerthsen assumed 
that the window was equivalent to the same thickness 6 of air at all 
energies, and calculated 5 by considering the celluloid to be composed 
of hydrogen and “ air-like’’ atoms. In this latter category he placed the 
C, N and O atoms which comprised 94:4 per cent by weight of the 
celluloid. 6 was then added to each observed value of 7 to give the total 
range R. 6 was about 10 per cent to 20 per cent of R, so that any error 
in its value could materially affect the final results. 

There exists an additional way of evaluating dE/dx from Gerthsen’s 
experiment, which he himself describes (1930 a). When the pressure in 
the ionization chamber was very low the ionization current was found to 
be a linear function of pressure. The slope of the current vs. pressure 
curve in this region gives directly the number of ion pairs, per mm. pressure 
of gas in the chamber, produced by the whole proton beam. Since the 
primary proton beam was simultaneously recorded, it is possible to state 
the number, », of ion pairs per cm. of path, in gas at 1 mm. pressure, 
produced by one proton of energy E’, where E’ is the energy of a proton 
after it has lost energy 4E in traversing the celluloid foil. With the 
assumed value of 8, JE and thence E’ could be calculated for various 
values of E. 

To translate n into a rate of energy loss it is necessary to know the 
mean energy, W, required to produce an ion pair in the gas at the 
energy EK’. Gerthsen could not determine W at a single energy, but by 
dividing the saturation current in the ion chamber by the current of 
primary protons, he could determine the total number N of ion pairs 
produced by a single proton of energy E’. The quotient E’/N then gave 
the average value of W over the energy range zero to E’. This ratio 
was found to be independent of the value of E’; it was therefore assumed 
that W itself was independent of energy and equal to the constant 
quotient E’/N. (From our evaluation of the data we find W=38 eV. 
Gerthsen gives W=35 eV., but we have cause to suspect the values of 
E on which this is based—see below.) 
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The rate of energy loss, —dE’/dx, of a proton of energy E’ is then 
given by the product nW (keV. per cm. of path at 1 mm. pressure). 

Support for the assumed constancy of W is to be found in a paper of 
Joos (1942). He quotes an experiment in which the ratio E’/N for 
protons in H, is found constant for energies from 4 to 14 keV. 

As Gerthsen has published his results, the values of dE/dx derived 
by this second method from his one experiment are very different from 
its values as derived from the first method mentioned earlier 
(dE/dx appears to be proportional to v, not to v! as required by equation (2)). 
Upon careful re-evaluation of the results we found that his values 
of 4K, and hence of HE’, were at variance with the integral energy-range 
relation (equation (1)). That is to say, the curve obtained by plotting 
(r+6) against E does not coincide with the plot of r against (E—4E) 
if one uses the values of 6 and 4E postulated by Gerthsen. It would 
appear that the discrepancy may have arisen through numerical errors ; 
in our own treatment of the data we have found that the values for 
dE/dx are essentially the same by either method of evaluation and 
satisfy the relation expressed by equation (2). (See fig. 2.) In 
determining 4E we took the mean of values obtained by three different 
methods :— 

(1) The method used by Gerthsen (1930 a.) (although our results differ 
from his and satisfy the integral energy-range relation). 

(2) We assume that the loss of velocity, dv, in the celluloid foil was 
independent of the energy of incidence E (a result found by Eckardt 
and Reusse, v. inf.) and accept Gerthsen’s value for the air equivalent 
6 of the foil. 

(3) We find directly from the data of Eckardt and Reusse the loss of 
energy suffered by a proton of various energies, EH, in passing through 
a celluloid foil of 80 my thickness. 

We consider that this point is sufficiently important to merit setting 
out the values obtained by the three methods, and this is done in 
Table I. below. One can see that the degree of concordance is quite 
satisfactory and gives some meaning to the average values of JE given 
in the last column : 


TABLE [. 
AE keV. 

0. eR ee ee eee 
E keV. Method 1 Method 2 Method 3 Average 
ee eS 
27-1 6:8 5:9 6-6 6-4 
34-6 2 6-7 or 7-4 
39-6 8-1 7-2 Bd 7-9 
45 8-0 7:8 8-9 8-2 
57 8-3 8-6 10-2 9-0 


pees fee ee ee ee 


The outcome of Gerthsen’s experiment would from this analysis 
appear to be that over the range 20 to 60 keV. the rate of energy loss of 
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protons in air or hydrogen is proportional to the square root of the 
velocity. The actual magnitude of this loss, in keV. per cm. of path in 
either gas at 1 mm. pressure, follows directly from the experimental 
results, and for air is given by 


we (x) =(5:90-+0-10) 10-5 v? keV. per cm., with v in cm./sec. 
alr 
§4. Tue ExpreRIMENTS OF ECKARDT AND REUSSE. 

A proton beam generated in a canal-ray tube was accelerated through 
an additional potential drop and was then magnetically analysed. Protons 
or Hj ions of a specified velocity were thus selected. The range of 
proton energies covered was from 4 to 50 keV. After passing through 
a thin celluloid film the protons were subjected to a further magnetic 
analysis in order to determine their velocity upon emergence. There 
was of course a spread of energy in the emergent beam, but the velocity 
at the peak of the distribution was measured as the significant one. 
The experiment was repeated for celluloid films of various thicknesses 
(from 20 to 330 my). 

It may be noted that the celluloid films were extremely thin, the 
thickest being only about 3000 A.U. Thus it was impossible to measure 
the thickness directly through the use of interference fringes. The 
method employed was (a) to note the change of interference colour 
when one of the thin films was placed over a relatively thick film of 
celluloid, or (b) to superpose portions cut from a single large foil until an 
interference colour was obtained. Unfortunately the account of the 
procedure as given by Eckardt (1930) is not very detailed, and its 
accuracy must remain open to question. 

The conclusions reached by Eckardt and Reusse are (a) that for a given 
incident velocity, the loss of velocity in passing through a foil is 
proportional to its thickness, and (b) that in passing through a given 
foil the energy loss is a linear function of the incident velocity v. These 
results may be expressed by the following equations : 


— Av=ade,) 21 2) 2) DE 
— AE=bv—c. el Te ne may a a 


Tt will be noted that equations (5) are not the same as the equations (3) 
and (4) which we have used in evaluating dE/da as a function of velocity 
The two sets of equations may, however, be readily related. Suppose 
that the rate of energy loss is given by 


1.€. —mvdv=f(v) dx, where m=mass of proton. 
Then ae peer y de oy (i ‘ 2 
a m I. Oe m p'(v) dv, say, oS) 
=-+m[¢(v)—d(v-+ Av)]. j 
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If we accept equation (5a), we have 


1 
Axc=— z Av, and hence the identity : 


1 
— 7 dv=m[ (0) —¢o(v+ Av)] 


=—mdg'(v)4dv+ . 
é il v 
Con: tl =) ae op oe 
sequently CO eas ta at Fo (from (7)) 
dE 
so that — game, STOUl( GO) eee! tee eke, teers) 
Differentiating (8), 
ak dv , 
ee ae OM, FOUL 0) eats emma) 
; dK d?E 
Now we may write AE~ Azx+4 —, Fred (4x), which by (8) and (9) becomes 
en Lae 
Ag ume 24 m Aa 
or _ Fe =0,0—a,Ae, Te ee ee oes AO) 
dE 
where «,v=— ——, and «, and «, are constants. 


dx 


Equation (10) may be seen at once to express the results of both 
equation (3) and equation (5b)—to the former when v is fixed and 
Az varies, to the latter when Ax is fixed and v varies. 

The following table sets out some values of «,v(=—dE/dz), «,v/E' and 
%, as we have calculated them by a least-squares analysis of the 
experimental data of Eckardt and Reusse: 


. 


TABLE IT. 


av keV/mp / 
Rcclidibid oa Og 


0-0395 + 0-0062 0-0164 +0-0026 | 0-00018 + 0-00008 


0-0621 + 0-0040 0-0195 + 0-0013 0-00019 + 0-00003 
0-0890 + 0-:0059 0:0160 + 0-0011 0-00009 + 0-00001 
0-1313 + 0-0061 0-0198 + 0-0009 0-00017 + 0-00001 


The fact that «, and a, are essentially independent of velocity may 
readily been seen. However, the value of (a2v*)/(«.E), which should 
be 4 (cf. equations (10)), is found from Table IT. to be 2-4-0°6 (approx.). 
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In comparing equation (8) with equation (4) we see that they are in 
agreement only if v9 is set equal to zero. The discrepancy probably 
arises because we have made a more careful analysis of the data, using 
least-squares solutions throughout, than did Eckardt and Reusse. In 
this analysis the parameter v, appeared, but it is doubtful that the 
experiments in themselves are accurate enough to give vy any significance. 
In any case it is clear that the energy loss cannot cease at and below v9, as 
equation (4) would demand (v, corresponds to about 0-2 keV. proton 


energy). 


The result of the experiments of Eckardt and Reusse is, therefore, that 


for protons traversing celluloid the rate of energy loss, —dE/dz, is a 
linear function of velocity. The quantitative expression of this relation 
is 


dx 
cm./sec. 

The result is of use to our problem only if we can derive from it the 
absolute value of the rate of energy loss in air or in hydrogen, and thus, 
using our initial assumptions, find the energy loss in a layer of ice. It 
is important to attempt this conversion, since the celluloid measurements 
extend to much lower energies (4 keV.) than do the measurements on 
air and H,. We will therefore consider this matter in the next section. 


~ (=) = (4:45 +0-11) x 10-2(v x 10-8— 0-18) keV. per mp with v in 
celluloid 


§ 5. CONVERSION OF CELLULOID Data To AIR. 


To convert (dE/d2)cetwoia into (dE/dz),;, two alternative methods are 
possible. The first method is briefly as follows: over the small range 
of velocities covered by Gerthsen, the plot of (dE/dz) vs. v (fig. 2) does 
not depart by more than about +10 per cent from the linear relation 
expressed by equation (4). If, therefore, one draws a straight line through 
these points, with an intercept at vy on the v axis, its slope is not likely 
to be in error by more than about +10 per cent. By comparing this slope 
with the slope of the corresponding line for the energy loss in celluloid 
(cf. fig. 4), one finds the number of ems. of air at 1 mm. pressure which 
are equivalent to 1 my of celluloid. The result of this comparison is 


dk 
— (Fe), =(0-40-£0-04)(0% 10-*—0-18), Bs ecard iN 


where v is in cm./sec. and (dE/dx) air is in keV. per em. of air at 
1 mm. pressure. 


The second method is to assume, with Gerthsen, that celluloid Bred: be 
considered as a combination of air and hydrogen. That is 


CAD) dK dK 
hie oe (=). any (i ik 


~a+oam). (2) 


i 
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The constants A and B may be evaluated from the composition and 
density of celluloid. For the density, which is not stated by Eckardt 
and Reusse, we have assumed a value of 1-48 gm./cem.’, which lies 
midway between the accepted limits of 1:35 and 1-60. We thus subject 
ourselves to a possible error of +10 per cent. Our estimate of (dB/dz) 
by this means is 


air 


dk 
= (=). = (0-37 £0-04)(vx 10-8—0-18). . . . (12) 
The striking agreement between (11) and (12) is heartening but probably 
fortuitous. We shall average them and so obtain, as the final outcome 
of the experiments of Eckardt and Reusse, the following relation for 
energy losses of protons in air between 4 and 50 keV. : 


_— (=). = (0-385 +0-04)(v x 10-8—0-18) keV. per cm. at 1 mm, 

We have now arrived at two expressions for the energy loss of protons 
in air, one from Gerthsen’s experiments (formula in §3) and the other 
from Eckardt and Reusse, above. Before converting these into energy 
losses in D,O, we will discuss Crenshaw’s work, on which the conversion 
will largely depend. 


§6. THE EXPERIMENTS OF CRENSHAW. 


We have already made several references to Crenshaw’s paper, which 
describes a very useful study of the energy losses of hydrogen ions in 
air, H, and D,, He and water vapour. 

A resolved beam from a Cockcroft-Walton accelerator was passed 
through a series of small diaphragms into a chamber which could be 
filled to a few mm. pressure with a gas without impairing the vacuum 
in the acceleration tube. After traversing the chamber the beam passed 
once more into high vacuum through a second series of diaphragms. 
It was then bent by a magnet, and fell upon a fluorescent screen. The 
shift of the beam spot caused by the admission of gas to the retarding 
chamber was compensated by altering the magnetic field. Pairs of 
values of magnetic field with and without gas in the chamber were 
obtained for numerous values of initial beam energy, mostly between 
60 and 340 keV. The initial energy was measured with a carefully 
calibrated resistance voltmeter. The energy loss was found from the 
change of magnetic field, this latter being itself calibrated in terms of the 
resistance voltmeter. The experiment was performed for several values 
of pressure in the retarding chamber, and was thus precisely analogous 
to the differential measurements on celluloid. The initial slope of a 
curve of JE vs. Ap led to the value of (dE/dx). End effects in the 
retarding chamber, due to pressure gradients in the diaphragm. systems, 
were eliminated by changing the distance between the entrance and 


exit systems. 
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The point most carefully studied in these experiments was the identity 
of energy losses. by protons and deuterons in H, and D,. Some earlier 
work by Crenshaw and others (1942) had indicated a difference of energy 
losses in H, and Dy, but this was not substantiated. Within the spread 
of the measurements (about +10 per cent at the highest energies and 
+5 per cent at the lowest) all four rates of energy loss were found to be 
equal for a given proton or deuteron velocity. After this had been 
established, further measurements on the energy losses of protons in 


Fig. 5. 


Energy losses in H, (Crenshaw) 


GE 

dx 
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H, and D, were carried out. The results of these two investigations are 
shown graphically in fig. 5. For convenience, \/E is used as abscissa 
and (dH/dx) is given in keV. per cm. at 1 mm. pressure. , 
The other aspect of Crenshaw’s work which interests us here is his 
study of the stopping powers of hydrogen and water vapour relative to 
air, In fig. 6 we have plotted the stopping power of hydrogen relative to 
air over the range 40 keV. to 5 MeV. proton energy. In this graph the 
results of Gerthsen, Crenshaw and Bethe* are combined. It may be 
seen that a good smooth curve can be drawn through them, so that the 
ratio S8.P.(H,)/S.P.(air) is well defined at all energies. We propose to 
ed 
* Livingston and Bethe, 1937, Rev. Mod. Phys., 9, 272. 
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use this curve to obtain S.P.(D,0)/S.P.(H,), and hence the energy loss in 
ice, in the following manner : 

(1) We read off values of S.P.(H,)/S.P. (air) for various energies. 
(2) We compute 8.P.(H,O)/S.P. (air) from the relation 


S.P.(H,0) S8.P.(H,) | 1S.P. (0,) 


S.P. (air) .P. (air) | 2S.P. (air) 
Rau 
265: 2. (air) aplanls 
Fig. 6. 
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0:0) : i O10 bo ine) 
The figure 0-535 for oxygen is valid for alpha-particles near the end of 
their range *, and cannot be very seriously in error at lower energies. The 
smooth curve resulting from this is shown as A in fig. 7. 

(3) We plot Crenshaw’s own values of S.P.(H,0)/S.P.(air) also in 
fig. 7, and draw through them a line (B) parallel to the curve A of (2) above. 
This we take to be the true representation of S.P.(D,0)/8.P.(air). We 
consider this procedure is preferable to using Crenshaw’s individual 
values of S.P.(H,O)/S.P.(air) as they stand, because it takes into account 
the probable trend of the stopping power with energy. With this curve 
we can now convert the results of Gerthsen, Eckardt and Reusse into 
energy losses in D,O. (We may comment, parenthetically, that Crenshaw’s 


i ee 
* Rutherford, Chadwick and Ellis, 1930, Radiations from Radioactive 
Substances (Cambridge: University Press), p. 97. 


550 A. P. French and F. G. P. Seidl on the 
Fig. 7. 


Stopping Power of H,O 


o8 


0-6 


Crenshow 
° = Bethe 


04 


02 


Proton energy (kev) 


10 100 1000 10000, 
' 


Fig. 8. 


Stopping Power ratios 
Aad Be be 


Proton energy (kev) 
loo 120 


80 


Energy Loss of Slow Deuterons in Heavy Ice 551 


values of the stopping power of D,O do seem surprisingly low, and cannot 
readily be joined to the computed curve at high energies. This, however, 
does not entitle us to ignore his results.) 

(4) Using the curve of fig. 6 and curve B of fig. 7, we find the ratio 
S.P.(D,0)/S.P. (air) 
S.P.(H,)/S.P. (air) 

The resulting curve is shown in fig. 8. For completeness our final curve 
of 8.P.(D,0)/S.P.(air) is also shown on this same graph. 


(5) With the aid of our curve from (4), we can now convert Crenshaw’s 
values (fig. 5) of energy losses in hydrogen into energy losses in D,O. 

The above procedure may seem unduly elaborate, but we consider that 
it extracts a maximum of significant data from Crenshaw’s results. To 
use the measurements on water vapour alone would be to discard the 
largest and most carefully conducted section of his work. 


—8.P.(D,0)/S.P.(H,). 


§7. THe Enerey Loss 1n Ice. 


We now have five sets of measurements which we can adapt to our 
needs in obtaining a final version of the energy loss curve. They are : 


1. Gerthsen Protons in air. 20 to 60 keV. 

2. Eckardt and Reusse Protons in celluloid. 4+to 50 keV. 

3. Crenshaw Protons and deuterons in H, and D,. 
(Equivalent proton energy : 30 to 340 keV.) 

4. Crenshaw Protons in H, and D,. 60 to 340 keV. 

5. Crenshaw. Deuterons in H,O. 


(Equivalent proton energy : 50 to 170 keV.) 
Using the results of the previous section, we have converted the first 
tour into energy losses in D,O vapour. The fifth is brought into line by 
plotting the points as though they were obtained for protons of half 
the energy. Fig. 9 shows the collected values, as usual in keV. per cm. path 
at 1 mm. pressure. The broken line is our attempt to integrate the various 
data into a single curve for the energy loss of protons from zero energy 
up to 340 keV. The maximum deviation from any of the parent curves 
may be seen to be about 10 per cent, which we regard as tolerably 
satisfactory in view of the numerous assumptions that have gone into the 
analysis. 

It remains only to make the final conversion into energy losses in ice, 
on the basis of our initial assumption (d) (§1). This we have done, with 
the supposition that both Gerthsen and Crenshaw conducted their 
experiments at a laboratory temperature of about 15°C. Under these 
conditions, one cm.? of D,O vapour weighs 1-11 micrograms. Table IIT. 
below lists the values of (dE/dm) in keV. per microgram at representative 
points over the whole range of energies covered in this review. 
Fig. 10 presents the results as a smooth curve over the limited range 
10-120 keV. deuteron energy. (Note that Table III. goes up to 700 keV. 
equivalent deuteron energy.) 
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$8. ConcLupIne REMARKS, 


Of the probable accuracy of the results given in Table III. we can 
say little. One might hope, taking an optimistic view, that they were 
good to +10 per cent, but we are after all confronted with complete 
uncertainty concerning the relative stopping powers of water vapour 
and ice. It is therefore of interest to mention a somewhat different 
approach to the problem. 

The analysis presented in this paper was undertaken in order to deduce 
the cross-section of the D(d, p) reaction, as a function of deuteron 
energy EK, from measured values of thick target yields N(E). In the 
original account (Bretscher and the authors 1948) of the experiments, 
values of both yields and cross-sections are given. The former are 
considered to be fairly accurate, but the latter were, by an oversight, 
obtained without reference to Crenshaw’s work on energy losses, and the 
present analysis should lead to more significant values of cross-sections. 


TABLE ITI. 
ee ee a es 


dE/dm oe (dE/dm) 

JE, | E, keV keV Jun VE, | Eykev. | yoy) 
2 4 0-211 ul 121 0-608 
3 9 0-325 12 144 0-556 
4 16 0-451 13 169 0-506 
5 25 0-580 14 196 0-455 
6 36 0-700 15 295 0-407 
7 49 0-753 16 256 0-356 
8 64 0-748 17 289 0-300 
9 81 0-712 18 324 0-262 
10 100 0-660 18-7 350 0-230 


This view has recently acquired support through the publication of 
results by Sanders, Moffatt and Roaf (1950). These workers obtained 
the D+D cross-section at a given energy directly by employing thin 
targets. Thus one can in fact use their values of cross-section, combined 
with values of dN/dE from Bretscher, French and Seidl (1948) to deduce 
values of dE/dx. If we do this, we find almost perfect agreement with the 
curve of fig. 10 at the lowest energy, and a maximum divergence of about 
-15-20 per cent at the highest energy for which a comparison is possible. 
(The values of Table IIT. rise more steeply than they should in this region), 
Unfortunately the measurements of Sanders e¢ al. do not go to a high 
enough energy to permit a comparison over the whole range of Table IIT., 
but the extent of agreement where the comparison is possible conforms 
rather well to the estimated accuracy of our values of dE/dz. 
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When more complete data on thin target yields are available, these can 
clearly be used to establish a more reliable energy-loss curve than that 
of fig. 10. In the meantime, however, it is hoped that our present range- 
energy relation may be of value. It should also be applicable to tritons 
in D,O, subject to our initial assumption (a) that hydrogen nuclei with 
equal velocities lose energy at equal rates. 
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SUMMARY. 


The photodisintegration of the deuteron has been observed in photo- 
graphic plates exposed, while wet with heavy water, to the y-rays arising 
from the bombardment of fluorine with 1-4 MeV. protons. Two groups 
of protons, produced by the known y-ray lines at 6-15 MeV. and 7-01 MeV., 
have been clearly resolved. It is shown that y-ray energies of this order 
of magnitude can be measured with a probable error of 0-05 MeV., and 
y-ray lines with an energy difference of 0-6 MeV. can be resolved. 
Suggestions are made for a further improvement of the technique whicb 
should increase the resolving power of the method. 

The angular distribution of the protons produced by each of the two 
y-ray lines has been measured, but statistical uncertainties must be 
further reduced before any useful comparison with theory can be made. 


§1. INTRODUCTION. 


THE photo-disintegration of the deuteron into a neutron and a proton, 
one of the simplest of nuclear processes, has been the subject of many 
experimental and theoretical investigations. The importance of the study 
lies, first, in the fact that any theory must be able to give a detailed 
account of the cross section and angular distribution of the products of 
the interaction for photons of different quantum energy, and accurate 
observations of these characteristics should provide a stringent test of 
the validity of different theories; and, second, because the process makes 
it possible to determine y-ray energies and thus provides a method for 
the investigation of y-ray spectra. 

At first sight it would appear that the reaction is particularly susceptible 
to detailed and accurate investigation. The absorption of the y-ray 
leads to the production of only two product particles, a neutron and 
proton, which recoil from one another. Suppose the quantum energy 
and the direction of motion of the y-ray to be known: its momentum is 
then defined, and the observation of the velocity of the recoiling proton 
is sufficient for the complete solution of the dynamics of the disinte- 


gration. 


* Communicated by Professor C. F. Powell, F.R.S. 
2, PZ, 
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In practice, however, it has proved to be very difficult to make 
observations of sufficient precision. This is largely due to the very low 
cross section for the reaction, which is of the order of 10-?? cm.” for 
y-rays of energy less than 20 MeV. As a result, competing nuclear 
processes can give rise, in some of the experimental conditions employed, 
to large numbers of protons which may be confused with those due to 
the disintegration of the deuteron ; and in other experiments the proton 
tracks have to be distinguished against a dense background of electrons 
produced by Compton recoil or pair production. 

In spite of these difficulties, the cross-section has been recently measured 
with a sufficient accuracy to provide a test of the validity of the various 
theories. Measurements have been made at quantum energies up to 
2-76 MeV. by Bishop, Collie et al. (1950), and at a series of energies, 
ranging from 4:45 MeV. to 17-6 MeV. (2:2 MeV. to 15:4 MeV. above the 
threshold), by Wilkinson (1950). Theoretical calculations for these and 
higher eneriges have been made by Hansson and Hulthen (1949) and by 
Marshall and Guth (1949). 

Information about the angular distribution ‘of the emitted protons is, 
however, very limited. Rough determinations, at a number of energies, 
have shown that there are no serious divergences from the predictions 
of the various theories, but the observations are not sufficiently precise 
to provide a crucial test. For such experiments Phillips, Lawson and 
Kruger (1950) have employed a Wilson chamber containing heavy 
methane, CD,, whilst other observers have detected the disintegration 
protons, emitted from solid or gaseous targets containing deuterium, by 
observing their tracks in neighbouring photographic plates. Most 
observations, however, have been made with the ‘“target’’ deuterium 
atoms incorporated in sensitive photographic emulsions, a method first 
suggested by Powell (1940). 

The use of emulsions “‘loaded”’ with deuterium has the advantage that, 
as in similar experiments with a Wilson chamber, the point of origin of 
a proton track can usually be accurately determined. Its length can 
therefore be measured with precision, and the inaccuracies which arise 
in experiments where the protons originate in a separate target of finite 
thickness can be eliminated. Further, the photographic method has an 
advantage over the expansion chamber, for these particular experiments, 
in being continuously sensitive. Whilst there are difficulties in applying 
the method to determinations of the absolute values of the cross-section 
for the reaction, it is well suited for making studies of the angular 
distributions. 

In applying the photographic method to the present problem, the most 
important technical difficulty to be overcome is that of incorporating 
the deuterium in the plates in sufficient concentration and in the form of a 
suitable chemical compound. The maximum y-ray flux which can be 
allowed to pass through the plate is limited by the maximum background 
fog, due to secondary electrons, which can be tolerated if the proton 
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tracks are to be distinguished. With this flux, and if the deuterium 
content of the plate is low, the number of disintegration protons per unit 
area is small and the work of searching the plate arduous. On the other 
hand, if the amount of loading material is large, there is a reduction in the 
quality of the proton tracks which leads to a reduction in the precision of 
the measurements of range ; and the inaccuracies are accentuated by the 
increase in the shrinkage factor of the emulsion, which leads to larger 
errors in the determination of the angles of dip of tracks. 

Several methods of “loading” an emulsion with deuterium have been 
employed :. Gibson, Green and Livesey (1947) described experiments in 
which calcium nitrate, with about 6 per cent by weight of heavy water 
of crystallization, was incorporated in the emulsion, but this method has 
various disadvantages ; the most important of these are that the calcium 
nitrate adversely effects the recording properties of the emulsion, and that 
the plates must be kept in hermetically sealed boxes to prevent water 
from the atmosphere taking the place of the heavy water in the emulsion. 

Tests have been made with emulsions containing about 50 per cent by 
volume of hexa-deutero-diacetin, but although these appeared promising, 
their properties do not yet appear to have been described in the literature. 

The method of loading most frequently employed has been to soak 
the emulsions in heavy water and to expose them to the y-radiation while 
wet (Goldhaber 1948, Waffler and Younis 1949, Hough 1950). The 
most important advantage of this procedure is that it gives a plate 
containing a large and reproducible concentration of deuterium, so that 
the number of proton tracks per unit area is relatively large for a given 
amount of y-ray blackening ; but there are the corresponding disadvan- 
tages, mentioned above, associated with the greatly increased volume of 
the emulsion during exposure and its large shrinkage when processed. 
The effect on the quality of the tracks can now be overcome, at least in ~ 
part, by the use of more sensitive emulsions, but the difficulty of the 
inaccuracies introduced in the measurements of the angles of dip remains. 

In the present paper it is shown that modern emulsions soaked in heavy 
water can be successfully employed for the measurement of y-ray energies 
in the region of 6 MeV. with an accuracy of 0-05 MeV., and allow y-ray 
lines with a separation of ~ 0-5 MeV. to be resolved. Secondly, results 
are described of experiments on the angular distribution of the protons 
- produced in the disintegration of deuterons by 6-15 MeV. and 7:01 MeV. 
y-rays, and methods for the further improvement of the measurements 
are discussed. 


§ 2. IRRADIATION OF PLATES. 


The y-rays employed for the irradiation of the plates were produced 
by bombarding a thick target of calcium fluoride with 1-4 MeV. protons. 
The latter were accelerated by the 1-7 MeV. Van der Graaf generator, 
built for the University of Bergen from designs by O. Dahl and B. Trumpy, 
which came into operation in March 1950. The target was 5mm. in 
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diameter, and eight plates were located with their centres 10 cm. from it 
and in a plane perpendicular to the proton beam. Four plates were 
orientated (see fig. 1) so that the y-rays entered their surfaces at, or 
near, normal incidence, and the other four so that the mean angle of 
incidence was 15°. This disposition of the plates was based on the 
following considerations :— 

The most accurate measurements of range are those made on tracks 
with small angles of dip. Further, for determining y-ray energies, the 
most accurate results are obtained from the measurements on protons 
ejected at nearly 90° to the line of motion of the parent photon; for in 
such cases the momentum of the photon has little influence on the 
calculated energy of the neutron, which is nearly equal to that of the 
proton. For a given exposure the plates placed with their surfaces 
normal to the y-ray flux yield the greatest number of tracks fulfilling these 
conditions, and they are therefore most suitable for the determination 
of the energy of the y-radiation. 


Fig. 1. 


He 


| 


Arrangement of plates around target : plates 1, 3, 5, 7 received y-rays approxi- 
mately perpendicularly, plates 2, 4, 6, 8 obliquely ; plates 1, 2, 3, 4, were 
soaked in D,O, 5, 6, 7, 8 in H,O. 


On the other hand, the best conditions for the investigation of the 
angular distribution of the protons are provided by plates in which all 
angles of emission are observed with equal probability. If observations 
are confined to tracks with small angles of dip, these conditions are 
provided by allowing the y-rays to pass nearly tangentially through the 
emulsion. For such experiments an angle of approach of 15° has been 
chosen to avoid appreciable reduction of the intensity of the radiation 
in passing obliquely through the walls of the containers and through the 
emulsions. 

The plates were contained in boxes of thin brass securely bolted to a 
base-plate which was attached to the target-holder, so that the geometry 


Photodisintegration of the Deuteron 559 


was well defined; there were no unnecessarily large masses of metal 
present, and the y-rays before reaching the emulsion had to traverse not 
more than a few millimetres of brass. 

An hour before the exposure, an Ilford C2 plate, 4in.x2in., with 
emulsion 200 «4 thick, was placed in each box, and half the boxes were 
filled with D,O and half with H,O. After exposure, during which 
the total flux of protons to the target corresponded to the passage of 
18 millicoulombs, the plates were dried to prevent fading of the latent 
image, and later processed by the “‘ temperature development ”? method. 
This procedure was found to be more convenient than one in which 
processing commenced immediately after exposure. 

Direct measurements with a micrometer showed that the ratio of the 
thicknesses of the emulsion during exposure and after processing and 
drying was 8-5+0-5. 


§ 3. EXAMINATION OF PLATES. 


A total of 943 tracks has been measured in one of the plates soaked in 
D,O and exposed with the y-rays entering at 15°. With some tracks it 
was impossible to be certain of the direction of motion of the proton, and 
these were rejected. It was confirmed that these tracks constituted a 
random sample, exceptional only because they were nearly rectilinear, 
the increased scattering commonly observed at the end of the range being 
not apparent. Their rejection therefore produced no effect on the 
measured distributions in energy or angle of the protons. For each of the 
943 tracks the following quantities were measured :— 


(i.) the coordinates, x and y, of the point of origin of the proton in 
the plate ; 

(ii.) the length, a, of its projection on a plane parallel to the surface of 
the emulsion, in units of length 0-85 p ; 

(iii.) the difference in depth, h, in microns, of the beginning and end 
points of the track: this quantity, multiplied by 10, gives, in 
units of 0-85 y, the original difference in depth of these points at 
the time the track was formed, since the shrinkage factor is 
equal to 8-5 ; 

(iv.) the angle « (see fig. 2) between the direction of the projection of 
the track on the surface of the emulsion and a standard direction 
OX, usually chosen to correspond to one of the edges of the plate. 


From these observations it is possible to calculate 
(v.) the angle of dip, 5, at the time of formation of the track: this 
angle is the mean inclination of the track, at, the time of its 
formation, to the surface of the emulsion; 
(vi.) the true range of the proton in the wet emulsion, in units of 0°85 p ; 
(vii.) two angles 7 and « (see fig. 2) defining the direction of motion of 
the parent y-ray and analogous to « and 6: 7 and « are functions 
of x and y ; 
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(viii.) the angle 6 between the direction of motion of the incident y-ray 


and the projected proton. 6 was directly determined from «, 
5,7, and e by a method of stereographic projection. 


we 


Diagram to show how a, 6, n, « define the proton and y-ray directions. 
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The results of the measurements for each track were plotted in a polar 
diagram showing R as a function of 0, those for which 8% 45° being plotted 
separately from those for §<45°. In the most unfavourable conditions, 
corresponding to large values of 8, the errors in 6 were sometimes as great 
as 5°, and in R five divisions. For §<45° the probable errors in 6 and R 
were 1-5° and one division respectively. A typical result for 425 tracks is 
shown in fig. 3. 
§ 4. CALCULATIONS. 

The distribution in energy of the y-rays, and the angular distribution 
of the disintegration protons in a coordinate system—the C-system or 
centre of mass system—moving with the deuteron after the absorption of 
the y-ray, were determined from the observed values of R and 6 by the 
following method : 

(a) Calculations of the stopping power to be expected for an emulsion 
saturated with D,O, based on the range-energy relation for a normal 
emulsion (Rotblat 1951) ,and the data of Livingston and Bethe (1937), give, 
for the range-energy relation between 1-5 MeV. and 3 MeV., 

R=19-3 E1® (in microns), ) 
or R=22-7 EY® (in units of 0-85 .), | 
where the energy E is measured in MeV. 

(6) The dynamics of the photodisintegration of the deuteron give, for 
the energy, E,, of the proton in the centre of mass system—to a very close 
approximation—the relation 

HeeWetl--O102 eos 6) 2 5% . (2) 
where E, is the observed energy of a proton emitted at an angle @ in the 
laboratory system of coordinates, the L-system. 

Equations (1) and (2) together give for the range Rg, of a proton of 
energy EK, 


(1) 


Roe=R (1—0-168 cos 8), 

where R is the observed range. R, is also the expected observed range of 
a proton emitted at right angles to the direction of motion of the parent 
photon. 

Now, to a close approximation, E, is related to the y-ray energy EK, by 
the relation 

: Eo=4 (E,—2-202), oe ere SoG AO I(4) 
all energies being measured in MeV. It follows that, for any particular 
y-ray energy, E, and the corresponding range R, may be calculated. 

By means of the above relations a map was constructed on tracing 
paper showing lines corresponding to given values of EK, on the polar 
diagram R, @. These lines are defined by the relation 


R (1—0-168 cos @)=constant, 


and with a very small (calculable) error they are circles, of radius Ro, 
with centre at a distance of 0-168 R, from the origin. They were drawn at 
intervals of 0-2 MeV. in the values of E,,. 


562 W.M. Gibson, T. Grotdal, J. J. Orlin and B. Trumpy on the 


The mechanics of the reaction give, to an adequate approximation, 
tan 0=sin ¢/(cos 6+0-051,) where ¢ is the angle of emission of the proton 
in the C-system. From this relation lines of constant ¢ were drawn on 
the map at 10° intervals. 

It may be noted that this map, reproduced in fig. 4, is completely 
independent of the measurements; it is based on theoretical calculation and 
well-established experimental facts, and no arbitrary parameters are 
involved. 

§ 5. Enercy DISTRIBUTIONS. 

When the map was superimposed on the plots of the results, it was seen 
that the points were distributed in two clearly marked bands, corres- 
ponding to two photon energies very close to the accepted values for the 
y-rays arising from the bombardment of fluorine with protons. The fact 
that in each plot the mean line through the points forming the more 
abundant band is so nearly parallel to the line for E,=6-15 MeV. gives a 


Fig. 4. 
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Map for obtaining y-ray energy E,, and centre of mass angle ¢ from R and @, 


very satisfactory check of the correctness of the calculations. The 
marked departure of the distribution of the experimental points in fig. 3 
from circular symmetry is a direct demonstration of the momentum of the 
photon, and the agreement of the experimental distribution with that 
anticipated is a confirmation of the generally accepted relationship between 
momentum and energy for a particle of zero rest-mass. 

Fig. 5 shows three plots of the distribution of the values of the energy of 
the y-rays as deduced from the observations. They were obtained, for 
0-2 MeV. intervals in the value of E,, by counting the numbers of 
experimental points between successive lines of constant E,. Fig. 5 (qa) 
is based on tracks with 55 45°, and shows poor resolution because of the 
inaccuracies, resulting from the large shrinkage factor, in the determin- 
ations of the range of these tracks. Fig. 5(b) is based on tracks with 
6<45°, and gives clear evidence for the presence of two unresolved groups 
of different intensity. Fig. 5(c), which shows two completely separated 
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groups with mean energies 6-15 +-0-05 MeV. and 7-06 -L0-07 MeV., has been 
deduced from the tracks for which § <45° and 60°<¢<120°. Curves (bd) 
and (c) presumably differ because the direction of motion of a small 
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intervals of E,, plotted against represent observed numbers of tracks 
E, ; (a) is based on tracks with with 8<45°, dots the numbers corrected 
angle of dip 45°, (b) on tracks for the effect of the limitation on the dip, 
with 86<45°, and (c) on tracks full lines the distributions predicted by 
with §<45° and 60°<$<120°. the simple theory, and broken lines the 


distributions predicted by the theory of 
Marshall and Guth. (a)| is for the 6-15 
MeV. group, (5) for the 7-01 MeV. group, 
and (c) for the two groups together. 
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proportion of the protons has been wrongly identified ; for values of ¢ near 
to 90° such an error makes little difference to the calculated value of E,, 
but the resulting error is much greater if the value of ¢ is near 0° or 180°. 

The y-ray energies 6-15+0-05 MeV. and 7-06+0-07 MeV. are in good 
agreement with the values 6-15 MeV. and 7-01 MeV. which the energies 
measured by Rasmussen, Hornyak and Lauritsen (1950), together with the 
relative yields obtained by Bennett, Bonner, Mandeville and Watt (1946) 
for the various resonances in the disintegration of fluorine by protons, 
give for the mean energies of the y-rays produced when 1-4 MeV. protons 
fall on a thick target containing fluorine. The close agreement for the 
6-15 MeV. line shows that systematic errors are slight, and that the 
uncertainty arising from “ straggling’ of the proton ranges, and from 
inaccuracies of measurement and geometry, is not greater than our estimate 
of 0:05 MeV. Less accuracy is claimed for the figure 7-06 Me\.., since it 
is based on relatively few tracks. 

The degree of resolution of the two groups shown in fig. 5 (c) indicates 
that it would be possible to observe as separate peaks two y-ray lines of 
which the energies differed by as little as 0-6 MeV., and to infer the 
existence of a doublet even if the separation were considerably less than 
this value. 


§ 6. ANGULAR DISTRIBUTIONS. 


The lines of constant ¢ in the map shown in fig. 4 are used to obtain 
distributions of the observed tracks, corresponding to each of the two 
y-ray lines, according to the angle of emission of the proton in the centre of 
mass system. The observed distribution with 6<45° is shown in fig. 6, 
together with the corrected results when account is taken of geometrical 
factors involved in limiting the tracks to those with angles of dip less than 
45°. The total number after correction is slightly greater than the total 
number of tracks observed with all values of 5, a result which is to be 
expected in view of the difficulty of observing very steeply dipping tracks, 
and of their greater chance of escaping from the emulsion. It was 
calculated that the effect of escape on the angular distribution obtained 
from tracks with 6<45° would be negligible in comparison with the 
statistical uncertainties, so no correction has been applied for it. 

Fig. 6 shows the expected angular distributions in the centre of mass 
system if the number of protons emitted per unit solid angle, at angle 4, 
were proportional to 

(a) sin? d; 

(b) a +sin® 4 (1428 cos 4), where f is the ratio of the velocity of the 
proton to that of light. (a) would result from a simple electric dipole 
interaction of the photon with the deuteron, and (6) if there were inter- 
ference between electric dipole and quadrupole terms (Marshall and Guth 
1949). The constant ‘‘a” is related to the nature of the neutron-proton 
interaction, but it has a magnitude of only about 0-02 at this y-ray energy, 
and makes very little difference to the shape of the curve. 
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Tt will be seen that our observed angular distributions are consistent 
with both (a) and (b) ; detailed examination of the results suggests that an 
increase of about a factor 4 in the number of measured tracks might 
reduce the statistical errors sufficiently to allow it to be shown whether or 
not this interference effect takes place. 


CONCLUSION. 


The results given above show that photographic plates soaked in heavy 
water are capable of giving useful information when used for measuring 
both y-rays energies and the angular distributions of protons from the 
photodisintegration of deuterons. It would therefore be valuable to increase 
the ease and accuracy with which the measurements can be made, and to 
this end we are studying the use of plates exposed at right angles to the 
incident y-rays and of more sensitive emulsions ; the former would make it 
possible to observe a given number of tracks suitable for the accurate 
determination of y-ray energy with less effort and calculation, and the 
latter would make the tracks themselves easier to find and would give more 
reliable measurements. This would probably more than offset any 
increase in the background due to electrons. Meanwhile, the statistical 
accuracy of the determinations of the angular distributions is being 
increased by the measurement of more tracks in plates exposed with the 
y-rays entering obliquely. 
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LIX. CORRESPONDENCE. 


Alpha-particles from the Proton Bombardment of Oxygen-18. 


By J. SEED, 
Cavendish Laboratory, Cambridge*. ( 


[Received February 28, 1951.] 


Tue yield of alpha-particles from the reaction 18O(px)*N has been studied 
using protons of energies 490 to 960 keV. having an energy homogeneity 
of 0-5 per cent. 

The alpha-particles emitted at an angle of 120° with the incident 
proton beam were analysed magnetically, and were detected by a thin 
zine sulphide screen and E.M.I. photomultiplier. 

The oxygen targets were in the form of thin layers of oxide on copper 
discs, and were prepared as follows: a polished copper target button 
was heated to redness in vacuo by means of an induction heater, and 
about 0:2 c.c. (at N.T.P.) of oxygen gas enriched about 30 times in 18O 
was allowed to enter, so that the surface of the disc was thinly oxidized. 
The oxide film was hardly visible, and from the observed width of the 
(px) group and known analyser resolution (4 per cent in energy), the 
target thickness was estimated to be about 5 keV. for 600 keV. protons. 

The yield of alpha-particles from these copper oxide targets was studied 
as a function of incident proton energy by observing the maximum 
intensity of the alpha-particle group obtained in the analyser: this is 
justified because the resolution width of the analyser was greater than 
the natural width of the alpha-particle group, calculated from target 
thickness and beam homogeneity. ; 

The excitation function, given in the figure, shows resonances at 
640 and 850 keV. The former resonance is close to that found at 680 keV. 
by Mileikowsky and Pauli (1950), who studied the alpha-particles from 
the 18O(pa) reaction at an angle of 90° with the incident protons. The 
difference is larger than can be accounted for by error in the proton 
energy in the present experiment, for which the voltage scale of the high 
tension equipment was calibrated directly at 441-4 and 873-5 keV. by 
observation of the thick target y-ray excitation functions of the reactions 
7Li(py)&Be (Fowler and Lauritsen 1949), and 19F(p, «y)!®O (Chao et al. 
1950). 

The magnetic analyser used in the present experiment was calibrated 
by observations on the alpha-particle and proton groups from the reactions 
*Be(dx)?Li, °Be(dx)7Li* (Q=7-150 --0-008 MeV. and Q—6-668 +0-009 MeV. 
respectively, Buechner and Strait 1949 and 1950 a) and #C(dp)8C 
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(Q=2-716+0-005), Buechner and Strait 1950 b). From measurement 
of the energy of the alpha-particle group from 18O(px) at the 850 keV. 
resonance, a Q value of 3:96-0:04 MeV. was obtained for the transition 
to the ©N ground state, in confirmation of the values reported in the 
literature (Burcham and Smith 1939, Freeman 1950). 

Taking this value for the energy release, the resonances at proton 
energies of 640 and 850 keV. would correspond to levels*at 8-59 and 
8-80 MeV. in the 1°F nucleus. No other reaction has yet given evidence 
of levels in 1°F in this region, though the reaction '8O(pn) shows closely 
spaced levels at about 10 MeV. excitation: by analogy with similar 
nuclei (1B, *N) having one particle less than an alpha-particle core, one 
would expect the 19F nucleus to have a fairly large number of low-lying 
excited states. 
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Excitation function for the reaction 48O(pa)!5N taken at 120°, and showing 
, the position of the 873-5 keV. resonance in the reaction °F (py). 


- A search for an alpha-particle group corresponding to transitions to 
excited states of “N was made using protons of energy 685, 837, and 
957 keV. and of homogeneity 30 keV. No group of intensity greater 
than one-tenth of the main group at alpba-energies greater than 1-2 MeV. 
was found. 
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The Isomeric State of Rak. 


By N. Featuer, F.R.5., 
The University, Edinburgh. 


| Received March 16, 1951.] 


Tue published results of Neumann, Howland and Perlman (1950) already 
make it very probable that *i{3RaE possesses an isomeric state which 
decays predominantly by «-disintegration, with disintegration energy 
5-12+0-05 MeV., to produce *{°Tl (4:2 m.). These authors concluded 
that the lifetime for «-emission from the isomeric state is greater than 
25 years. They also indicated some of the difficulties in accounting for 
the facts, which they established, that RaF is not produced to the extent 
of more than 1 in 2000 disintegrations of RaK* (either directly by 
B-disintegration or indirectly through isomeric transition to RaE followed 
by f-disintegration), and that RaD is not produced to the extent of more 
than about 1 in 200 disintegrations (by electron capture). On the 
assumption that the «-disintegration of RaE* is to the ground state of 
206T1—which leads to the smallest permissible value for the excitation 
energy of the isomeric state—they suggested that the chief difficulty lies 
in explaining the slowness of the isomeric transition RaE* — Rak. 
Even if this transition is of the fifth order they concluded that the 
excitation energy must lie very close to the lower limit of the experi- 
mentally deduced value of 0-18+0-11 MeV. Moreover, in this connection 
it should not be forgotten that the nucleus *{{RaE is already known to 
possess at least three non-isomeric states in the range of excitation 
energies from 0 to 0-05 MeV. (Feather 1949). 

Accepting a lifetime of the order of 100 years for x-emission, it is clear 
that the difficulty in explaining the slowness of the isomeric transition 
Rak* — Rak is more acute than that regarding the competing f-disinte- 
gration Rak*— RaF. It being generally accepted that the normal 
B-disintegration of Rak is characterized by the spin change 2->0, reference 
to the revised Sargent diagram of Feather and Richardson (1948) indicates 
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that £-disintegration from the isomeric state need not be characterized 
by a spin change of more than 5-—>0 in order to explain its non-obser- 
vation. This would only require the change 52 (with or without 
change of parity) for the isomeric transition, 

There are, however, some difficulties regarding the important assumption 
that the «-disintegration of RaE* is to the ground state of oleate 
spin of this B-active state is likely to be 0 or 1, for the 6-disintegration .of 
“Sr Ll is clearly allowed (Feather and Richardson 1948) and *88Pb is 
generally regarded as having spin 0. On the above assumption, then, 
and accepting the conclusions which follow from it, we should expect 
the «-disintegration *{gRaE*— *8°T] to be considerably less favoured than 
the parallel disintegration from the ground state Rak. The latter 
disintegration would be characterized by the spin change 2— 0(1), the 
former by the change >6—0(1). Plainly contrary to this conclusion 
is the evidence from the Geiger—Nuttall diagram for Z=83 (Broda and 
Feather 1947). On the basis of this diagram the «-disintegration of RaK* 
(r ~ 100 y.) appears more favoured, by a factor possibly as great as 10 
(rather than less favoured by a factor of the order of 103) than the 
x-disintegration of RaK, due allowance being made. for the difference in 
energy. It would appear that there are three possibilities of resolving 
this seeming contradiction : 

(i) the lifetime of RaK* for x-emission may in fact be of the order of 

10° to 10° y.* 
(ii) the £-disintegration of RaE* may be to an isomeric state of 72°T1; 
(iii) the «-disintegration of RaE* may be abnormal in the sense that 
it violates the general regularities exhibited in the Geiger—Nuttall 
diagram. 

No doubt direct evidence concerning (i) will not long continue lacking, 
but, judging the present position it would appear surprising that the 
activity of RaK* has been detected at all, if its lifetime is really as long as 
105y. Also, the longer the lifetime for «-emission, the smaller the experi- 
mental upper limit to the transition probability for the isomeric transition 
Rak* — RaE and the larger the spin which has to be assigned to the 
isomeric state in order to explain its properties. Possibility (i), in fact, 
does not alone provide a satisfactory basis of resolution. 

Possibility (ii) introduces similar difficulties. Since its acceptance 
leads inevitably to the conclusion that the excitation energy of the 
isomeric state RaK* is greater than previously supposed, it likewise 
makes more difficult the problem of explaining the relative slowness of all 
processes, except «-disintegration, which lead to the de-excitation of this 
state. Again, an isomeric state of *{/Tl is hardly likely to have escaped. 
detection in the experiments of Fajans and Voigt (1940) and Neumann 
et al. (1950). 

Concerning possibility (iii), this would be realized formally if the 
effective nuclear radius of the isomeric state RaK* were greater than that 
of the ground state RaE. Physically, however, the smallness of the 
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x-disintegration constants which characterize all ground-state nuclei 
having Z—83 (when these are compared with the «- -disintegration constants 
of an for which Z>84) is most probably connected not so much 
with a real decrease in nuclear radius as Z decreases from 84 to 83 as with 
the difficulty of formation of the nascent «-particle when only one of its 
constituent protons is available in a loosely bound state (‘ outside” the 
closed shell of 82 protons). If possibility (iii) were to be interpreted in 
this way it would obviously imply that one of the properties of Rak*, the 
isomeric state of the nucleus containing (126+1) neutrons and (82-1) 
protons, is that the intranuclear formation of an «-particle is less difficult 


% 


in this state than it is in the ground stat of the same system—or in any - 


of the known non-isomeric states of any nucleus of proton number 83 
and neutron number 128 or greater. This is a possibility of a type not 
hitherto considered in theories of «-disintegratiou. 
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An Example of the (n, p; m~) Reaction in the Photographic Emulsion. 


By 8. J. Gotpsack and N. Page, 
University, Manchester *. 


| Received March 9, 1951.] 


[Plate XX.] 


We have recently found the event shown in Pi. XX. in a 400 micron thick 
Ilford G5 Nuclear Research Emulsion, exposed vertically for 21 days at 
an altitude of 2860 m. in a magnetic field of 34,000 gauss (Dilworth et al. 
1950). The event is interpreted as the production of a 7~-meson (track A) 
in the collision of a high-energy, cosmic-ray neutron with a proton in 
the emulsion, according to the following scheme : 


Doct De Deb tenn ee aa et eal an Te 


The technical details of the event are as follows. Track A is almost 
continuous throughout its length and shows marked scattering especially 
towards the right-hand end of the track where there are two sharp 
deflections. The grain density of the track increases noticeably from the 
ns ee See 
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centre of the star up to the second sharp deflection, but at this point there 
is a decrease in grain density and the track becomes quite straight (D). 
The grain density and scattering of track A are consistent with its having 
been produced by a slow meson (presumably a 7 -meson) moving to the 
right and producing a single-prong star (D) at the end of its range. The 
small scattering and high grain density of track D exclude the possibility 
of its being a «*+-meson arising from the decay’of a positive 7-meson. 
The range of the meson is 370 corresponding to a kinetic energy of 
3-6 MeV. 

The meson is created in a star which has two lightly-ionizing particles 
(B and C) but no evaporation or recoil particles. The track B is only 
180 « long and is therefore too short to be identified, but its grain density 
is near minimum. Minimum ionization in this plate corresponds to 
34:2+0-8 grains/100: track B has 37-:0+6-0 grains/100y. Track C, 
which is 1-1 cm. long, has a grain density of 46-8 -+0-7 grains/100 and is 
therefore slightly above minimum. The small angle multiple scattering 
of this track has been measured by the method of Fowler (1950), and the 
arithmetic mean angle of scattering is 0-048 +0-006°/100 1 corresponding 
to a value of pB=700+85 MeV./c. 

The grain density of a proton showing this scattering should be 
1:5+0:15 times minimum. The particle is therefore identified as a 
proton, the kinetic energy being 400-+50 MeV., and the momentum 
990 +120 MeV./c. This identification is confirmed by the sign of the 
magnetic deflection. The total deviation between the ends of the track 
is 1-23°. To this must be added a small correction for the distortion 
of the emulsion : measured by the method of Cosyns and Vanderhaeghe 
(1950), this proves to be +0-05°. Assuming that the particle is moving 
downwards from the star, it is positive with a probability of 97 per cent. 

It is clear that the equation (1) is the simplest interpretation of the 
event which is consistent with conservation of charge, the third particle 
being a knock-on proton, initially at rest in the emulsion. 

If the particle B is a proton, its energy must be at least 650 MeV. and 
hence the kinetic energy of the incident neutron which produced the star 
must have been at least 1200 MeV. 

Finally it is of interest to review the evidence for possible simple 
schemes for single meson production in nucleon-nucleon interactions. 
These are collected in Table I. | 

It is seen from the Table that at least six such reactions are possible 
all of which might be observed in the Ilford G5 emulsion. Events | and 2 
can occur in collisions with hydrogen nuclei in the emulsion which, it is 
estimated, account for 4 per cent of the total geometrical cross section 
of the emulsion. All six events could occur in collisions with nucleons 
in the nuclei of heavier atoms in the emulsion, but in such cases an 
evaporation star or a nuclear recoil fragment would usually be observed. 
It is possible, however, that the two interacting nucleons might leave 
the nucleus without imparting an appreciable amount of energy to it. 
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Event 2 has been observed in the cloud chamber (Armenteros et al. 
1951 a); two examples of type 3 occurring in Ilford G5 emulsion have been 
reported by Cosyns et al. (1949), and a similar event was produced in 
Ilford C2 emulsion using artificially produced neutrons from the Berkeley 
cyclotron (Smith et al. 1950). Armenteros et al. (1951 b) have pointed 
out that the event described by Hopper and Biswas (1950) and inter- 
preted by them as an example of the decay of a neutral V particle similar 
to those found by Rochester and Butler (1947) and by Seriff et al. (1950) 
may be another example of this process. 

At present it appears impossible to distinguish, in the emulsion, between 
high energy events of the types 3 and 4 b and the decay of V particles. 


TABLE I. 


Types of simple star involving the production of one meson. 


No. of 
Products | prongs per References 
star 


Incident} Target 
Particle | Particle 


Present letter 


Neutron| Proton | p+p+a- 3 
Neutron] Proton | n+n-+-a* if 
3 


Armenteros ef al. 
5 195la 

Neutron| Neutron] »+p+a- y Cosyns, et al .1949 
Smith, et al. 1950 


Proton | Proton | n+p-+at 


Proton | Neutron| p+p+-a- 
Proton | Neutron} n--n-+-at 
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